- CLN-list - ginac.de

[cppcheck] Member variable 'gcd_result::gcd_degree' is not initialised in the constructor.
by Martin Ettl 23 Jan '11

23 Jan '11

Hi, during a check with the static code analysis tool cppcheck of the current git head of cln, the tool discoverd a not initialized member variable of class gcd_result::gcd_degree. Please refer the attached patch that fixes the issue. Best regards from the cppcheck team Ettl Martin

2 1

exponants of floating points (or do I need it ?)
by Franck Sabatié 14 Jan '11

14 Jan '11

Hello, I have a hard time seeing how to best implement the following : I have a complicated equation of 4 variables that I need evaluated with a certain number of digits in precision (30 or so) : dummy = (15/4)*(1/(pow(xi,5)*(1+i-k)*(2+i-k)*(3+i-k)*(4+i-k)*(5+i-k)))*( ( (3*pow(3+i-k,2)-(1+i-k)*(5+i-k))*pow(xi*xi-xb,2) -2*(2+i-k)*(4+i-k)*(1-xi*xi)*(xb*xb-xi*xi) ) * ( pow((xb+xi)/(1+xi),3+i-k)- pow((xb-xi)/(1-xi),3+i-k) ) + 6*(3+i-k)*xi*(1-xb)*(xi*xi-xb)*(pow((xb+xi)/(1+xi),3+i-k)+ pow((xb-xi)/(1-xi),3+i-k)) ); where the variables are i, k, xb and xi. Although most of the equation is relatively easy to implement in CLN, I have a hard time figuring out how to deal with the "pow". There appears to be no such exponant function for floats ? Also, it is still unclear to me how to enter this equation painlessly into my code. Do I have to put cl_float(1,30) when I use 1 for instance ? (same for other numbers). Should I evaluate this in floating precision or with "real numbers" ? And finally how do I evaluate expressions such as (something)^(something else) with the (something else) NOT being an integer. Sorry about the very basic questions but most of the demo examples deal with integers and the few that do not, have very little information in them. Thank you for your time, Franck

2 1

Power function
by Joshua Friedman 14 Jan '11

14 Jan '11

I think you need to compute exp(k ln(x) ) for x^k Sent from my iPad On Jan 14, 2011, at 6:00 AM, cln-list-request(a)ginac.de wrote: > Send CLN-list mailing list submissions to > cln-list(a)ginac.de > > To subscribe or unsubscribe via the World Wide Web, visit > https://www.cebix.net/mailman/listinfo/cln-list > or, via email, send a message with subject or body 'help' to > cln-list-request(a)ginac.de > > You can reach the person managing the list at > cln-list-owner(a)ginac.de > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of CLN-list digest..." > > > Today's Topics: > > 1. exponants of floating points (or do I need it ?) (Franck Sabati?) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 14 Jan 2011 09:40:23 +0100 > From: Franck Sabati? <fsabatie(a)gmail.com> > To: cln-list(a)ginac.de > Subject: [CLN-list] exponants of floating points (or do I need it ?) > Message-ID: > <AANLkTinG7GBxWw0We=8tGv_DNW7CoHK3L-UcFNqaU7iF(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Hello, > > I have a hard time seeing how to best implement the following : I have a > complicated equation of 4 variables that I need evaluated with a certain > number of digits in precision (30 or so) : > > dummy = (15/4)*(1/(pow(xi,5)*(1+i-k)*(2+i-k)*(3+i-k)*(4+i-k)*(5+i-k)))*( > ( (3*pow(3+i-k,2)-(1+i-k)*(5+i-k))*pow(xi*xi-xb,2) > -2*(2+i-k)*(4+i-k)*(1-xi*xi)*(xb*xb-xi*xi) ) * ( pow((xb+xi)/(1+xi),3+i-k)- > pow((xb-xi)/(1-xi),3+i-k) ) + > 6*(3+i-k)*xi*(1-xb)*(xi*xi-xb)*(pow((xb+xi)/(1+xi),3+i-k)+ > pow((xb-xi)/(1-xi),3+i-k)) ); > > > where the variables are i, k, xb and xi. Although most of the equation is > relatively easy to implement in CLN, I have a hard time figuring out how to > deal with the "pow". There appears to be no such exponant function for > floats ? Also, it is still unclear to me how to enter this equation > painlessly into my code. Do I have to put cl_float(1,30) when I use 1 for > instance ? (same for other numbers). Should I evaluate this in floating > precision or with "real numbers" ? And finally how do I evaluate expressions > such as (something)^(something else) with the (something else) NOT being an > integer. > > Sorry about the very basic questions but most of the demo examples deal with > integers and the few that do not, have very little information in them. > > Thank you for your time, > Franck >

1 0

Re: [CLN-list] CLN-list Digest, Vol 55, Issue 2
by Joshua Friedman 10 Jan '11

10 Jan '11

Is there an easy to to tell if I succeeded in compiling cln to 64 bit on OS X? On Sat, Dec 11, 2010 at 6:00 AM, <cln-list-request(a)ginac.de> wrote: > Send CLN-list mailing list submissions to > cln-list(a)ginac.de > > To subscribe or unsubscribe via the World Wide Web, visit > https://www.cebix.net/mailman/listinfo/cln-list > or, via email, send a message with subject or body 'help' to > cln-list-request(a)ginac.de > > You can reach the person managing the list at > cln-list-owner(a)ginac.de > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of CLN-list digest..." > > > Today's Topics: > > 1. Re: Does someone see a way to speed this up (Bruno Haible) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 11 Dec 2010 00:50:47 +0100 > From: Bruno Haible <bruno(a)clisp.org> > To: cln-list(a)ginac.de > Cc: Joshua Friedman <crowneagle(a)gmail.com> > Subject: Re: [CLN-list] Does someone see a way to speed this up > Message-ID: <201012110050.49085.bruno(a)clisp.org> > Content-Type: text/plain; charset="utf-8" > > Hi Joshua, > > > I am using a lot of function calls, would passing constant pointers speed > up > > my program significantly? ... > > I don't expect people to read the code, just look at the way I am using > cln. > > Is this an efficient way? > > Based on past experience, not on your program in particular, here's a list > of suggestions, starting with the most promising ones, and the least > promising ones at the end. > > 1) Use a profiler that displays the bottlenecks and relative consumption > of CPU time. On Linux/x86, the tool of choice is > valgrind <http://valgrind.org/docs/manual/cl-manual.html>, together > with kcachegrind for the visualization of the results. > > 2) Are you using excess precision in large parts of the computations? > Often, in order to get 20 digits of a result, you need 40 digits > throughout the computation if you choose a constant precision. > But sometimes you can actually do part of the computation with > 22 digits and part with 40 digits, and the result will still have > 20 digits of accuracy. > For analyzing this, you need to have a certain understanding of the > numerical stability of your algorithm. For example, when you are > evaluating integrals along a path around a pole, you use 40 > digits in the computation and get 20 digits accuracy in the result; > but when you are evaluating an integrals along a path around a > region with no singularity, then 40 digits in the computation normally > lead to slightly less than 40 digits accuracy in the result. > > 3) When computing integrals, use higher-order integration methods than > simple averaged summation, if the integrand has continuous derivates. > If your integrand is infinitely differentiable, then Runge-Kutta > (the formula with 1/6, 4/6, 1/6) is a big win, and higher-order > Runge-Kutta methods are even better. > > 4) Pull calls to transcendental functions out of loops if possible. > In CLN, elementary operations and sqrt are fast, whereas sin, exp, > log, etc. are significantly slower. > > 5) Pull common computations out of loops. The compiler won't do it > for it. Things like cl_float(2,prec)*PI are quite fast to compute, > but you may save a couple of percent by moving all these "trivial" > common subexpressions outside the loops. > > 6) Then only, start to worry about whether you pass numbers by > "copy" (e.g. cl_N x) or by reference (e.g. const cl_N& x). > CLN numbers are references anyway, so by doing this you can save > the increment and decrement of a reference count, nothing else. > If you had many function calls and many small-sized bignums > and floats, you could save maybe at most 20% by using this trick. > It's more a micro-optimization than the other five points. > > Bruno > > > ------------------------------ > > _______________________________________________ > CLN-list mailing list > CLN-list(a)ginac.de > https://www.cebix.net/mailman/listinfo/cln-list > > > End of CLN-list Digest, Vol 55, Issue 2 > *************************************** > -- Joshua Friedman PhD CrownEagle(a)gmail.com http://www.math.sunysb.edu/~joshua

2 1

Does someone see a way to speed this up
by Joshua Friedman 11 Dec '10

11 Dec '10

I am using a lot of function calls, would passing constant pointers speed up my program significantly? I am doing computations with modular forms (2 variable functions, written as a Fourier series in one variable). I am basically computing integrals of these functions and searching for maximums. The program takes a long time to run, I think much longer than it should. I don't expect people to read the code, just look at the way I am using cln. Is this an efficient way? #include <iostream.h> #include <fstream.h> #include <stdio.h> #include <stdlib.h> #include <string.h> #include <cln/real.h> #include <cln/integer.h> #include <iostream.h> #include <fstream.h> #include <cln/io.h> #include <cln/integer_io.h> #include <cln/real_io.h> #include <cln/float_io.h> #include <cln/float_io.h> #include <cln/integer.h> #include <cln/real.h> #include <cln/complex.h> // We do I/O. #include <cln/io.h> #include <cln/integer_io.h> #include <cln/real_io.h> #include <cln/complex_io.h> // We use the timing functions. #include <cln/timing.h> //using namespace std; using namespace cln; typedef cln::cl_I bigI; typedef cln::cl_F bigF; typedef cln::cl_N bigC; typedef cln::cl_R bigR; bigF h(bigF y, bigF a[], int k, int M, int pr); bigF ff(bigF x, bigF y, bigF a[], int k, int M, int pr); bigF v(bigF y, bigF a[], int k, int M, int pr); bigF integralH(bigF a[], int k, int M, int n, int pr); bigF integralL(bigF a[], int k, int M, int n, int pr); bigF max(bigF a[][100],bigF b[], int k, int D, int M, bigF e, int pr); //for each weight k (2k gives actual weight), make 2d arrary of coefficients //write routine to search FD domain (less than 2k/pi, check the 2), first do //sup for indiv forms, than the average one. int main(int argc, char *argv[]) { //the number of coefficients to use int pbits= 128; int nPartitions = 100; int divisor = 4; float fbinsize = 0.1; int startK = 22; cout << "\n Enter prec, numPartitions, intervalSize, divisor to reduce coef, startK "; cin >> pbits >> nPartitions >> fbinsize >> divisor >> startK; std::string snum; char w; int K,P,D; bigI S, H; std::ifstream myfile ("full_output.txt"); if (myfile.is_open()) { while (! myfile.eof() ) { D = 0; myfile >> w >> K >> P >> D; bigI bigInt; int pr = pbits; cln::float_format_t prec = cln::float_format_t(pr); bigF binSize = cl_float(fbinsize,prec); bigF a[D][100]; bigF b[D]; { for (int n = 0; n < D; n++) { for (int u = 0; u < P; u++) { myfile >> S >> bigInt; a[n][u] = cl_float(bigInt,prec); } if (K >= startK) { b[n] = integralH(a[n],K, P/divisor, nPartitions, pr) + integralL(a[n],K, P/divisor, nPartitions, pr); } //b[n] = cl_float(1,prec); //cout << "\nLoaded Coeff"; } if (K >= startK) { if (D>0) { bigF sNorm = max(a,b,K,D,P, binSize, pr); cout << "\n" << K << "," << sNorm; } } } } myfile.close(); } else cout << "Unable to open file"; return 0; } //Used to do the integral over cylindar y > 1 bigF h(bigF y, bigF a[], int k, int M, int pr) { cln::float_format_t prec = cln::float_format_t(pr); const bigF PI = pi(prec); bigF s = cln::cl_float(0,prec); for (int n = 1; n < M; n++) if (cl_float(4,prec)*PI*cl_float(n,prec)*y < cl_float(300,prec)) s = s + a[n]*a[n]*exp(cl_float(-4,prec)*PI*cl_float(n,prec)*y); return exp(ln(y)*cl_float(2*k-2, prec))*s; } //ff is equiv to g, use it as a test and see which is more effiecient bigF ff(bigF x, bigF y, bigF a[], int k, int M, int pr) { cln::float_format_t prec = cln::float_format_t(pr); const bigF PI = pi(prec); bigF q = cl_float(0.0,prec); int j = 1; bigF w = cl_float(0.0,prec); int r = 0; bigF hh = cl_float(0,prec); for (j = 1; j <= M; j++) { if (cl_float(4,prec)*PI*cl_float(j,prec)*y < cl_float(50,prec)) w = w + a[j]*a[j]*exp(-cl_float(4,prec)*PI*cl_float(j,prec)*y); } for (r = 3; r <= M; r++) { q = cl_float(0.0,prec); for (j = 1; 2*j < r; j++) q = q + a[j]*a[r-j]*cos(cl_float(2,prec)*PI*(cl_float(2,prec)*cl_float(j,prec) - cl_float(r,prec))*x); if (cl_float(2,prec)*PI*cl_float(r,prec)*y < cl_float(50,prec)) q = q*exp(cl_float(-2,prec)*PI*cl_float(r,prec)*y); else q = cl_float(0,prec); w = w + 2*q; } bigF tt = cln::cl_float(2*k,prec); return exp(ln(y)*tt)*w; } bigF max(bigF a[][100], bigF b[], int k, int D, int M, bigF e, int pr) { bigF x,y,half; cln::float_format_t prec = cln::float_format_t(pr); bigF one = cl_float(1,prec); x = cl_float(0,prec); y = cl_float(1,prec); half = cl_float(0.5,prec); bigF K = cl_float(k,prec); bigF mxT = cl_float(0,prec); bigF sum = cl_float(0,prec); bigF w = sqrt(cl_float(3,prec))/cl_float(2,prec); while ( x < half) { y = cl_float(1,prec); while (y < K/cl_float(6,prec)) { sum = cl_float(0,prec); for (int i = 0; i < D; i++) sum = sum + ff(x,y,a[i],k,M,pr)/b[i]; if (sum > mxT) mxT = sum; y = y + e; } x = x + e; } bigF mxB = cl_float(0,prec); x = cl_float(0,prec); while ( x < half) { y = sqrt(cl_float(1) - x*x); while (y < one) { sum = cl_float(0,prec); for (int i = 0; i < D; i++) sum = sum + ff(x,y,a[i],k,M,pr)/b[i]; if (sum > mxB) mxB = sum; y = y + e; } x = x + e; } if (mxB > mxT) return mxB; return mxT; } //twice the Integral f(z) square in x direction from 0 to 1/2; //y is variable from sqrt(3)/2 to 1 bigF v(bigF y, bigF a[], int k, int M, int pr) { cln::float_format_t prec = cln::float_format(pr); const bigF PI = pi(prec); bigF q = cl_float(0.0, prec); bigF w = cl_float(0.0, prec); bigF one = cl_float(1, prec); bigF two = cl_float(2, prec); bigF three = cl_float(3, prec); bigF four = cl_float(4, prec); for (int j = 1; j <= M; j++) { bigF J = cl_float(j,prec); if (four*PI*J*y < cl_float(50,prec)) w = w + a[j]*a[j]*exp(-four*PI*J*y); } w = two*(cl_float(0.5,prec) - sqrt(one-y*y))*w; for (int r = 3; r <= M; r++) { bigF R = cl_float(r,prec); q = cl_float(0,prec); for (int j = 1; 2*j < r; j++) { bigF J = cl_float(j,prec); q = q + a[j]*a[r-j]/(two*PI*(two*J - R))* (sin(PI*(two*J - R)) - sin(two*PI*(two*J - R)*sqrt(one-y*y))); } //q = q + a[j]*a[r-j]*cos(2*PI*(2*j - r)*x); if (two*PI*R*y < cl_float(50,prec)) q = q*exp(-two*PI*R*y); else q = cl_float(0,prec); w = w + two*q; } //for (int n = 1; n < h; n++) s = s + a[n]*exp(2*PI*I*z*n); bigF tt = cln::cl_float((float)(2*k-2),prec); return exp(ln(y)*tt)*w; } bigF integralH(bigF a[], int k, int M, int n, int pr) { const bigF PI = "3.14159265358979323846264338327950288"; cln::float_format_t prec = cln::float_format_t(pr); bigF q = cl_float(0.0, prec); bigF w = cl_float(0.0, prec); bigF aa = cl_float(0.0, prec); bigF b = cl_float(0.0, prec); bigF d = cl_float(0.0, prec); bigF one = cl_float(1, prec); bigF two = cl_float(2, prec); bigF three = cl_float(3, prec); bigF four = cl_float(4, prec); bigF five = cl_float(5, prec); aa = one; //modify this later b = cl_float(k,prec)/two; d = (b-aa)/cl_float(n,prec); //f(x) is the dummy function q = h(aa,a,k,M, pr) + h(b,a,k,M, pr); for (int j = 1; 2*j <= n-2; j++) w = w + h(aa + (two*cl_float(j,prec))*d,a,k,M, pr); q = q + two*w; w = cl_float(0,prec); for (int j = 1; 2*j <= n; j++) w = w + h(aa + (two*cl_float(j,prec)-one)*d,a,k,M, pr); q = q + four*w; q = q*d/three; return q; } bigF integralL(bigF a[],int k, int M, int n, int pr) { cln::float_format_t prec = cln::float_format(pr); bigF q = cl_float(0.0, prec); bigF w = cl_float(0.0, prec); bigF aa = cl_float(0.0, prec); bigF b = cl_float(1.0, prec); bigF d = cl_float(0.0, prec); bigF one = cl_float(1, prec); bigF two = cl_float(2, prec); bigF three = cl_float(3, prec); bigF four = cl_float(4, prec); bigF five = cl_float(5, prec); aa = sqrt(three)/two; b = one; d = (b-aa)/cl_float(n,prec); //f(x) is the dummy function q = v(aa,a,k,M, pr) + v(b,a,k,M, pr); for (int j = 1; 2*j <= n-2; j++) w = w + v(aa + two*cl_float(j,prec)*d,a,k,M, pr); q = q + two*w; w = cl_float(0,prec); for (int j = 1; 2*j <= n; j++) w = w + v(aa + (two*cl_float(j,prec)-one)*d,a,k,M, pr); q = q + four*w; q = q*d/three; return q; } -- Joshua Friedman PhD CrownEagle(a)gmail.com http://www.math.sunysb.edu/~joshua

2 1

Re: [CLN-list] CLN
by Richard B. Kreckel 08 Nov '10

08 Nov '10

Dear Kåre, On 10/29/2010 09:33 PM, Kåre Olaussen wrote: > Dear Bruno and Richard. > > I tried to subscribe to the CLN mailing list, but the subscription page does not work. Indeed. There is a problem in openSUSE 11.3 packaging leading to spurious error messages after subscribing. (The version of MailMan in openSUSE 11.3 is less than the one in previous versions of openSUSE and not up-to-date for Python 2.6 shipped with openSUSE 11.*.) This will be fixed some time. Anyway, you are subscribed now, according to the members list. I'm CC'ing you just to be sure. > We want to use CLN to make world record calculations in various directions, > but this requires effective use of parallel (or multi-thread) computations. > > Our experience with (careful) multithread computations using CLN is spotless. > I have read your reservation in the mailing list archive that CLN is not 'thread-safe'. > > I don't understand what you mean by that! All CLN objects are reference-counted and this code is susceptible to race-conditions where different threads manipulate the reference counts without paying attention to concurrent threads. > Is that like declaring that "no car is collision safe" (in which sense I think no software > package is 'thread-safe'), or are there specific problem areas you know of? > We have already made many-many-million CLN algebraic operations both in > serial and threaded (using openMP), and always obtained the same results. That is entirely possible if you don't share CLN objects between threads or do so only very infrequently. The most likely problem you'll ever hit is probably a segmentation fault while accessing memory another thread as prematurely freed. > We also want to run our (CLN) programs on large computer clusters -- using MPI. > This could be made much more efficient if we could transfer for instance CL_F > numbers between different computers in binary form. > > Could you please help me understand how this can be done? Instead of hacking something low-level, let me make a suggestion: convert your floating-point numbers to a triple (sign, mantissa, exponent), serialize it to a buffer in base 32, send that across your cluster and let the receiver reverse this operation. Check the manual for the type decoded_float. The point is that conversion to/from a base that is a power of two is O(N) with little overhead. It is significantly more efficient than conversion to any other base which is O(N*log(N)*log(log(N))). And 32 is the largest supported base. This way, the overhead compared to the internal representation is only 1-log2(32)/log2(256)=3/8. Unless your program is extremely communication-intense (in which case you should think about redesigning it), this is almost optimal. Best wishes -richy. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>

1 0

Linker error: multiple definition of...
by Lorne McIntosh 06 Oct '10

06 Oct '10

Hi all, I'm attempting to build CLN 1.3.1 on Windows 7 32-bit using MinGW (and NetBeans IDE 6.9.1). The build seemed to be going well until linking, when I got the following error: ../src/.libs/libcln.a(cl_F_decode.o):cl_F_decode.cc:(.text$_ZN3cln12decode_floatERKNS_5cl_LFE[cln::decode_float(cln::cl_LF const&)]+0x0): multiple definition of `cln::decode_float(cln::cl_LF const&)' ../src/.libs/libcln.a(cl_LF_decode.o):d:\School\Grad\Thesis\Prototype09\cln-1.3.1\src/float/lfloat/misc/cl_LF_decode.cc:19: first defined here

2 1

Difficulty building rpm of CLN
by galen 22 Sep '10

22 Sep '10

I am working on an updated rpm of your library for PCLinuxOS. I was hoping for some assistance, as I am running out of ideas. I appreciated any attention you may offer to this problem. Thanks for this library. First of all, I can build the library successfully using ./configure, make, and make install. However, when I try to build it with a spec file, it comes up with all kinds of errors. I started with the spec PCLinuxOS spec file for 1.1.13, then looked for clues from Mandriva, Fedora and the included cln.spec.in. Here seems to be the environment created for building the rpm: + CFLAGS='-O2 -pipe -Wp,-D_FORTIFY_SOURCE=2 -fomit-frame-pointer -march=i586 -mtune=generic -fasynchronous-unwind-tables' + export CFLAGS + CXXFLAGS='-O2 -pipe -Wp,-D_FORTIFY_SOURCE=2 -fomit-frame-pointer -march=i586 -mtune=generic -fasynchronous-unwind-tables' + export CXXFLAGS + FFLAGS='-O2 -pipe -Wp,-D_FORTIFY_SOURCE=2 -fomit-frame-pointer -march=i586 -mtune=generic -fasynchronous-unwind-tables' Here are the errors when I build without any other modifications: checking whether build environment is sane... /usr/share/automake-1.11/am/depend2.am: am__fastdepCCAS does not appear in AM_CONDITIONAL /usr/share/automake-1.11/am/depend2.am: The usual way to define `am__fastdepCCAS' is to add `AM_PROG_AS' /usr/share/automake-1.11/am/depend2.am: to `configure.ac' and run `aclocal' and `autoconf' again. src/Makefile.am: Preprocessed Assembler source seen but `CCAS' is undefined src/Makefile.am: The usual way to define `CCAS' is to add `AM_PROG_AS' src/Makefile.am: to `configure.ac' and run `aclocal' and `autoconf' again. src/Makefile.am: Preprocessed Assembler source seen but `CCASFLAGS' is undefined src/Makefile.am: The usual way to define `CCASFLAGS' is to add `AM_PROG_AS' src/Makefile.am: to `configure.ac' and run `aclocal' and `autoconf' again. yes checking for a thread-safe mkdir -p... /bin/mkdir -p checking for gawk... gawk checking whether make sets $(MAKE)... yes checking for gcc... gcc checking whether the C compiler works... yes checking for C compiler default output file name... a.out checking for suffix of executables... checking whether we are cross compiling... no checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc accepts -g... yes checking for gcc option to accept ISO C89... none needed checking for style of include used by make... GNU checking dependency style of gcc... gcc3 checking how to run the C preprocessor... gcc -E checking for g++... g++ checking whether we are using the GNU C++ compiler... yes checking whether g++ accepts -g... yes checking dependency style of g++... gcc3 checking how to run the C++ preprocessor... g++ -E ./configure: line 4759: CL_AS_UNDERSCORE: command not found ./configure: line 4760: CL_AS_NOEXECSTACK: command not found ./configure: line 4762: syntax error near unexpected token `autoconf' ./configure: line 4762: ` CL_CANONICAL_HOST(autoconf)' make: *** [config.status] Error 2 make: *** Waiting for unfinished jobs.... make: *** [Makefile.in] Error 1 error: Bad exit status from /home/galen/src/tmp/rpm-tmp.68019 (%build) RPM build errors: Bad exit status from /home/galen/src/tmp/rpm-tmp.68019 (%build) From this output, I've tried adding: echo 'AM_PROG_AS' >> configure.ac aclocal autoconf just before %configure. This leads to other errors which I'll gladly outline if necessary. Thank you for taking the time to read this, and any feedback would be appreciated. Galen

2 1

Conversion bug?
by Jens Vollinga 21 Sep '10

21 Sep '10

Hi, is there maybe a bug in the implementation of double_approx()? If I do something like #include <iostream> #include <cln/cln.h> using namespace std; using namespace cln; int main() { cl_F c = cl_F("1.23456789012345678901234567890_200"); cout << c << endl; cout << double_approx(c) << endl; return 0; } the second output gives only single precision (aka float_approx)! Regards, Jens

2 1

CLN patch 2
by Jan private 11 Sep '10

11 Sep '10

Hello Alexei, here is another hairy one. 1. Problem: MSVC does not tag "extern cl_private_thing cl_I_constructor_from_UL (uint32 wert);" in cln/number.h with the cln namespace because it is declared inside the body of an inlined function (symbol created is >>?cl_I_constructor_from_UL@@YAPAUcl_anything@cln@@K@Z<<) But cl_I_constructor_from_UL is defined inside cln namespace in cl_I_from_UL.cc (symbol created is >>?cl_I_constructor_from_UL@cln@@YAPAUcl_anything@1@K@Z<<) So we get unresolved external symbols on every assignment and constructor operation. 2. Solution (rather, work-around): Create an auxiliary function cln::cl_private_thing cl_I_constructor_from_UL (uint32 wert); outside of namespace cln which references cln::cl_private_thing cln::cl_I_constructor_from_UL(wert); Note: The attached patch is exemplary only. As soon as the method to solve this issue is agreed on, similar patches will have to be made for other constructors Best regards, Jan

1 0