Ask Your Question
2

newton's method for multiple variables / arbitrary precision

asked 2014-05-20 18:25:59 +0200

mf gravatar image

I am trying to find a numerical approximation with arbitrary precision to a real solution to a system of multivariate polynomial equations. I start out with an approximation which is somewhat close to solving the system, up to an precision of about 1e-05. (Meaning that the equations that I try to evaluate are not zero, but smaller than 1e-05 for my starting value)

In this question is it recommended to use scipy's fmin_tnc method, which is what I did. This works out very nicely and it quickly gave a new solution which now solves my system with precision 1e-07. In the Scipy doc it is stated that one can set the "epsilon" parameter, but not smaller than machine precision. So it seems like I can't get much more precision with this method?!

Let's say I want to solve my system with precision 1e-250. My questions are:

  1. Can I use the fmin_tnc function to find solutions with higher precision?
  2. I there another way in sage to find real solutions to polynomial systems locally (e.g. with the pari/gp)?
edit retag flag offensive close merge delete

Comments

Could you provide the system? I do not know whether pari/gp can help, but it is not hard to setup a Newton method. From what I saw the method fmin_trunc looks much more powerful than what you need.

vdelecroix gravatar imagevdelecroix ( 2014-06-29 12:26:32 +0200 )edit

@vdelecroix: I figured out how to do this, see my answer below

mf gravatar imagemf ( 2014-07-03 09:07:10 +0200 )edit

@mf: cool! Could you select your answer (that way your question will be seen as answered in the list).

vdelecroix gravatar imagevdelecroix ( 2014-07-03 10:00:36 +0200 )edit

@vdelecroix: will do as soon as I have enough karma to do so..

mf gravatar imagemf ( 2014-07-03 10:14:08 +0200 )edit

1 Answer

Sort by ยป oldest newest most voted
2

answered 2014-07-03 09:06:39 +0200

mf gravatar image

updated 2014-07-03 14:21:08 +0200

calc314 gravatar image

One easy way to this is to use mpmath.findroot!

edit flag offensive delete link more

Comments

Nice! I've added a link to the documentation.

calc314 gravatar imagecalc314 ( 2014-07-03 14:22:08 +0200 )edit
slelievre gravatar imageslelievre ( 2018-11-23 18:14:35 +0200 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2014-05-20 18:25:59 +0200

Seen: 698 times

Last updated: Jul 03 '14