Changeset bd0aeb in git for Singular/LIB/lejeune.lib


Ignore:
Timestamp:
Jun 23, 2005, 3:52:12 PM (18 years ago)
Author:
Nadine Cremer <cremer@…>
Branches:
(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
Children:
a868907fc6dd9b59a0d27bd87056410c92a84223
Parents:
e7ff6bef3e40732ca87472519b69d1fc7a29d922
Message:
*cremer: nothing new except for comments


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

Legend:

Unmodified
Added
Removed
  • Singular/LIB/lejeune.lib

    re7ff6b rbd0aeb  
    11
    22//-*- mode:C++;-*-
    3 // $Id: lejeune.lib,v 1.4 2005-06-23 13:40:16 cremer Exp $
     3// $Id: lejeune.lib,v 1.5 2005-06-23 13:52:12 cremer Exp $
     4
    45
    56info="
     
    1415    ringchange(i);       changes the ring to the one needed in ith step     
    1516    plugin_coeffs(i,f)   plugs tpolys into f, up to power i
    16     maxidealstep(i,N);   returns ideal needed for contraction in ith step;
     17    maxidealstep(i,N);   returns ideal needed for contraction in ith step
    1718                         N is number of variables of input f
    1819    formaldiff(f,k);     computes the formal derivatives D_I with |I|<k
     
    5556     resultdiff=resultdiff,J;
    5657   }
    57   resultdiff;~
    58   export(resultdiff);
    59   return(R);
     58  export(resultdiff);                // exports the result
     59  return(R);                         // return the ring
    6060 
    6161}
     
    6666{
    6767  def r=basering;
    68   def R=ringchange(i);
    69   setring R;
     68  def R=ringchange(i);              // changes the ring
     69  setring R;                        // makes it new basering;
    7070  ideal I=tpolys(i,startvar);
    71   poly g=imap(r,f);
    72   export(g);
    73   map h=r,I;
    74   ideal J=h(f);
    75   export(h);
    76   matrix result=coeffs(J[1],t);
    77   export result;
    78   return(R);
     71  poly g=imap(r,f);                 // maps f to new basering
     72  export(g);                        // export it
     73  map h=r,I;                        // define map according to our purpose
     74  ideal J=h(f);                     // gives f with power series plugged in
     75  export(h);                   
     76  matrix result=coeffs(J[1],t);     // gives the t-coefficients
     77  export result;                    // export it i.o. to use it later on
     78  return(R);                        // return ring (ring change!)
    7979}
    8080
     
    8484{
    8585  int startvar=nvars(basering);
    86   export(startvar);
     86  export(startvar);         
    8787  string str=variables(startvar,i);
    88   def R=changevar(""+varstr(r)+",t,"+variables(startvar,i)+"");
    89   return(R);
     88  def R=changevar(""+varstr(r)+",t,"+variables(startvar,i)+"");  // change
     89  return(R);                 // return the ring, needed in future proc
    9090}
    9191
     
    101101     for (s=1;s<=i;s++)
    102102     {
    103        str=""+a_z(u)+"("+string(s)+")";    
    104        l[(u-1)*i+s]=str;
     103       str=""+a_z(u)+"("+string(s)+")"; // creates new variables   
     104       l[(u-1)*i+s]=str;                // saves them in a list
    105105     }
    106106   }
    107107  //l=insert(l,"t");
    108   string str1=string(l);
    109   return(str1);
     108  string str1=string(l);                // makes the list into a string,
     109  return(str1);                         // (needed for ring change)
    110110}
    111111
    112112
    113 proc a_z (int n)                      // returns nth letter of the alphabet
     113proc a_z (int n)                        // returns nth letter of the alphabet
    114114{
    115115  if((n<1)||(n>26))                     // input admissible?
     
    127127proc tpolys (int i,int k)             // constructs polynomials a(1)*t+...
    128128{                                     // has to be called from tpolys
    129   int s,t;                             // loop variables
     129  int s,t;                            // loop variables
    130130  int v;           
    131131  poly sum;
     
    147147
    148148
    149 proc maxidealstep (int i,int N)   
    150 {
     149proc maxidealstep (int i,int N)       // returns ideal needed for
     150{                                     // differentiation in ith step
    151151  ideal I=var(N+1+i);
    152152  int j;
Note: See TracChangeset for help on using the changeset viewer.