Changeset fb40dfd in git


Ignore:
Timestamp:
Jan 5, 2015, 4:26:36 PM (9 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
71e7cb88a8d50e24a3fe14ea10cfc63d1885daa0
Parents:
b8245236db9ce8035a9e4bbb28ef632e7ca8e4d8
Message:
optimized Primdec::AnnExt_R
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/primdec.lib

    rb824523 rfb40dfd  
    41244124
    41254125  if(n < 0)
    4126     {
    4127       ideal ann = ideal(1);
    4128       return(ann);
    4129     }
     4126  {
     4127    return(ideal(1));
     4128  }
    41304129  int l = size(re);
    41314130
    41324131  if(n < l)
    4133     {
    4134       matrix f = transpose(re[n+1]);
    4135       if(n == 0)
    4136         {
    4137           matrix g = 0*gen(ncols(f));
    4138         }
    4139       else
    4140         {
    4141           matrix g = transpose(re[n]);
    4142         }
    4143       module k = syz(f);
    4144       ideal ann = quotient1(g,k);
    4145       return(ann);
    4146     }
     4132  {
     4133    matrix f = transpose(re[n+1]);
     4134    if(n == 0)
     4135    {
     4136      matrix g = matrix(0,1,ncols(f));
     4137    }
     4138    else
     4139    {
     4140      matrix g = transpose(re[n]);
     4141    }
     4142    module k = syz(f);
     4143    return(quotient1(g,k));
     4144  }
    41474145
    41484146  if(n == l)
    4149     {
    4150       ideal ann = Ann(transpose(re[n]));
    4151       return(ann);
    4152     }
    4153 
    4154   ideal ann = ideal(1);
    4155   return(ann);
     4147  {
     4148    return(Ann(transpose(re[n])));
     4149  }
     4150
     4151  return(ideal(1));
    41564152}
    41574153///////////////////////////////////////////////////////////////////////////////
     
    70477043
    70487044  op@P = option(get);
    7049   execute("ring gnir = ("+charstr(basering)+"),("+varstr(basering)+"),(C,lp);");
     7045  def gnir=changeordTo(basering,"lp");
     7046  setring gnir;
    70507047
    70517048  op=option(get);
Note: See TracChangeset for help on using the changeset viewer.