Changeset 565d9b in git for Singular/LIB/lejeune.lib
 Timestamp:
 Jun 24, 2005, 2:01:09 PM (18 years ago)
 Branches:
 (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
 Children:
 8761f2832a06650bd36ab1c7aeff785ae639c7bd
 Parents:
 ee17bf95ed1666f2faa03fb99a49b4518f2e3473
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

Singular/LIB/lejeune.lib
ree17bf r565d9b 1 1 2 2 //* mode:C++;* 3 // $Id: lejeune.lib,v 1. 9 20050624 09:25:08cremer Exp $3 // $Id: lejeune.lib,v 1.10 20050624 12:01:09 cremer Exp $ 4 4 5 5 … … 18 18 N is number of variables of input f 19 19 formaldiff(f,k); computes the formal derivatives D_I with I<k 20 f_set (f,H); returns the set F corresponding to H as described by20 f_setstep(f,H); iterates the steps given by H, saved in f_set f_set(f,H); returns the set F corresponding to H as described by 21 21 M. Lejeune 22 22 "; … … 28 28 29 29 30 proc f_set (poly f,intvec H) 30 proc f_set (poly f,intvec H) // puts together the single steps from 31 { // f_setstep 32 def r=basering; 33 int b=size(H); 34 int i; 35 def R=ringchange(b1); 36 setring R; 37 list l; 38 for(i=2;i<=b;i++) 39 { 40 setring r; 41 def tmp=f_setstep(f,intvec(H[1..i])); 42 setring R; 43 ideal I=imap(tmp,resultf_set); 44 l[i1]=I; 45 } 46 l; 47 return(R); 48 } 49 50 51 proc f_setstep (poly f,intvec H) 31 52 { 32 53 int p; // loop variable … … 47 68 } 48 69 } 49 def r=basering;50 for(p=b1;p>1;p) // iterating steps51 {52 if(p==2)53 {54 def tmp=f_set(f,intvec(H[1..p]));55 setring tmp;56 def tmp1=resultf_set;57 export(tmp1);58 //def tmp1=resultf_set;59 //setring r;60 }61 70 def R=formaldiff(f,b1,c,H[b]); // actual step 62 71 setring R; 63 def T=resultdiff; 64 ideal resultf_set; 65 //if(b==2) 66 // { 67 // ideal resultf_set; 68 // resultf_set=T; 69 // export(resultf_set); 70 // } 71 //else 72 // { ideal resultf_set=imap(tmp,tmp1); 73 // resultf_set=resultf_set,T; 74 // export(resultf_set); 75 // } 76 77 resultf_set; 72 def resultf_set=resultdiff; 73 export(resultf_set); 78 74 return(R); 79 75 }
Note: See TracChangeset
for help on using the changeset viewer.