Hi there, As far as I know there is no arbitrary precision within scipy/numpy. On the other hand scipy/numpy is shipped with Sage and if you do for example
sage: m = matrix(RDF, [[2,0],[1,0]]) the matrix m is (in the backend) a numpy matrix stored under m._numpy_matrix. Sadly you can not access this attribute because it is at C level (Cython). You may have a look at sage/matrix/matrix_double_dense.pyx The first thing to do would be to do the list of existing softwares and libraries for numerical things (knowing that ATLAS, numpy/scipy are already there) and start implementing interfaces. As far as I know, there is one nice numerical thing in Sage : root finding of polynomials. The code is in sage/rings/polynomials/real_roots.pyx and sage/rings/polynomial/complex_roots.pyx Best, Vincent 2014-02-13 2:31 UTC+01:00, ref...@uncg.edu <ref...@uncg.edu>: > Oh, FEM brings me back to my master's days. Programming the FEM can > sometimes be painful. Perhaps I'm not the best one to ask about this, but > you can use implement in python using Numpy and Scipy...and then run it in > Sage. > > On Wednesday, February 12, 2014 6:52:37 PM UTC-5, Chris Gorman wrote: >> >> kcrisman and Rick, >> >> I guess that I am wondering if there is a group that is devoted to >> contributing to the numerical aspects of Sage and, if so, what they are >> focusing on. I know that I would be interested in helping implement is FEM >> >> and arbitrary-precision numerical integration. I would also like to work >> on >> the linear algebra, but I do not know whether it would be better to focus >> >> on numpy or scipy at that point. I guess I just want to help improve >> Sage's >> numerical capabilities. >> >> Chris >> >> >> >> On Wednesday, February 12, 2014 11:20:33 AM UTC-8, ref...@uncg.edu wrote: >>> >>> Chris, >>> >>> I'm not sure what you are looking for. But, I'm planning on submitting >>> some personal algorithms to sage. Sometimes I had to make some >>> improvements to built-in functions to suit my purpose. For example, >>> numerical integration built into sage is not arbitrary precision...which >>> I >>> had to make arbitrary for my algorithm because of the highly oscillatory >>> >>> nature of the function I was integrating. >>> >>> Rick >>> >>> On Wednesday, February 12, 2014 11:26:50 AM UTC-5, Chris Gorman wrote: >>>> >>>> Does anyone have know who is working on improving the numerical methods >>>> >>>> in Sage? I am beginning my graduate program in numerical analysis and >>>> would >>>> like to use Sage for my work and research. >>>> >>> > > -- > You received this message because you are subscribed to the Google Groups > "sage-devel" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to sage-devel+unsubscr...@googlegroups.com. > To post to this group, send email to sage-devel@googlegroups.com. > Visit this group at http://groups.google.com/group/sage-devel. > For more options, visit https://groups.google.com/groups/opt_out. > -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To post to this group, send email to sage-devel@googlegroups.com. Visit this group at http://groups.google.com/group/sage-devel. For more options, visit https://groups.google.com/groups/opt_out.
