[BACK]Return to bfct CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / lib

Diff for /OpenXM_contrib2/asir2000/lib/bfct between version 1.7 and 1.9

version 1.7, 2000/12/14 01:38:37 version 1.9, 2000/12/14 09:36:17
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/lib/bfct,v 1.6 2000/12/13 05:37:31 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/lib/bfct,v 1.8 2000/12/14 09:13:37 noro Exp $
 */  */
 /* requires 'primdec' */  /* requires 'primdec' */
   
Line 54 
Line 54 
 def ann(F)  def ann(F)
 {  {
         V = vars(F);          V = vars(F);
         W = append([y1,y2,t],V);  
         N = length(V);          N = length(V);
         B = [1-y1*y2,t-y1*F];          D = newvect(N);
   
           for ( I = 0; I < N; I++ )
                   D[I] = [deg(F,V[I]),V[I]];
           qsort(D,compare_first);
           for ( V = [], I = N-1; I >= 0; I-- )
                   V = cons(D[I][1],V);
   
         for ( I = N-1, DV = []; I >= 0; I-- )          for ( I = N-1, DV = []; I >= 0; I-- )
                 DV = cons(strtov("d"+rtostr(V[I])),DV);                  DV = cons(strtov("d"+rtostr(V[I])),DV);
   
           W = append([y1,y2,t],V);
         DW = append([dy1,dy2,dt],DV);          DW = append([dy1,dy2,dt],DV);
   
           B = [1-y1*y2,t-y1*F];
         for ( I = 0; I < N; I++ ) {          for ( I = 0; I < N; I++ ) {
                 B = cons(DV[I]+y1*diff(F,V[I])*dt,B);                  B = cons(DV[I]+y1*diff(F,V[I])*dt,B);
         }          }
Line 79  def ann(F)
Line 89  def ann(F)
   
 def indicial1(F,V)  def indicial1(F,V)
 {  {
         V = vars(F);  
         W = append([y1,t],V);          W = append([y1,t],V);
         N = length(V);          N = length(V);
         B = [t-y1*F];          B = [t-y1*F];
Line 241  def weyl_gennf(G,TL,V,O)
Line 250  def weyl_gennf(G,TL,V,O)
                 NF = remove_cont(NF);                  NF = remove_cont(NF);
                 H = cons(NF,H);                  H = cons(NF,H);
         }          }
         print("");          if ( dp_gr_print() ) print("");
         TNF = time()[0]-T0;          TNF = time()[0]-T0;
         if ( dp_gr_print() )          if ( dp_gr_print() )
                 print("gennf(TAB="+rtostr(TTAB)+" NF="+rtostr(TNF)+")");                  print("gennf(TAB="+rtostr(TTAB)+" NF="+rtostr(TNF)+")");
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.9
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>