Changeset 034ce1 in git for Singular/LIB/finvar.lib


Ignore:
Timestamp:
May 12, 2000, 2:25:47 PM (24 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
7bb71fad3a9d45a2e30c8458373b89a01eecf500
Parents:
73d7b5cf051d2bb349542039ed3ac2b513f7c109
Message:
* hannes: execute


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

Legend:

Unmodified
Added
Removed

  • Singular/LIB/finvar.lib

    r73d7b5 r034ce1  
    1 // $Id: finvar.lib,v 1.26 2000-05-02 14:52:48 Singular Exp $
     1// $Id: finvar.lib,v 1.27 2000-05-12 12:25:43 Singular Exp $
    22// author: Agnes Eileen Heydtmann, email:agnes@math.uni-sb.de
    33// last change: 98/11/05
    44///////////////////////////////////////////////////////////////////////////////
    5 version="$Id: finvar.lib,v 1.26 2000-05-02 14:52:48 Singular Exp $"
     5version="$Id: finvar.lib,v 1.27 2000-05-12 12:25:43 Singular Exp $"
    66info="
    77LIBRARY:  finvar.lib       LIBRARY TO CALCULATE INVARIANT RINGS OF FINITE GROUPS
    … …  
    765765    if (minpoly==0)
    766766    { if (i>size(chst))
    767       { execute "ring "+newring+"=0,("+varstr(br)+"),("+ordstr(br)+")";
     767      { execute("ring "+newring+"=0,("+varstr(br)+"),("+ordstr(br)+")");
    768768      }
    769769      else
    770770      { chst=chst[i..size(chst)];
    771         execute "ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")";
     771        execute
     772        ("ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")");
    772773      }
    773774    }
    … …  
    776777      minp=minp[2..size(minp)-1];
    777778      chst=chst[i..size(chst)];
    778       execute "ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")";
    779       execute "minpoly="+minp;
     779      execute("ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")");
     780      execute("minpoly="+minp);
    780781    }
    781782    matrix I=diag(1,n);
    … …  
    813814      }
    814815      setring `newring`;
    815       execute "matrix G(i)["+string(n)+"]["+string(n)+"]="+stM(i);
     816      execute("matrix G(i)["+string(n)+"]["+string(n)+"]="+stM(i));
    816817      p=det(I-v1*G(i));                // denominator of new term -
    817818      M[1,1]=M[1,1]*p+M[1,2];          // expanding M[1,1]/M[1,2] + 1/p
    … …  
    12021203    if (minpoly==0)
    12031204    { if (i>size(chst))
    1204       { execute "ring "+newring+"=0,("+varstr(br)+"),("+ordstr(br)+")";
     1205      { execute("ring "+newring+"=0,("+varstr(br)+"),("+ordstr(br)+")");
    12051206      }
    12061207      else
    12071208      { chst=chst[i..size(chst)];
    1208         execute "ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")";
     1209        execute
     1210        ("ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")");
    12091211      }
    12101212    }
    … …  
    12131215      minp=minp[2..size(minp)-1];
    12141216      chst=chst[i..size(chst)];
    1215       execute "ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")";
    1216       execute "minpoly="+minp;
     1217      execute("ring "+newring+"=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")");
     1218      execute("minpoly="+minp);
    12171219    }
    12181220    poly v1=var(1);                    // the Molien series will be in terms of
    … …  
    12341236    }
    12351237    setring `newring`;
    1236     execute "matrix G(1)["+string(n)+"]["+string(n)+"]="+stM(1);
     1238    execute("matrix G(1)["+string(n)+"]["+string(n)+"]="+stM(1));
    12371239    matrix A(2)[1][2];                 // A(2) will contain the Molien series -
    12381240    A(2)[1,1]=1;                       // A(2)[1,1] will be the numerator
    … …  
    12711273        }
    12721274        setring `newring`;
    1273         execute "matrix G(i)["+string(n)+"]["+string(n)+"]="+stM(i);
     1275        execute("matrix G(i)["+string(n)+"]["+string(n)+"]="+stM(i));
    12741276        p=det(I-v1*G(i));              // denominator of new term -
    12751277        A(2)[1,1]=A(2)[1,1]*p+A(2)[1,2]; // expanding A(2)[1,1]/A(2)[1,2] +1/p
    … …  
    13301332          }
    13311333          setring `newring`;
    1332           execute "matrix G(g)["+string(n)+"]["+string(n)+"]="+stM(g);
     1334          execute("matrix G(g)["+string(n)+"]["+string(n)+"]="+stM(g));
    13331335          p=det(I-v1*G(g));            // denominator of new term
    13341336          A(2)[1,1]=A(2)[1,1]*p+A(2)[1,2];
    … …  
    14731475    poly A(1)=M[1,2];                  // denominator of Molien series (for now)
    14741476    string mp=string(minpoly);
    1475     execute "ring R=("+charstr(br)+"),("+varstr(br)+"),ds;";
    1476     execute "minpoly=number("+mp+");";
     1477    execute("ring R=("+charstr(br)+"),("+varstr(br)+"),ds;");
     1478    execute("minpoly=number("+mp+");");
    14771479    poly A(1)=0;                       // A(1) will contain the sum of n terms -
    14781480    poly min;                          // min will be our smallest term -
    … …  
    15961598  }
    15971599 //----------------------------------------------------------------------------
    1598   execute "ring T=("+charstr(br)+"),("+varstr(br)+",p(1..m)),lp;";
     1600  execute("ring T=("+charstr(br)+"),("+varstr(br)+",p(1..m)),lp;");
    15991601  // p(1..m) are the general coefficients of the general polynomial of degree g
    1600   execute "ideal vars="+varstr(br)+";";
     1602  execute("ideal vars="+varstr(br)+";");
    16011603  map f;
    16021604  ideal mon=imap(br,mon);
    … …  
    30453047          }
    30463048          if (not(flag))
    3047           { execute "test_dim="+string(answer[1..size(answer)]);
     3049          { execute("test_dim="+string(answer[1..size(answer)]));
    30483050            if (test_dim<=max)
    30493051            { flag=1;
    … …  
    31543156            }
    31553157            if (not(flag))
    3156             { execute "test_dim="+string(answer[1..size(answer)]);
     3158            { execute("test_dim="+string(answer[1..size(answer)]));
    31573159              if (test_dim<=max)
    31583160              { flag=1;
    … …  
    51165118  string mp=string(minpoly);           // generating a ring where we can do
    51175119                                       // elimination
    5118   execute "ring R=("+charstr(br)+"),(x(1..n),y(1..n),h),dp;";
    5119   execute "minpoly=number("+mp+");";
     5120  execute("ring R=("+charstr(br)+"),(x(1..n),y(1..n),h),dp;");
     5121  execute("minpoly=number("+mp+");");
    51205122  map f=br,maxideal(1);                // canonical mapping
    51215123  matrix M[k][m+k*n];
    … …  
    56405642    int m=ncols(F);
    56415643    string mp=string(minpoly);
    5642     execute "ring R=("+charstr(br)+"),("+varstr(br)+",y(1..m)),dp;";
    5643     execute "minpoly=number("+mp+");";
     5644    execute("ring R=("+charstr(br)+"),("+varstr(br)+",y(1..m)),dp;");
     5645    execute("minpoly=number("+mp+");");
    56445646    ideal I=ideal(imap(br,F));
    56455647    for (int i=1;i<=m;i++)
    … …  
    56475649    }
    56485650    I=elim(I,1,n);
    5649     execute "ring "+newring+"=("+charstr(br)+"),(y(1..m)),dp(m);";
    5650     execute "minpoly=number("+mp+");";
     5651    execute("ring "+newring+"=("+charstr(br)+"),(y(1..m)),dp(m);");
     5652    execute("minpoly=number("+mp+");");
    56515653    ideal vars;
    56525654    for (i=2;i<=n;i++)
    … …  
    56585660    kill emb, vars, R;
    56595661    keepring `newring`;
    5660     // execute "keepring "+newring+";";
    56615662    return();
    56625663  }
    … …  
    57015702    int m=ncols(F);
    57025703    string mp=string(minpoly);
    5703     execute "ring R=("+charstr(br)+"),("+varstr(br)+",y(1..m)),lp;";
    5704     execute "minpoly=number("+mp+");";
     5704    execute("ring R=("+charstr(br)+"),("+varstr(br)+",y(1..m)),lp;");
     5705    execute("minpoly=number("+mp+");");
    57055706    ideal J=ideal(imap(br,F));
    57065707    ideal I=imap(br,I);
    … …  
    57275728    }
    57285729    G=compress(G);
    5729     execute "ring "+newring+"=("+charstr(br)+"),(y(1..m)),lp;";
    5730     execute "minpoly=number("+mp+");";
     5730    execute("ring "+newring+"=("+charstr(br)+"),(y(1..m)),lp;");
     5731    execute("minpoly=number("+mp+");");
    57315732    ideal vars;
    57325733    for (i=2;i<=n;i++)
    … …  
    57385739    kill vars, emb;
    57395740    keepring `newring`;
    5740     // execute "keepring "+newring+";";
    57415741    return();
    57425742  }
    … …  
    57755775    string newring="E";
    57765776    relative_orbit_variety(I,F,newring);
    5777     execute "ring R=("+charstr(br)+"),("+varstr(br)+","+varstr(E)+"),lp;";
     5777    execute("ring R=("+charstr(br)+"),("+varstr(br)+","+varstr(E)+"),lp;");
    57785778    ideal F=imap(br,F);
    57795779    for (int i=1;i<=n;i++)
Note: See TracChangeset for help on using the changeset viewer.