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

Diff for /OpenXM_contrib2/asir2000/engine/dalg.c between version 1.10 and 1.11

version 1.10, 2005/08/02 07:16:42 version 1.11, 2005/08/24 06:28:39
Line 1 
Line 1 
 /*  /*
  * $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.9 2005/07/11 00:24:02 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/engine/dalg.c,v 1.10 2005/08/02 07:16:42 noro Exp $
 */  */
   
 #include "ca.h"  #include "ca.h"
Line 675  NODE inv_or_split_dalg(DAlg a,DAlg *c)
Line 675  NODE inv_or_split_dalg(DAlg a,DAlg *c)
                                 break;                                  break;
                 if ( j >= 0 ) {                  if ( j >= 0 ) {
                         dp_subd(m,mb[j],&d);                          dp_subd(m,mb[j],&d);
                         muld(CO,d,simp[j]->nm,&u);                          if ( simp[j] ) {
                         MKDAlg(u,simp[j]->dn,t);                                  muld(CO,d,simp[j]->nm,&u);
                         simpdalg(t,&simp[i]);                                  MKDAlg(u,simp[j]->dn,t);
                                   simpdalg(t,&simp[i]);
                           } else
                                   simp[i] = 0;
                 } else {                  } else {
                         MKDAlg(m,ONE,t);                          MKDAlg(m,ONE,t);
                         muldalg(t,a0,&simp[i]);                          muldalg(t,a0,&simp[i]);
                 }                  }
                 gcdn(NM(simp[i]->dn),ln,&gn); divsn(ln,gn,&qn);                  if ( simp[i] ) {
                 muln(NM(simp[i]->dn),qn,&ln);                          gcdn(NM(simp[i]->dn),ln,&gn); divsn(ln,gn,&qn);
                           muln(NM(simp[i]->dn),qn,&ln);
                   }
         }          }
         initd(current_spec);          initd(current_spec);
         NTOQ(ln,1,dn);          NTOQ(ln,1,dn);
Line 691  NODE inv_or_split_dalg(DAlg a,DAlg *c)
Line 696  NODE inv_or_split_dalg(DAlg a,DAlg *c)
         mat = (Q **)BDY(mobj);          mat = (Q **)BDY(mobj);
         mulq(dn,a->dn,&mat[0][dim]);          mulq(dn,a->dn,&mat[0][dim]);
         for ( j = 0; j < dim; j++ ) {          for ( j = 0; j < dim; j++ ) {
                 divq(dn,simp[j]->dn,&mul);                  if ( simp[j] ) {
                 for ( i = dim-1, mp = BDY(simp[j]->nm); mp && i >= 0; i-- )                          divq(dn,simp[j]->dn,&mul);
                         if ( dl_equal(n,BDY(mb[i])->dl,mp->dl) ) {                          for ( i = dim-1, mp = BDY(simp[j]->nm); mp && i >= 0; i-- )
                                 mulq(mul,(Q)mp->c,&mat[i][j]);                                  if ( dl_equal(n,BDY(mb[i])->dl,mp->dl) ) {
                                 mp = NEXT(mp);                                          mulq(mul,(Q)mp->c,&mat[i][j]);
                         }                                          mp = NEXT(mp);
                                   }
                   }
         }          }
         get_eg(&eg0);          get_eg(&eg0);
         rank = generic_gauss_elim_hensel(mobj,&sol,&dnsol,&rinfo,&cinfo);          rank = generic_gauss_elim_hensel(mobj,&sol,&dnsol,&rinfo,&cinfo);
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>