Changeset da408f in git for Singular/LIB/solve.lib

Timestamp:

Jul 8, 1999, 12:18:13 PM (24 years ago)

Author:

Moritz Wenk <wenk@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

b719a30005f8bcbcca9d638b7908dda038ffee56

Parents:

4deddb979be065737d63e926fc05f264a6b1078a

Message:

*wenk: changed uressolve CMD_3 -> CMD_M (4)
	       laguerre  CMD2_ -> CMD_3
       removed "setFloatDigits" in extra.cc
       fixed output (2.2e33 -> 2.2e+33)
       adapted solve.lib to uressolve, laguerre, extended examples


git-svn-id: file:///usr/local/Singular/svn/trunk@3247 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• Singular/LIB/solve.lib (modified) (6 diffs)

Singular/LIB/solve.lib

@@ -1 +1 @@
 ///////////////////////////////////////////////////////////////////////////////
 
-version="$Id: solve.lib,v 1.11 1999-07-07 16:38:28 obachman Exp $";
+version="$Id: solve.lib,v 1.12 1999-07-08 10:18:13 wenk Exp $";
 info="
 LIBRARY: solve.lib     PROCEDURES TO SOLVE POLYNOMIAL SYSTEMS
@@ -67 +67 @@
   }
 
-  int digits= system("setFloatDigits",prec);
-
-  return(uressolve(gls,typ,polish));
+  return(uressolve(gls,typ,prec,polish));
 
 }
@@ -81 +79 @@
   // result is a list (x,y)-coordinates as strings
 
-  // now with complex coefficient field, precision is 10 digits
-  ring rsc= (real,10,I),(x,y),lp;
+  // now with complex coefficient field, precision is 20 digits
+  ring rsc= (real,20,I),(x,y),lp;
   ideal i = (2+3*I)*x2 + (0.35+I*45.0e-2)*y2 - 8, x2 + xy + (42.7)*y2;
-  ures_solve(i);
+  list l= ures_solve(i);
   // result is a list of (x,y)-coordinates of complex numbers
+  l;
+  // check the result
+  subst(subst(i[1],x,l[1][1]),y,l[1][2]);
 }
 ///////////////////////////////////////////////////////////////////////////////
@@ -126 +127 @@
   }
 
-  int digits= system("setFloatDigits",prec);
-
-  return(laguerre(f,polish));
+  return(laguerre(f,prec,polish));
 
 }
@@ -147 +146 @@
   list l = laguerre_solve(f);
   l;
+  // check result, value of substituted poly should be near to zero
+  subst(f,x,l[1]);
+  subst(f,x,l[2]);
 }
 ///////////////////////////////////////////////////////////////////////////////
@@ -244 +246 @@
   ideal p = 2,3;
   ideal v= 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16;
-  interpolate( p,v,3 );
-}
-///////////////////////////////////////////////////////////////////////////////
+  poly ip= interpolate( p,v,3 );
+  ip;
+  // compute poly at point 2,3, result must be 2
+  subst(subst(ip,x,2),y,3);
+  // compute poly at point 2^15,3^15, result must be 16
+  subst(subst(ip,x,2^15),y,3^15);
+}
+///////////////////////////////////////////////////////////////////////////////