Changeset a23a7fd in git for Singular/LIB/gaussman.lib
 Timestamp:
 Dec 6, 2000, 3:56:26 PM (23 years ago)
 Branches:
 (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
 Children:
 ac31058626d6128287fd40e307e34fc22c7c0beb
 Parents:
 a79bc624a96888a1bce8c52f3c8f6a33640de235
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

Singular/LIB/gaussman.lib
ra79bc62 ra23a7fd 1 1 /////////////////////////////////////////////////////////////////////////////// 2 2 3 version="$Id: gaussman.lib,v 1.1 3 20001206 12:34:50mschulze Exp $";3 version="$Id: gaussman.lib,v 1.14 20001206 14:56:26 mschulze Exp $"; 4 4 info=" 5 5 LIBRARY: gaussman.lib GAUSSMANIN CONNECTION OF A SINGULARITY … … 820 820 } 821 821 822 dbprint(printlevelvoice+2, 823 "//gaussman::vjacob: compute multiplication of Jacobian algebra"); 822 824 list M; 823 825 for(i=ncols(m);i>=1;i) … … 836 838 while(r0<number(s[ncols(s)]s[1])) 837 839 { 840 dbprint(printlevelvoice+2, 841 "//gaussman::vjacob: compute next difference r0 of spectral numbers"); 838 842 r1=number(s[ncols(s)]s[1]); 839 843 for(j=ncols(s);j>=1;j) … … 856 860 } 857 861 r0=r1; 862 dbprint(printlevelvoice+2,"//gaussman::vjacob: r0="+string(r1)); 858 863 859 864 l=ideal(); … … 863 868 } 864 869 865 j0=mud[ncols(s)]+1; 866 for(j=ncols(s);j>=1;j) 870 dbprint(printlevelvoice+2,"//gaussman::vjacob: compute conditions"); 871 j=ncols(s); 872 j0=mu; 873 while(j>=1) 867 874 { 868 875 i0=1; … … 883 890 if(i0>1) 884 891 { 885 l[k]=l[k],M0[1..i01,j0 ..j0+d[j]1];892 l[k]=l[k],M0[1..i01,j0d[j]+1..j0]; 886 893 } 887 894 } 888 895 j0=j0d[j]; 889 } 890 896 j; 897 } 898 899 dbprint(printlevelvoice+2,"//gaussman::vjacob: compute condition matrix"); 891 900 L=transpose(module(l[1])); 892 901 for(k=2;k<=ncols(m);k++) … … 895 904 } 896 905 906 dbprint(printlevelvoice+2, 907 "//gaussman::vjacob: compute kernel of condition matrix"); 897 908 v0=v0+list(syz(L)); 898 909 s0=s0,r0; 899 910 } 900 911 912 dbprint(printlevelvoice+2,"//gaussman::vjacob: compute graded parts"); 901 913 option(redSB); 902 914 for(i=1;i<size(v0);i++) … … 906 918 } 907 919 920 dbprint(printlevelvoice+2, 921 "//gaussman::vjacob: remove trivial graded parts"); 908 922 i=1; 909 923 while(size(v0[i])==0)
Note: See TracChangeset
for help on using the changeset viewer.