- CLN-list - ginac.de

Re: [CLN-list] how to locate arrays in CLN library
by Richard B. Kreckel 19 Jun '09

19 Jun '09

Hi! Breta, please direct qustions regarding CLN to the mailing list. Thanks. sopik(a)fzu.cz wrote: > first of all, I'd like to thank you for being interested. You're right, it > wasn't very clever to ask without giving any example code. Here it is > > #include <cln/complex.h> > #include <cln/complex_io.h> > #include <cln/float.h> > #include <iostream> > > using namespace std; > using namespace cln; > float_format_t default_float_format = float_format(40); > > int main(void) > { > int i, N = 10; > cl_N *A; > > A = (cl_N*)malloc(sizeof(cl_N)*N); > for(i=0; i<N; i++) { > *(A+i) = complex(cl_float(i),cl_float(2)); > } > for(i=0; i<N; i++) { > cout << *(A+i) << endl; > } > cout << "See, it works!" << endl; > } > > The understandable error message is: > Unauthorized access to memory (SIGSEGV) > > However I don't know, how to do things properly. > Thank you for any advice. Remember that malloc gives you a chunk of *uninitialized* memory. And inside the first loop you are assigning one cl_N to another one, but the first one does not exist yet. That results in undefined behavior. You could either read up on the C++ feature "placement new" or, better, simply use C++ standard containers instead of malloc'ed memory. Try pushing your cl_N into a std::vector<cln::cl_N>. This way, you don't have to worry about freeing your memory later. Hope this helps -rbk. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

2 1

Re: [CLN-list] how to locate arrays in CLN library
by sopik＠fzu.cz 19 Jun '09

19 Jun '09

Hi, first of all, I'd like to thank you for being interested. You're right, it wasn't very clever to ask without giving any example code. Here it is #include <cln/complex.h> #include <cln/complex_io.h> #include <cln/float.h> #include <iostream> using namespace std; using namespace cln; float_format_t default_float_format = float_format(40); int main(void) { int i, N = 10; cl_N *A; A = (cl_N*)malloc(sizeof(cl_N)*N); for(i=0; i<N; i++) { *(A+i) = complex(cl_float(i),cl_float(2)); } for(i=0; i<N; i++) { cout << *(A+i) << endl; } cout << "See, it works!" << endl; } The understandable error message is: Unauthorized access to memory (SIGSEGV) However I don't know, how to do things properly. Thank you for any advice. Breta > Hi, > > sopik(a)fzu.cz wrote: >> I'd like to to work with fields of complex numbers of high precision using >> the CLN library. I can't successfuly allocate the memory using memalloc however. I tired few possibilities, but nothing worked. There is no notice >> about this topic in manual or in mail discussion, so it seems to be straightforward. But I have no idea how to do it. > > I'm not sure what you are trying to do, but working with complex numbers is straightforward. This is C++ and you should rely on the power of constructors/destructors and, if necessary, use new/delete for memory management, instead of malloc/free. > > Here is a working mini-example: > > #include <cln/cln.h> > #include <iostream> > using namespace std; > using namespace cln; > > int main() > { > cl_N x = -4; > cl_N sqrt_x = sqrt( x ); > if ( sqrt_x == complex( 0, 2 ) ) { > cout << "See, it works!" << endl; > } > } > > If this doesn't help, please post some example code. > > Best wishes > -rbk. > -- > Richard B. Kreckel > <http://www.ginac.de/~kreckel/> >

1 0

how to locate arrays in CLN library
by sopik＠fzu.cz 18 Jun '09

18 Jun '09

Hi, I'd like to to work with fields of complex numbers of high precision using the CLN library. I can't successfuly allocate the memory using memalloc however. I tired few possibilities, but nothing worked. There is no notice about this topic in manual or in mail discussion, so it seems to be straightforward. But I have no idea how to do it. Thanks for a help Breta Sopik

2 1

Re: [CLN-list] msvc support patch
by Bruno Haible 10 May '09

10 May '09

Michael Goffioul wrote: > This link will hopefully clarify this. > > http://www.gnu.org/licenses/gpl-faq.html#WindowsRuntimeAndGPL This FAQ does not make the distinction between whether the installer contains the run-time libraries or not. Based on the verbiage in the GPL (explained in the other mail) I therefore think that this FAQ assumes that the installer does not contain the run-time libraries. > there's no point in discussing the rest. There is, because you or the Octave people have other ways to build an installer that fulfils the letter of the GPL: - Build with an older version of MSVC++, whose runtime libraries are contained in the minimum version of Windows that your package supports, or - Distribute the the MSVC++ runtime libraries in a separate installer (then it's mere "aggregation" in the terms of the GPL), or - Arrange for the installer to pop up a browser window to the Microsoft VC++ redistributables download page. The last two options are still not following the spirit of the GPL (because they encourage the installation of closed-source libraries), but legally not a GPL violation. In the second option, there's also the Microsoft license that you have to look at. I find it funny that the license in the "Microsoft Visual C++ 2008 Redistributables" (vcredist_x86.exe) expressly forbids redistributing its contents. Bruno

2 4

Re: [CLN-list] [GiNaC-list] msvc support patch (Was: use of STL iterators in GiNaC)
by Michael Goffioul 21 Apr '09

21 Apr '09

On Sun, Apr 19, 2009 at 9:53 PM, Richard B. Kreckel <kreckel(a)ginac.de> wrote: > Hi! > > Michael Goffioul wrote: >>> >>> (How) Did you manage to compile CLN with msvc? >> >> Magic. >> >> VC++ per-se is quite (not fully) standard compliant. Most of the >> problems are in the autotool based build framework. With a bit >> of coding, you can make VC++ look like GCC. > > Interesting. Could you write up a summary of what you did and report it on > <cln-list(a)ginac.de>? This way, others will have a chance to pick up on your > work. CLN/GiNaC compilation with VC++ is part of the making of Windows package for octave. All compilation scripts and patches are available under SVN on octave-forge http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/admin/Wi… (run_compilation.sh and libs/) The main problem when compiling packages with VC++ is not the compiler itself, which is pretty C/C++ compliant, but the build toolchain. Most packages in the UNIX world use autoconf/automake/libtool toolchain and those tools do not support VC++ very well (at all). But you can work around that problem by using wrappers around VC++ compiler to make it look like GCC, by translating command flags. You can use such wrappers to build tools equivalent to gcc, g++, ar, ranlib... but based on VC++. You can also provide some missing headers like stdint.h, inttypes.h or unistd.h (you can find some on the web), to further ease porting packages. Once you have those tools, it becomes easier to compile packages using the normal "./configure && make" way, using the MSYS shell. But you still need to take care of some specific issues: 1) libtool: using some post-processing magic, it's possible to make libtool to support VC++ (through the gcc-like wrappers) 2) symbol export: VC++ does not export symbols automatically, so you need to take of it; either by using a .def file (some packages provide their own), or by using dllexport/dllimport qualifiers (again some packages support that as well), or by parsing yourself the object files to generate a .def file automatically (this works well especially for plain C code); of course, when you build a static library (which is what I do or CLN/GiNaC), you don't have to care about that 3) export of data: even with the tricks of point 2), you still need to decorate exported data with dllimport, as VC++ does not have the auto-import feature of MinGW; for packages that do not provide this decoration, you have to write your own patch. For the specific case of CLN, there were a couple of issues I also had to take care of: - struct/class, const/non-const: at some places, CLN foward declares objects as being "class", while there are actually "struct"; VC++ generates different mangled names for struct and class, so you end up with unresolved references. The same happen for const/non-const (forward declare as non-const, but define as const). - assembly code: here I cheated a little bit and used gcc to compile the .S file; this required some hacks in the Makefile - some namespace issues: apparently, declaring a function inside the body of another method, which is parted of the cln namespace, does not make VC++ to tag the declared function with the cln namespace; for instance in number.h, I had to move the "extern" statements of CL_DEFINE_LONG_CONSTRUCTOR out of the method bodies. Bye, Michael.

4 6

CLN repository, the fork, and all that
by Richard B. Kreckel 04 Feb '09

04 Feb '09

Hello, Here is some information about the immediate future. I've traded the CVS repository for a git repository. The CLN web page informs about how to obtain the current sources and how to build them. Actually, that is a clone of the CVS-to-git conversion done by Alexei Sheplyakov, who announced it as a fork on this list a while ago. Alexei has explained that he intends to continue to contribute to CLN. Technically, he will do so by sending requests to pull from his repository. After all, that is how git is designed to work. In this spirit, all his patches to make CLN less dependent from GCC are already in the current repository. Kudos to Alexei for making this possible. Everyone is invited to give the current version a test. Please report any problems you may find. Please send reports to this list. Better yet, send patches! ;-) My immediate plan is to prepare a release at around the time GCC-4.4 is released. Due to some minor header problems, CLN-1.2.2 won't compile on that new compiler version. Since Alexei's work made changes to the ABI, the new release will have to increment the soname version. Cheers -richy. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

1 0

modular polynomials
by Jens Vollinga 30 Sep '08

30 Sep '08

Hi, using cl_UP_MI I get very strange effects that look like problems with the memory management (the contents of one polynomial is suddenly vanishing while operating on another for example). So far, I figured out that having the leading coefficient set to zero via set_coeff() causes trouble even if I call finalize(). The manual is kind of sparse on that topic so I am not quite sure whether this is a bug or a feature. But even when I try to avoid that case, call finalize everywhere, and only do assignments (hope they are safe ...), then polynomials aside from the ones involved in the operations change. I will investigate this case further but it would be helpful to know whether the cl_UP_MI has problems/bugs? Regards, Jens

4 10

FYI: CLN fork
by Alexei Sheplyakov 27 Aug '08

27 Aug '08

Hello! CLN is a great library, but it has a number of problems. 1. CLN can be built only with GNU C++ compiler. Even worse -- some bits of CLN (in particular, include/cln/modules.h) work only with specific versions of GCC, so CLN can break at any moment (this actually happens with every new [major] release of GCC). 2. The code is a bit obfuscated. 3. The build scripts are extremely user- and hacker-hostile. Since CLN developers keep [silently] refusing patches addressing these issues, I've decided to fork CLN. The code can be obtained from git repository: git://ffmssmsc.jinr.ru:443/varg/cln.git To view the repository on-line point your browser to http://ffmssmsc.jinr.ru/~varg/git/gitweb.cgi?p=varg/cln.git Best regards, Alexei -- All science is either physics or stamp collecting.

2 1

floating point overflow
by David van der Spoel 03 Jun '08

03 Jun '08

Hi, I get the following message when trying to do a rather complex piece of code: terminate called after throwing an instance of 'cln::floating_point_overflow_exception' what(): floating point overflow. I am running on a Mac in 64/Intel bit mode as per the instructions in previous posts, with GMP. This is using cl_F with 100 decimals precision throughout. Any clue what causes this and how to prevent it? -- David van der Spoel, Ph.D., Professor of Biology Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University. Box 596, 75124 Uppsala, Sweden. Phone: +46184714205. Fax: +4618511755. spoel(a)xray.bmc.uu.se spoel(a)gromacs.org http://folding.bmc.uu.se

2 1

CLN 1.2.2 released
by Richard B. Kreckel 07 Apr '08

07 Apr '08

Hi! CLN 1.2.2 has been released. This release re-establishes binary compatibility with CLN-1.2.0 and fixes the ARM build, both inadvertently broken in the previous release. Cheers -richy. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

3 7