Changeset a39a07 in git for Singular/LIB/primdec.lib


Ignore:
Timestamp:
May 13, 1998, 12:54:07 PM (26 years ago)
Author:
Olaf Bachmann <obachman@…>
Branches:
(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
Children:
edf0e9681a3d776622b3fefcdd261407d1a25780
Parents:
546de146b5808486dd0067fe60ba45a307acd96b
Message:
* changed varname lres to lres0 -- lres is now a command


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

Legend:

Unmodified
Added
Removed
  • Singular/LIB/primdec.lib

    r546de1 ra39a07  
    1 // $Id: primdec.lib,v 1.15 1998-05-13 07:38:35 obachman Exp $
     1// $Id: primdec.lib,v 1.16 1998-05-13 10:54:07 obachman Exp $
    22///////////////////////////////////////////////////////
    33// primdec.lib
     
    1111//////////////////////////////////////////////////////
    1212
    13 version="$Id: primdec.lib,v 1.15 1998-05-13 07:38:35 obachman Exp $";
     13version="$Id: primdec.lib,v 1.16 1998-05-13 10:54:07 obachman Exp $";
    1414info="
    1515LIBRARY: primdec.lib: PROCEDURE FOR PRIMARY DECOMPOSITION (I)
     
    844844  int nva = nvars(basering);
    845845  int @k,@s,@n,@k1,zz;
    846   list primary,lres,lres1,act,@lh,@wh;
     846  list primary,lres0,lres1,act,@lh,@wh;
    847847  map phi,psi,phi1,psi1;
    848848  ideal jmap,jmap1,jmap2,helpprim,@qh,@qht,ser1;
     
    10821082       
    10831083       
    1084        kill lres;
    1085        list lres;
     1084       kill lres0;
     1085       list lres0;
    10861086       if(size(@lh)==2)
    10871087       {
    10881088          helpprim=@lh[2];
    1089           lres[1]=primary[2*@k-1];
     1089          lres0[1]=primary[2*@k-1];
    10901090          ser1=psi(helpprim);
    1091           lres[2]=psi1(ser1);
    1092           if(size(reduce(lres[2],lres[1],1))==0)
     1091          lres0[2]=psi1(ser1);
     1092          if(size(reduce(lres0[2],lres0[1],1))==0)
    10931093          {
    10941094             primary[2*@k]=primary[2*@k-1];
     
    11051105                @f=act[1][@n]^act[2][@n];
    11061106                ser1=psi(@f);
    1107                 lres[2*@n-1]=interred(primary[2*@k-1]+psi1(ser1));
     1107                lres0[2*@n-1]=interred(primary[2*@k-1]+psi1(ser1));
    11081108                helpprim=@lh[2*@n];
    11091109                ser1=psi(helpprim);
    1110                 lres[2*@n]=psi1(ser1);
     1110                lres0[2*@n]=psi1(ser1);
    11111111             }
    11121112          }
     
    11141114          {
    11151115             lres1=psi(@lh);
    1116              lres=psi1(lres1);
     1116             lres0=psi1(lres1);
    11171117          }
    11181118       }
     
    11301130       execute(@ri);
    11311131       list @lvec;
    1132        list @lr=imap(@P,lres);
     1132       list @lr=imap(@P,lres0);
    11331133       ideal @lr1;
    11341134
     
    11961196       kill @lvec;
    11971197       setring @P;
    1198        lres=imap(@Phelp1,@lr);
     1198       lres0=imap(@Phelp1,@lr);
    11991199       kill @Phelp1;
    1200        for(@n=1;@n<=size(lres);@n++)
     1200       for(@n=1;@n<=size(lres0);@n++)
    12011201       {
    1202           lres[@n]=clearSB(lres[@n]);
    1203           attrib(lres[@n],"isSB",1);
     1202          lres0[@n]=clearSB(lres0[@n]);
     1203          attrib(lres0[@n],"isSB",1);
    12041204       }
    12051205
    1206        primary[2*@k-1]=lres[1];
    1207        primary[2*@k]=lres[2];
     1206       primary[2*@k-1]=lres0[1];
     1207       primary[2*@k]=lres0[2];
    12081208       @s=size(primary)/2;
    1209        for(@n=1;@n<=size(lres)/2-1;@n++)
     1209       for(@n=1;@n<=size(lres0)/2-1;@n++)
    12101210       {
    1211          primary[2*@s+2*@n-1]=lres[2*@n+1];
    1212          primary[2*@s+2*@n]=lres[2*@n+2];
     1211         primary[2*@s+2*@n-1]=lres0[2*@n+1];
     1212         primary[2*@s+2*@n]=lres0[2*@n+2];
    12131213       }
    12141214       @k--;
Note: See TracChangeset for help on using the changeset viewer.