[BACK]Return to sp CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / lib

Diff for /OpenXM_contrib2/asir2000/lib/sp between version 1.16 and 1.17

version 1.16, 2010/07/14 04:48:14 version 1.17, 2016/02/04 04:17:21
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/lib/sp,v 1.15 2006/06/23 08:57:47 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/lib/sp,v 1.16 2010/07/14 04:48:14 noro Exp $
 */  */
 /*  /*
         sp : functions related to algebraic number fields          sp : functions related to algebraic number fields
Line 211  def flatmf(L) {
Line 211  def flatmf(L) {
 def af(P,AL)  def af(P,AL)
 {  {
         RESTIME=UFTIME=GCDTIME=SQTIME=0;          RESTIME=UFTIME=GCDTIME=SQTIME=0;
       V = var(P);
       LC = coef(P,deg(P,V),V);
       if ( ntype(LC) != 1 )
         P = simpalg(1/LC*P);
         S = reverse(asq(P));          S = reverse(asq(P));
         for ( L = []; S != []; S = cdr(S) ) {          for ( L = []; S != []; S = cdr(S) ) {
                 FM = car(S); F = FM[0]; M = FM[1];                  FM = car(S); F = FM[0]; M = FM[1];
Line 227  def af_sp(P,AL,HINT)
Line 231  def af_sp(P,AL,HINT)
 {  {
         if ( !P || type(P) == NUM )          if ( !P || type(P) == NUM )
                 return [P];                  return [P];
       V = var(P);
       LC = coef(P,deg(P,V),V);
       if ( ntype(LC) != 1 )
         P = simpalg(1/LC*P);
         P1 = simpcoef(simpalg(P));          P1 = simpcoef(simpalg(P));
         return af_spmain(P1,AL,1,HINT,P1,[]);          return af_spmain(P1,AL,1,HINT,P1,[]);
 }  }
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>