Changeset 9173792 in git for Singular/LIB/ntsolve.lib
- Timestamp:
- Feb 19, 2001, 3:17:50 PM (23 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38dfc5131670d387a89455159ed1e071997eec94')
- Children:
- 25c431d26faee1a17b854e49114f755c02aaf1ef
- Parents:
- 584d1f1ad4878b9855984324f19d6383aa60632c
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/ntsolve.lib
r584d1f1 r9173792 1 1 //(GMG, last modified 16.12.00) 2 2 /////////////////////////////////////////////////////////////////////////////// 3 version="$Id: ntsolve.lib,v 1.1 1 2001-01-16 13:48:35 SingularExp $";3 version="$Id: ntsolve.lib,v 1.12 2001-02-19 14:17:47 lossen Exp $"; 4 4 category="Symbolic-numerical solving"; 5 5 info=" … … 9 9 10 10 PROCEDURES: 11 nt_solve(G, ..); find one real root of 0-dimensional ideal G12 triMNewton(G, ..); find one real root for 0-dim triangular system G11 nt_solve(G,ini,[..]); find one real root of 0-dimensional ideal G 12 triMNewton(G,a,[..]); find one real root for 0-dim triangular system G 13 13 "; 14 14 … … 16 16 17 17 proc nt_solve (ideal gls, ideal ini, list #) 18 "USAGE: nt_solve(gls,ini[,ipar]); gls ideal, ini ideal, ipar list or intvec 19 gls: contains the equations, for which a solution will be computed 20 ini: ideal of initial values (approximate solutions to start with) 21 ipar: control integers (default: ipar = 100,10) 22 ipar[1] - max. number of iterations 23 ipar[2] - accuracy (we have the l_2-norm ||.||): accept solution 24 sol if ||gls(sol)|| < eps0*(0.1^ipar[2]) where eps0 = ||gls(ini)|| 25 is the initial error 18 "USAGE: nt_solve(gls,ini[,ipar]); gls,ini= ideals, ipar=list/intvec,@* 19 gls: contains the equations, for which a solution will be computed 20 ini: ideal of initial values (approximate solutions to start with),@* 21 ipar: control integers (default: ipar = 100,10) 22 @format 23 ipar[1]: max. number of iterations 24 ipar[2]: accuracy (we have the l_2-norm ||.||): accept solution @code{sol} 25 if ||gls(sol)|| < eps0*(0.1^ipar[2]) 26 where eps0 = ||gls(ini)|| is the initial error 27 @end format 26 28 ASSUME: gls is a zerodimensional ideal with nvars(basering) = size(gls) (>1) 27 29 RETURN: ideal, coordinates of one solution (if found), 0 else 28 NOTE: printlevel >0: displays comments (default =0)30 NOTE: if printlevel >0: displays comments (default =0) 29 31 EXAMPLE: example nt_solve; shows an example 30 32 " … … 254 256 255 257 proc triMNewton (ideal G, ideal a, list #) 256 "USAGE: triMNewton(G,a[,ipar]); G ideal, a ideal, ipar list orintvec257 ASSUME: G: g1,..,gn, a triangular system of n equations in n vars, 258 i.e gi=gi(var(n-i+1),..,var(n))259 260 zero of G to start with (with a[i] to be substituted in var(i))261 262 ipar[1] - max. number of iterations263 ipar[2] - accuracy (we have as norm |.| absolute value ):264 accept solution sol if |G(sol)| < |G(a)|*(0.1^ipar[2])265 itb: int, the maximal number of iterations performed266 err: number, an error bound258 "USAGE: triMNewton(G,a[,ipar]); G,a= ideals, ipar=list/intvec 259 ASSUME: G: g1,..,gn, a triangular system of n equations in n vars, i.e. 260 gi=gi(var(n-i+1),..,var(n)),@* 261 a: ideal of numbers, coordinates of an approximation of a common 262 zero of G to start with (with a[i] to be substituted in var(i)),@* 263 ipar: control integer vector (default: ipar = 100,10) 264 @format 265 ipar[1]: max. number of iterations 266 ipar[2]: accuracy (we have as norm |.| absolute value ): 267 accept solution @code{sol} if |G(sol)| < |G(a)|*(0.1^ipar[2]). 268 @end format 267 269 RETURN: an ideal, coordinates of a better approximation of a zero of G 268 270 EXAMPLE: example triMNewton; shows an example
Note: See TracChangeset
for help on using the changeset viewer.