Changeset 93085a in git


Ignore:
Timestamp:
Mar 20, 2001, 3:18:03 PM (22 years ago)
Author:
Mathias Schulze <mschulze@…>
Branches:
(u'spielwiese', '8d54773d6c9e2f1d2593a28bc68b7eeab54ed529')
Children:
8960eca892efa27138b5030e016a1aede9a12bea
Parents:
8c4269a2193d8ebf4d512ac0304fe2a3644f8533
Message:
*mschulze: moved procedures for spectra to gaussman.lib


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

Legend:

Unmodified
Added
Removed

  • Singular/LIB/spectrum.lib

    r8c4269a r93085a  
    11///////////////////////////////////////////////////////////////////////////////
    2 version="$Id: spectrum.lib,v 1.12 2001-02-02 16:21:37 mschulze Exp $";
     2version="$Id: spectrum.lib,v 1.13 2001-03-20 14:18:03 mschulze Exp $";
    33category="Singularities";
    44info="
    … …  
    77
    88PROCEDURES:
    9  spectrumnd(poly[,1]);   spectrum of a nondegenerate isolated singularity
    10  semicont(s1,s2[,1]);    tests if s2 is semicontinous for s1
    11  semicontqh(s1,s2);      semicontinuity test using open and half open intervals
    12  spadd(s1,s2);           sum of two spectra s1 and s2
    13  spmul(s,k);             product of spectrum s with integer k
     9 spectrumnd(poly f[,1]);  spectrum of nondegenerate isolated singularity f
    1410";
    1511
    1612///////////////////////////////////////////////////////////////////////////////
    1713
    18 proc spectrumnd (poly f, list #)
    19 "USAGE:   spectrumnd(f[,1]);  f poly
    20 ASSUME:  f has nondegenerate principal part
    21 RETURN:  a list containing the spectrum of f
    22 NOTE     if a second argument 1 is given,
    23          no test for a degenerate principal part will be done
    24 EXAMPLE: example spectrumnd; shows examples
     14proc spectrumnd (poly f,list #)
     15"USAGE:    spectrumnd(f[,1]); poly f
     16ASSUME:   basering has local ordering,
     17          f has isolated singularity at 0 and nondegenerate principal part
     18RETURN:
     19@format
     20list S:
     21  ideal S[1]: spectral numbers
     22  intvec S[2]:
     23    int S[2][i]: multiplicity of spectral number S[1][i]
     24@end format
     25NOTE:     if a second argument 1 is given,
     26          no test for a degenerate principal part will be done
     27SEE_ALSO: gaussman_lib
     28KEYWORDS: singularities; Gauss-Manin connection; spectrum
     29EXAMPLE:  example spectrumnd; shows an example
    2530"
    2631{
    2732  if(size(#)==0)
    28   {  return(system("spectrum",f)); }
    29   return (system("spectrum",f,#[1]));
     33  {
     34    list l=system("spectrum",f);
     35  }
     36  else
     37  {
     38    list l=system("spectrum",f,#[1]);
     39  }
     40  ideal s=number(l[4][1])/l[5][1]-1;
     41  int i;
     42  for(i=l[3];i>1;i--)
     43  {
     44    s[i]=number(l[4][i])/l[5][i]-1;
     45  }
     46  return(list(s,l[6]));
    3047}
    3148example
    … …  
    3754///////////////////////////////////////////////////////////////////////////////
    3855
    39 proc semicont (list s, list ss, list #)
    40 "USAGE:   semicont(s,ss[,1]);  s, ss list
    41 RETURN:  1 if ss is semicontinous for s using half open intervals, 0 otherwise
    42 NOTE:    if a third argument 1 is given, open intervals are used
    43 EXAMPLE: example semicont; shows examples
    44 "
     56proc semicont(list #)
    4557{
    46   if(size(#)==0)
    47   {  return(system("semic",s,ss)); }
    48   return (system("semic",s,ss,#[1]));
    49 }
    50 example
    51 { "EXAMPLE:"; echo = 2;
     58  ERROR("semicont was replaced by spsemicont in gaussman.lib");
    5259}
    5360///////////////////////////////////////////////////////////////////////////////
    54 proc semicontsqh (list s, list ss)
    55 "USAGE: semicontsqh(s,ss);  s, ss list
    56 RETURN: 1 if ss is semicontinous for s using open and half open intervals,
    57         0 otherwise
    58 EXAMPLE: example semicontsqh; shows examples
    59 "
     61
     62proc semicontqh(list #)
    6063{
    61   return (system("semic",s,ss,1));
    62 }
    63 example
    64 { "EXAMPLE:"; echo = 2;
     64  ERROR("semicontqh was replaced by spsemicont in gaussman.lib");
    6565}
    6666///////////////////////////////////////////////////////////////////////////////
    67 proc spadd (list s1, list s2)
    68 "USAGE:  spadd(s1,s2);  s1, s2 list
    69 RETURN:  a list containing the sum of the two spectra s1 and s2
    70 EXAMPLE: example spadd; shows examples
    71 "
    72 {
    73   return (system("spadd",s1,s2));
    74 }
    75 example
    76 { "EXAMPLE:"; echo = 2;
    77 }
    78 ///////////////////////////////////////////////////////////////////////////////
    79 proc spmul (list s, int k)
    80 "USAGE:  spmul(s,k);  s list, k int
    81 RETURN:  a list containing the product of the spectrum s with the integer k
    82 EXAMPLE: example spmul; shows examples
    83 "
    84 {
    85   return (system("spmul",s,k));
    86 }
    87 example
    88 { "EXAMPLE:"; echo = 2;
    89 }
    90 ///////////////////////////////////////////////////////////////////////////////
     67
Note: See TracChangeset for help on using the changeset viewer.