[BACK]Return to algnum.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / builtin

Diff for /OpenXM_contrib2/asir2000/builtin/algnum.c between version 1.12 and 1.14

version 1.12, 2005/08/02 07:16:41 version 1.14, 2013/11/17 17:34:59
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/builtin/algnum.c,v 1.11 2005/07/11 00:24:02 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/algnum.c,v 1.13 2007/02/15 02:06:20 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 146  void Pdalgtodp(NODE arg,LIST *r)
Line 146  void Pdalgtodp(NODE arg,LIST *r)
   
 void Pdptodalg(NODE arg,DAlg *r)  void Pdptodalg(NODE arg,DAlg *r)
 {  {
         DP d;          DP d,nm,nm1;
           MP m;
           Q c,a;
         DAlg t;          DAlg t;
   
         d = (DP)ARG0(arg);          d = (DP)ARG0(arg);
         MKDAlg(d,ONE,t);          if ( !d ) *r = 0;
         simpdalg(t,r);          else {
                   for ( m = BDY(d); m; m = NEXT(m) )
                           if ( !INT((Q)m->c) ) break;
                   if ( !m ) {
                           MKDAlg(d,(Q)ONE,t);
                   } else {
                           dp_ptozp(d,&nm);
                           divq((Q)BDY(d)->c,(Q)BDY(nm)->c,&c);
                           NTOQ(NM(c),SGN(c),a);
                           muldc(CO,nm,(P)a,&nm1);
                           NTOQ(DN(c),1,a);
                           MKDAlg(nm1,a,t);
                   }
                   simpdalg(t,r);
           }
 }  }
   
 void Pdalgtoup(NODE arg,LIST *r)  void Pdalgtoup(NODE arg,LIST *r)
Line 203  void Pinv_or_split_dalg(NODE arg,Obj *rp)
Line 219  void Pinv_or_split_dalg(NODE arg,Obj *rp)
 {  {
         NODE gen,t,nd0,nd;          NODE gen,t,nd0,nd;
         LIST list;          LIST list;
         int l,i,j,k,n;          int l,i,j,n;
         DP *ps,*ps1,*psw;          DP *ps,*ps1,*psw;
         NumberField nf;          NumberField nf;
         DAlg inv;          DAlg inv;
Line 239  void Pinv_or_split_dalg(NODE arg,Obj *rp)
Line 255  void Pinv_or_split_dalg(NODE arg,Obj *rp)
                 for ( ; i >= 0; i-- ) {                  for ( ; i >= 0; i-- ) {
                         NEXTNODE(nd0,nd); BDY(nd) = (pointer)psw[i];                          NEXTNODE(nd0,nd); BDY(nd) = (pointer)psw[i];
                 }                  }
                 for ( ; t; t = NEXT(t), k++ ) {                  for ( ; t; t = NEXT(t) ) {
                         NEXTNODE(nd0,nd); BDY(nd) = BDY(t);                          NEXTNODE(nd0,nd); BDY(nd) = BDY(t);
                 }                  }
                 NEXT(nd) = 0;                  NEXT(nd) = 0;
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.14
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>