- CLN-list - ginac.de

We are moving...
by Richard B. Kreckel 08 Oct '05

08 Oct '05

Dear all, Due to ongoing technical problems at Mainz University (*), it was decided to move GiNaC's and CLN's web server / CVS server / mailing lists to a rented server. This should allow us to continue communicating via mailing lists, doing development and updating the web pages, things that have become difficult during the last 10 weeks or so. We are going to move the DNS entries ASAP. The hosted machine's IP-address is going to be 81.169.130.167. We'll certainly see some outage while the transfer occurs, but hopefully not for long time. Meanwhile, you can all make sure the address you use as mailing list addresses is ...(a)ginac.de, not ...(a)thep.physik.uni-mainz.de or ...(a)zino.physik.uni-mainz.de because these will cease to work. We apologize for the inconveninece -richy. (*) Case in point: right now the web server is dead. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

1 0

DLL for mingw?
by David Bateman 05 Oct '05

05 Oct '05

I'm trying to build a DLL of CLN for mingw. I absolutely must have a DLL as my use of CLN is with octave (www.octave.org) and with the lack of symlinks under windows if I use a static library, CLN with be linked in to all of the seperate octave functions, and the amount of disk used with be enormous... However, I have a problem of many many resolved symbols in CLN when linking the DLL. I imagine the same problem must exist under cygwin as well. A sample of the type of error message I'm getting are exam.o(.text+0x133): In function `Z41__static_initialization_and_destruction_0ii': exam.cc:50: undefined reference to `global constructors keyed to cl_module__cl_prin_globals__firstglobalfun' exam_I.o(.text+0x96):e:/octave/build/bld/cln-1.1.9/tests/: undefined reference to `cln::cl_I::cl_I(char const*)' exam_RA.o(.text+0x96):e:/octave/build/bld/cln-1.1.9/tests/: undefined reference to `cln::cl_RA::cl_RA(char const*)' etc... Though I have hundreds of such unresolved symbols... Is there something basic that I'm missing? Or are these really unresolved and the problem only appears in the windows builds? Regards David -- David Bateman David.Bateman(a)motorola.com Motorola Labs - Paris +33 1 69 35 48 04 (Ph) Parc Les Algorithmes, Commune de St Aubin +33 1 69 35 77 01 (Fax) 91193 Gif-Sur-Yvette FRANCE The information contained in this communication has been classified as: [x] General Business Information [ ] Motorola Internal Use Only [ ] Motorola Confidential Proprietary

1 0

Types in cl_float
by Isidro Cachadiña Gutiérrez 12 Sep '05

12 Sep '05

Hi: I miss a cl_float(long double x, cl_F ) definition to make the code below to compile. The code compiles when A is defined as double. ----------------------------------------------------------------------------- #include <iostream> #include <cln/float.h> #include <cln/float_io.h> using namespace cln; using namespace std; int main(int argc, char **argv) { long double A; A=1; cl_F a; float_format_t precision; precision=float_format(100); a=cl_float(A,precision); cerr << a << endl; } --------------------------------------------------------------------------------------- -- *************************************************************************** * Dr. Isidro Cachadiña Gutiérrez * * Departamento de Física * * Facultad de Ciencias * * Universidad de Extremadura * * 06071 Badajoz ( SPAIN ) * * email: icacha(a)unex.es * * Teléfono: +34 924 289 300 Ext. 6826 Fax: +34 924 289 651 * *************************************************************************** * Usuario Linux: 8569 * * Para clave pública GnuPG: http://onsager.unex.es/firma.pub.asc * ***************************************************************************

1 0

CLN on MINGW
by David Robinow 23 Aug '05

23 Aug '05

The following patch allows building on Windows under MSYS. make check succeeds. --- /c/home/cln_cvs/cln/src/base/random/cl_random_from.cc 2004-06-23 17:11:21.000000000 -0400 +++ cl_random_from.cc 2005-08-23 11:31:13.938737700 -0400 @@ -50,6 +50,14 @@ return seed_lo + tmsbuf.tms_utime + tmsbuf.tms_stime; } +#elif defined(__MINGW32__) +#include <sys/timeb.h> +inline uint32 get_seed (void) +{ + struct timeb tm; + ftime(&tm); + return cln::highlow32(tm.time,tm.millitm * 1000); +} #endif #endif

1 0

Python Interface to Class Library for Numbers
by Ray Buvel 20 Aug '05

20 Aug '05

I have just released a Python interface for the Class Library for Numbers. For more information see http://calcrpnpy.sourceforge.net/clnum.html Ray Buvel

1 0

CLN on MinGW
by John Pye 12 Aug '05

12 Aug '05

Hi again Richy, Just wanted to ask if anyone's succeeded in building CLN on MinGW / MSYS ? Is this something that other people have shown an interest in? I get the following error when I try to ./configure && make: > g++ -g -O2 -I../include -I../include -I./base -Ibase -c > ./base/random/cl_random_from.cc -DDLL_EXPORT -DPIC -o > .libs/cl_random_from.o > ./base/random/cl_random_from.cc: *In constructor > `cln::random_state::random_state()':* > ./base/random/cl_random_from.cc:68: *error: `::get_seed' has not been > declared* > make[3]: *** [cl_random_from.lo] Error 1 The problem could be as simple as needing to alter the #if stuff like (defined(_WIN32) && defined(__GNUC__)) -- I think that MinGW may use slightly different macros eg WIN32 in place of _WIN32. Or it could be more complex, I know that MinGW doesn't implement glibc so you might need more complicated #ifdefs and provide some aliases to MFC library calls. My experience with MinGW is that it runs a lot more quickly -- both compiling and then running the compiled code. The Cywin POSIX layer seems to slow things down quite a bit in computationally intensive stuff. Cheers JP -- John Pye School of Mechanical and Manufacturing Engineering The University of New South Wales Sydney NSW 2052 Australia t +61 2 9385 5127 f +61 2 9663 1222 mailto:john.pye_AT_student_DOT_unsw.edu.au http://pye.dyndns.org/

2 2

Options for FP errors in CLN?
by John Pye 29 Jul '05

29 Jul '05

Hi there, I've been doing a lot of work with GiNaC recently but just now I think I have a problem which relates more specifically with CLN, so hopefully I'm asking this question in the right place. I have had a few problems now with CLN's error handling: it tends to just abort when something goes wrong. Is there a reason why it doesn't just throw exceptions that I can catch? >void cl_error_floating_point_nan (void) >{ > fprint(std::cerr, "floating point NaN occurred.\n"); > cl_abort(); >} > > I can see that I could rewrite the 'cl_abort' function so that it does throw exceptions, but I just wanted to point out that CLN is a big job to compile and it might be nice if it gives this sort of functionality out-of-the-box. I can't use the suggested method from the manual either, > The library does type checks, range checks, consistency checks at many > places. When one of these fails, the function |cl_abort()| is called. > Its default implementation is to perform an |exit(1)|, so you won't > have a core dump. But for debugging, it is best to set a breakpoint at > this function: > > >(gdb) break cl_abort > > > When this breakpoint is hit, look at the stack's backtrace: > > >(gdb) where > > because my pre-compiled CLN doesn't seem to have debug symbols: > [john@jdpipe nla-ginac]$ gdb pipesystem > GNU gdb Red Hat Linux (6.1post-1.20040607.43rh) ... > This GDB was configured as "i386-redhat-linux-gnu"...Using host > libthread_db library "/lib/tls/libthread_db.so.1". > > (gdb) break cl_abort > Function "cl_abort" not defined. > Make breakpoint pending on future shared library load? (y or [n]) y > Breakpoint 1 (cl_abort) pending. > (gdb) run > Starting program: /home/john/nla-ginac/pipesystem > Reading symbols from shared object read from target memory...done. > Loaded system supplied DSO at 0xb7f9f000 ... > floating point NaN occurred. > > Program exited with code 01. > (gdb) where > No stack. So, I'd like to suggest replacing cl_abort with either a SIGFPE or a standard exception of some sort. But please, if there are other issues at hand, or reasons why you don't do things that way, I'd be keen to know. Cheers JP

3 3

floating point error
by Elvis Dieguez 12 Jul '05

12 Jul '05

Hi, This a bit difficult to explain but let me try... The boost matrix library uses templates extensively in its code. I am using the library in the following way (for example) matrix<cl_F> G(3,3), T(3,3); G(0,0) = "-0.3e-5_5000"; G(0,1) = "0.5e-5_5000"; G(0,2) = "0.8e-5_5000"; ... etc When I try the following T(0,0) = G(0,0); for(int i = 0; i < 10; i++) T(0,0) = T(0,0)*G(0,0); cout << endl << T << endl; It works just fine, which is no surprise since I specified a 5000 precision. However, if I try the following T = prod(G,G); for(int i = 0; i < 10; i++) T = prod(T,G); cout << endl << T << endl; I get a floating point underflow error. I think the problem is the call to prod(*,*) creates a matrix where the cl_F matrix entries are defined with some smaller, default precision (probably equal to that of a long float from what I can tell). Is there a way to set the default precision for all created cln variables? If not, I will have to write my own matrix multiplication function... something I would not rather do if at all possible. I've looked at the cln documentation and it appears |float_format_t default_float_format| Global variable: the default float format used when converting rational numbers to floats. there is a way but I do not quite understand how one might change the default precision. Thank you for your help. Elvis Dieguez -- ...how one lives is so far distant from how one ought to live, that he who neglects what is done for what ought to be done, sooner effects his ruin than his preservation; for a man who wishes to act entirely up to his professions of virtue soon meets with what destroys him among so much that is evil. -- Nicolo Machiavelli, The Prince (1515)

2 1

converting cl_F to double?
by Leo Liberti 11 Jul '05

11 Jul '05

Dear CLN List, could anyone please tell me how to convert a cl_F to a double? I would like to write something like cl_F var = 1.0; double d = var; but the error is "`class cln::cl_F' used where a floating point value was expected". Thanks, Leo ___________________________________ Yahoo! Mail: gratis 1GB per i messaggi e allegati da 10MB http://mail.yahoo.it

2 1

Floating point underflow.
by Elvis Dieguez 01 Jul '05

01 Jul '05

Hi, I wrote some code using the Boost matrix library and CLN. But I get 'floating point underflow' errors every time. I thought by using CLN, I would not have problems with very small and very large numbers. Is there anything I can do.. I really need this to work! Sincerely, Elvis Dieguez

2 1