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

Diff for /OpenXM_contrib2/asir2000/engine/nd.c between version 1.140 and 1.142

version 1.140, 2006/06/12 21:51:33 version 1.142, 2006/08/09 02:43:38
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.139 2006/06/12 00:46:48 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.141 2006/06/17 10:12:06 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 1548  int do_diagonalize(int sugar,int m)
Line 1548  int do_diagonalize(int sugar,int m)
         Q q,num,den;          Q q,num,den;
         union oNDC dn;          union oNDC dn;
   
         for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) {          for ( i = nd_psn-1; i >= 0; i-- ) {
                   if ( SG(nd_psh[i]) != sugar ) continue;
                 if ( nd_demand )                  if ( nd_demand )
                         nfv = ndv_load(i);                          nfv = ndv_load(i);
                 else                  else
Line 1641  again:
Line 1642  again:
                         }                          }
                         nfv = ndtondv(m,nf); nd_free(nf);                          nfv = ndtondv(m,nf); nd_free(nf);
                         nh = ndv_newps(m,nfv,0);                          nh = ndv_newps(m,nfv,0);
                         if ( ishomo && ++diag_count == diag_period ) {                          if ( !m && (ishomo && ++diag_count == diag_period) ) {
                                 diag_count = 0;                                  diag_count = 0;
                                 stat = do_diagonalize(sugar,m);                                  stat = do_diagonalize(sugar,m);
                                 if ( !stat ) {                                  if ( !stat ) {
Line 1676  int do_diagonalize_trace(int sugar,int m)
Line 1677  int do_diagonalize_trace(int sugar,int m)
         Q q,den,num;          Q q,den,num;
         union oNDC dn;          union oNDC dn;
   
         for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) {          for ( i = nd_psn-1; i >= 0; i-- ) {
                   if ( SG(nd_psh[i]) != sugar ) continue;
                 /* for nd_ps */                  /* for nd_ps */
                 s = ndvtond(m,nd_ps[i]);                  s = ndvtond(m,nd_ps[i]);
                 s = nd_separate_head(s,&head);                  s = nd_separate_head(s,&head);

Legend:
Removed from v.1.140  
changed lines
  Added in v.1.142

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>