Re: [sage-devel] Fwd: [mpir-devel] New MPIR-related project

2010-05-05 Thread Sergey Bochkanov
Hello, Gonzalo. You wrote 3 мая 2010 г., 20:46:41: Although I'm more of a GPL (v2) type, and I wouldn't write BSD code unless I had a very good reason, I could consider relicensing some interface code like this. I've actually relicensed this code to LGPL v2+ (for gmpy), but I'm not sure I

Re: [sage-devel] Re: Fwd: [mpir-devel] New MPIR-related project

2010-05-05 Thread Sergey Bochkanov
Hello, rjf. You wrote 4 мая 2010 г., 19:21:42: Note that converting (in applications of these Lisp programs) to gmp integers and back is not something you do very often. It depends on problems you solve. It may be one of the short operands and long accumulator kind. Or it just

Re: [sage-devel] Support Vector Machine

2010-05-05 Thread Sergey Bochkanov
Hello, William. You wrote 4 мая 2010 г., 22:34:58: Could you write a little more to sage-devel about why mathematicians might care about support vector machines -- it's possible that most people reading this have never heard of them. Being sage-devel reader for a while I may say that

Re: [sage-devel] Support Vector Machine

2010-05-06 Thread Sergey Bochkanov
Hello, William. You wrote 5 мая 2010 г., 16:43:39: Yes, numerical analysis is definitely really, really needed in Sage. In fact, it is likely to be by far the most important part of Sage in the long run. Ithinkthat I should look in this direction... I hope to implement Python

Re: [sage-devel] Support Vector Machine

2010-05-08 Thread Sergey Bochkanov
Hello, William by the way, what do you think about using in SAGE both double and multiple precision versions of ALGLIB? I have fully functional multiple precision version of ALGLIB which relies on MPFR for floating point computations. It includes multiple precisionlinear

Re: [sage-devel] Support Vector Machine

2010-05-08 Thread Sergey Bochkanov
Hello, William. My first thought is that it would be *awesome* to have _efficient_ multiprecision linear algebra, etc., in Sage. This is the sort of thing people expect from Sage, but are often disappointed to find that we don't have. As for the efficiency, I may say that

Re: [sage-devel] Support Vector Machine

2010-05-09 Thread Sergey Bochkanov
Hello, Francesco. Just FYI, you might want to look into C++ expression templates - they might allow you to retain the convenient notation with no overhead. I know about expression templates, but even with templates sometimes creation of temporary is the only way to calculate value of a

Re: [sage-devel] Support Vector Machine

2010-05-09 Thread Sergey Bochkanov
Hello, Francesco. You wrote 9 мая 2010 г., 20:06:24: If the expression is to be intended like this: for i in (0,size): x[i] =((a[i] + b[i]) * (c[i] + d[i]) - e[i]) / (f[i] + 3) (i.e., element-wise), then I'm pretty confident expression templates could eliminate all temporaries. I am

Re: [sage-devel] Re: GMP web sites says Sage relies on GMP

2010-05-17 Thread Sergey Bochkanov
Hello, William. You wrote 17 мая 2010 г., 22:13:46: I have written to them about 20 times to get this changed. Did they ever notice that license changed to LGPL 3, so there is nothing to complain now? :( BTW, MPIR team is thinking on developing whole new BSD-licensed library.

Re: [sage-devel] Re: GMP web sites says Sage relies on GMP

2010-05-18 Thread Sergey Bochkanov
Hello, Bill. You wrote 18 мая 2010 г., 15:37:56: The FSF notified us that it was their intention to try and get Torbjorn to remove the outdated statements from his website about MPIR. He hasn't done so. And he never mentioned that license was changed to LGPL 3 mostly due to his objections

Re: [sage-devel] Re: GMP web sites says Sage relies on GMP

2010-05-18 Thread Sergey Bochkanov
Hello, just want to say that at least I can sign such petition (if one day you decide it is needed). Or to write this letter. Good luck! I hope this issue will be resolved once and for all. -- With best regards, Sergey mailto:sergey.bochka...@alglib.net -- To

Re: [sage-devel] Re: GMP web sites says Sage relies on GMP

2010-05-20 Thread Sergey Bochkanov
Hello, Bill. One eventually gets fed up, and I'm just announcing what I already planned to do months ago. I think you took the wrong hint. You did impressive and important job which is appreciated by both SAGE and bignum communities. I know that people sometimes get tired, but I hope that

Re: [sage-devel] multiprecision linear algebra

2010-05-21 Thread Sergey Bochkanov
Hello, Jason. You wrote 21 мая 2010 г., 3:48:03: A while ago someone talked about incorporating ALGLIB into Sage to provide, among other things, multiprecision linear algebra. It was me :) I am working on this issue. I expect that interface between double precision ALGLIB and Python

Re: [sage-devel] multiprecision linear algebra

2010-05-21 Thread Sergey Bochkanov
correction to my previous message: It leads to unnecessary allocation of temporaries and some performance penalty (about 25-30% for 128-bit precision). Sorry, I made a mistake when estimating MPLAKACK performance. 25-30% are ALGLIB's, not MPLAPACK's. I've never measured MPLAPACK

[sage-devel] Re: [mpir-devel] Re: MPIR documentation and the GMP website

2010-05-25 Thread Sergey Bochkanov
Hello, Bill. Now, we are looking into why this isn't appearing on the title page of the pdf documentation. That appears to be a technical issue! We will also add The GMP Development Team as per the latest GMP documentation. I think that when you fix _everything_ someone still will be able

[sage-devel] Re: [mpir-devel] Re: MPIR documentation and the GMP website

2010-05-25 Thread Sergey Bochkanov
Hello, Bill. So it is established that we are being accused of academic fraud for removing the author information WHEN IT WAS DONE BY SOMEONE ON THE GMP TEAM in the actual version of GMP we forked from!! Confirmed. I think what is currently on the GMP website is FAR more damaging to the

Re: [sage-devel] Re: [mpir-devel] Re: MPIR documentation and the GMP website

2010-05-25 Thread Sergey Bochkanov
Hello, Bill. In fact, unless I am mistaken refresh your browser! Now let's all forget this ever happened and let's get on with our jobs. +1! -- With best regards, Sergey mailto:sergey.bochka...@alglib.net -- To post to this group, send an email to

Re: [sage-devel] Make Cephes a standard spkg

2010-05-26 Thread Sergey Bochkanov
Hello, David. Is license info available? I've contacted Stephen Moshier several years ago and he made it clear that anyone can use Cephes under any license (as long as copyright is preserved). I think that someone can contact him and ask for the explicit permission to license Cephes

[sage-devel] ANSI C vs. C99

2010-06-08 Thread Sergey Bochkanov
Hello, SAGE developers! I have a question about C standard which SAGE packages should conform to. Is it C99? Plain ANSI C? Modern ANSI C, i.e. ANSI C which borrowed several features from C99? I am working on Python=ALGLIB interface with intention to make it usable from SAGE. It

[sage-devel] Re: ANSI C vs. C99

2010-06-09 Thread Sergey Bochkanov
Hello, thanks to all who've replied to this discussion! Taking into account what was said above, I've decided to target pure ANSI C, i.e. C without newer constructs like // comments and other stuff from newer standards. I don't want to use C99 because it isn't supported by MSVC (de facto

Re: [sage-devel] Re: ANSI C vs. C99

2010-06-10 Thread Sergey Bochkanov
Hello, Robert. You wrote 10 июня 2010 г., 7:58:49: I'd be glad to know exactly what is special about writing ALGLIB that's more difficult to do with using Cython. I can't say that it is difficult to do with Cython. I just want to say that for me - and just for me - Cython has no benefits over

Re: [sage-devel] Re: ANSI C vs. C99

2010-06-10 Thread Sergey Bochkanov
Hello, David You wrote 10 июня 2010 г., 15:25:03: It would be great if you could check your code with the Sun Studio compiler, which you can get for free for Linux or Solaris. It is not great, it is must-have. Building ALGLIB on every platform where SAGE builds and with every popular

Re: [sage-devel] Re: ANSI C vs. C99

2010-06-10 Thread Sergey Bochkanov
Hello, Having checked again, __sparc is preferable. That is supported by gcc, g++, cc and CC (the Sun compilers) Thanks a lot! I've also got machines running AIX, tru64, HP-UX and IRIX. None of the CPUs in those sorts of machines would support SSE. Clearly there are many system which

Re: [sage-devel] MPIR secrets and bsdnt

2010-06-24 Thread Sergey Bochkanov
Hello, Bill. You wrote 22 июня 2010 г., 3:30:30: BSDNT .. scratch. One idea I had was to write all the low level routines in LLVM instead of assembly, then just build suitable optimisations into LLVM so that nice assembly is automatically produced. This is entirely feasible and would

[sage-devel] ALGLIB spkg is released (first beta)

2010-07-08 Thread Sergey Bochkanov
Hello! I've talked about making Sageable interface to ALGLIB two months ago. It took twice as much time as I previously thought, but now it is ready to go public. You can download first beta of ALGLIB double precision interface (as .spkg or .zip file) from

Re: [sage-devel] Re: ALGLIB spkg is released (first beta)

2010-07-09 Thread Sergey Bochkanov
Hello, Jason. You wrote 9 июля 2010 г., 15:34:07: On OSX 10.6.4, I get an error: ... ... This looks like it might be useful, but I don't know: http://lists.apple.com/archives/unix-porting/2003/Oct/msg00032.html I think that it is easy to fix, thanks for the link you gave me. The only

Re: [sage-devel] ALGLIB spkg is released (first beta)

2010-07-10 Thread Sergey Bochkanov
Hello, David! Thanks for your reply, it was very useful. You wrote 10 июля 2010 г., 1:47:31: 1) It does not work on Solaris 10 on SPARC - see below. I suspect there are some GNU-specific linker options, which will fail if the Sun linker is used on Solaris 10. Looks like it is actually

[sage-devel] spkg-check

2010-07-16 Thread Sergey Bochkanov
Hello, I have a question about spkg-check file. When does it executed - after spkg-install call or independently of it? -- With best regards, Sergey mailto:sergey.bochka...@alglib.net -- To post to this group, send an email to sage-devel@googlegroups.com To

Re: [sage-devel] spkg-check

2010-07-16 Thread Sergey Bochkanov
Hello, William. You wrote 16 июля 2010 г., 17:12:53: I have a question about spkg-check file. When does it executed - after spkg-install call or independently of it? After. So if spkg-check is executed, I can be sure that spkg was successfully installed into system? -- With best regards,

Re: [sage-devel] spkg-check

2010-07-16 Thread Sergey Bochkanov
Hello, Dr.. You wrote 16 июля 2010 г., 17:37:25: So if spkg-check is executed, I can be sure that spkg was successfully installed into system? Not exactly, as some packages unfortunately ignore errors, so if the package fails to install, you will not know it. I am talking about my own

[sage-devel] ALGLIB spkg is released (second version)

2010-07-26 Thread Sergey Bochkanov
Hello! Second beta of ALGLIB.spkg is ready. It can be downloaded from http://www.alglib.net/share/2010-07-26-alglib-for-sage/ Following issues were solved in this release: * compilation on SPARC/Solaris * missing support for SAGE64 * trashing of SAGE_LOCAL with additional dirs Issues

Re: [sage-devel] ALGLIB spkg is released (second version)

2010-07-29 Thread Sergey Bochkanov
Hello, Carl. Here is some documentation about Python-ALGLIB interface (called X-interface in this document): http://www.alglib.net/share/2010-07-26-alglib-for-sage/x-interface.pdf It is very short, but should give a hint about meaning of all these functions, expectations,

Re: [sage-devel] ALGLIB spkg is released (second version)

2010-08-03 Thread Sergey Bochkanov
Hello, Carl. You wrote 3 августа 2010 г., 6:39:49: I've attached a patch to make alglib allow input in the following formats: * boolean vector/matrix: GF(2), RDF, numpy.bool, numpy.float64 * integer vector/matrix: RDF, numpy.int32, numpy.float64 * real: RDF, numpy.float64 * complex: CDF,

Re: [sage-devel] doctest failure with leading zero not printing

2010-08-05 Thread Sergey Bochkanov
Hello, David You wrote 4 августа 2010 г., 14:49:47: 8: sage: maxima('asinh(1.0)') Expected: 0.881373587019543 Got: .8813735870195429 Clearly there is some numerical noise issues, which need to be investigated, but are probably trivial. I think that it is more than

Re: [sage-devel] ALGLIB spkg is released (second version)

2010-08-05 Thread Sergey Bochkanov
Hello, Carl. You wrote 3 августа 2010 г., 6:39:49: On Wed, Jul 28, 2010 at 11:38 PM, Sergey Bochkanov sergey.bochka...@alglib.net wrote: My proposal is to make * boolean vector/matrix = GF(2), RDF (non-zero = True) * integer vector/matrix = RDF * real = RDF * complex = CDF One more

Re: [sage-devel] ALGLIB spkg is released (second version)

2010-08-05 Thread Sergey Bochkanov
Hello, Carl. You wrote 3 августа 2010 г., 6:39:49: On Wed, Jul 28, 2010 at 11:38 PM, Sergey Bochkanov sergey.bochka...@alglib.net wrote: My proposal is to make * boolean vector/matrix = GF(2), RDF (non-zero = True) * integer vector/matrix = RDF * real = RDF * complex = CDF Forgot

Re: [sage-devel] Re: doctest failure with leading zero not printing

2010-08-05 Thread Sergey Bochkanov
Hello, Jason. You wrote 5 августа 2010 г., 20:17:11: (depending on problem and its stability properties). However I don't know how to do it using doctest framework. Python tries to output numbers with full precision,and there is no way to tell doctest framework to compare decimal

Re: [sage-devel] Re: doctest failure with leading zero not printing

2010-08-05 Thread Sergey Bochkanov
Hello, David. You wrote 5 августа 2010 г., 20:44:53: With all the numerical noise issues I've seen in Sage, the three dots solves its. So if we expect 1.00 but get 1.01 we can change that to 1.0... and the test will pass. However, what if we get

[sage-devel] Next version of alglib.spkg (and some thoughts about debugging)

2010-08-09 Thread Sergey Bochkanov
Hello! New version of ALGLIB spkg is released. It can be downloaded from http://www.alglib.net/share/2010-08-09-alglib-for-sage/ New version of spkg: * compiles under OS X, SPARC/Solaris, Linux * was tested wuth 2.7 and 3.1.2 branches of Python * includes patch by Carl Witty

Re: [sage-devel] Next version of alglib.spkg (and some thoughts about debugging)

2010-08-09 Thread Sergey Bochkanov
Hello, Carl. You wrote 9 августа 2010 г., 19:58:12: Add this line at the top of check.py: import sage.all Thanks! I've added it on top of _alglib.py because a) check.py, as every other file, explicitly imports _alglib.py b) I prefer to modify one file instead of two Apparently, it works

Re: [sage-devel] Linear algebra reviewer(s) needed

2010-08-11 Thread Sergey Bochkanov
Hello, Rob. You wrote 11 августа 2010 г., 1:10:25: I have an undergraduate student who is funded this summer to develop some routines for Sage which create random matrices with nice properties for beginning students of linear algebra to practice with (or for faculty to put on exams, or to put

Re: [sage-devel] alglib

2010-10-12 Thread Sergey Bochkanov
Hello, Jason. You wrote 13 октября 2010 г., 8:16:15: I'm curious what the current status of the alglib spkg is---anyone know? Works were paused for a while (I've spent all August and September working on ALGLIB 3.0 for C++ and C#), but I don't want to abandon this project. I plan to