Changeset 2bbd4a in git


Ignore:
Timestamp:
May 19, 2014, 4:40:21 PM (10 years ago)
Author:
Anne Frühbis-Krüger <anne@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'fc741b6502fd8a97288eaa3eba6e5220f3c3df87')
Children:
46c7c3543716f62ef5bb7e5989d0c0876634829a
Parents:
4bc25c236b0d155b3b1e4ab40855eacaaa27c3eb
Message:
fixed corner case (minimal resolution very short) in AnnExt_R
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/primdec.lib

    r4bc25c r2bbd4a  
    39953995  if(n<nvars(basering))
    39963996  {
    3997      matrix f=transpose(re[n+1]);      //Hom(_,R)
    3998      module k=nres(f,2)[2];            //the kernel
    3999      matrix g=transpose(re[n]);        //the image of Hom(_,R)
    4000 
    4001      ideal ann=quotient1(g,k);           //the anihilator
    4002   }
    4003   else
    4004   {
    4005      ideal ann=Ann(transpose(re[n]));
    4006   }
     3997    if(size(re[n+1])>0)
     3998    {
     3999        matrix f=transpose(re[n+1]);      //Hom(_,R)
     4000        module k=nres(f,2)[2];            //the kernel
     4001        matrix g=transpose(re[n]);        //the image of Hom(_,R)
     4002
     4003        ideal ann=quotient1(g,k);         //the anihilator
     4004        return(ann);
     4005    }
     4006  }
     4007  // remaining case: re[n+1] is trivial
     4008  // either n is at least number of variables or
     4009  // resolution happens to be shorter
     4010  ideal ann=Ann(transpose(re[n]));
    40074011  return(ann);
    40084012}
Note: See TracChangeset for help on using the changeset viewer.