Changeset 8de0ce in git for Singular/LIB/lejeune.lib


Ignore:
Timestamp:
Jul 4, 2005, 9:58:23 AM (19 years ago)
Author:
Nadine Cremer <cremer@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
2f6f3ee6e6cc26011f65dd0589eebe383ae68f81
Parents:
0540757afbdc88eb339dd7280833cf529b3d7010
Message:
*cremer:removed another bug in g_set


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

Legend:

Unmodified
Added
Removed
  • Singular/LIB/lejeune.lib

    r054075 r8de0ce  
    11
    22//-*- mode:C++;-*-
    3 // $Id: lejeune.lib,v 1.15 2005-07-01 10:32:24 cremer Exp $
     3// $Id: lejeune.lib,v 1.16 2005-07-04 07:58:23 cremer Exp $
    44
    55
     
    109109  matrix coe=Coe;                     // gives the t-coeff. after plugging in
    110110  poly fkv;                           // need this stuff for the following
    111   ideal m;                            // loops...
    112   ideal m1,m2,J,resultdiff;
     111  ideal m;
     112  ideal m1=maxidealstep(i,startvar); // computes the corresponding ideal
     113  ideal m3,m2,J,resultdiff;
    113114  for(v=1;v<=k;v++)                   // consider the different t-coeff.
    114115   {
     
    118119     for(s=1;s<=k-v;s++)
    119120       {
    120          m1=maxidealstep(i,startvar); // computes the corresponding ideal
    121          m1=m1^s;
    122          u=size(m1);
     121         m3=m1^s;
     122         u=size(m3);
    123123         for(t=1;t<=u;t++)
    124124          {
    125             m2=contract(m1[t],m);     // actual differentiation
     125            m2=contract(m3[t],m);     // actual differentiation
    126126            J=J,m2;
    127127          }
     
    312312  ideal m;
    313313  ideal m1=maxidealstep(i,startvar);
    314   ideal m2,J,g_resultdiff;
     314  ideal m3,m2,J,g_resultdiff;
    315315  for(v=1;v<=k;v++)                   // consider the different t-coeff.
    316316    {
     
    320320     for(s=1;s<=k-v+1;s++)           // remark: "s<=k-v+1" special!for G!!!
    321321       {
    322          m1=m1^s;
    323          u=size(m1);
     322         m3=m1^s;
     323         u=size(m3);
    324324         for(t=1;t<=u;t++)
    325325          {
    326             m2=contract(m1[t],m);     // actual differentiation
     326            m2=contract(m3[t],m);     // actual differentiation
    327327            J=J,m2;
    328328          }
Note: See TracChangeset for help on using the changeset viewer.