/[ascend]/trunk/base/generic/solver/linsol.c
ViewVC logotype

Diff of /trunk/base/generic/solver/linsol.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 707 by ben.allan, Tue Apr 4 07:18:41 2006 UTC revision 708 by johnpye, Tue Jun 27 07:34:31 2006 UTC
# Line 818  static void invert(linsol_system_t sys) Line 818  static void invert(linsol_system_t sys)
818    
819     length = sys->rng.high - sys->rng.low + 1;     length = sys->rng.high - sys->rng.low + 1;
820     tmp_array_origin = length > 0 ?     tmp_array_origin = length > 0 ?
821        (real64 *)ascmalloc( length*sizeof(real64) ) : NULL;        ASC_NEW_ARRAY(real64,length) : NULL;
822     tmp = tmp_array_origin != NULL ?     tmp = tmp_array_origin != NULL ?
823        tmp_array_origin - sys->rng.low : NULL;        tmp_array_origin - sys->rng.low : NULL;
824    
# Line 879  static real64 *raise_capacity(real64 *ve Line 879  static real64 *raise_capacity(real64 *ve
879  /*  /*
880     real64 *newvec;     real64 *newvec;
881     newvec = newcap > 0 ?     newvec = newcap > 0 ?
882        (real64 *)ascmalloc( newcap * sizeof(real64) ) : NULL;        ASC_NEW_ARRAY(real64,newcap ) : NULL;
883     if( vec != NULL ) {     if( vec != NULL ) {
884        mem_move_cast(vec,newvec,oldcap*sizeof(real64));        mem_move_cast(vec,newvec,oldcap*sizeof(real64));
885        ascfree( (POINTER)vec );        ascfree( (POINTER)vec );
# Line 892  static real64 *raise_capacity(real64 *ve Line 892  static real64 *raise_capacity(real64 *ve
892      newvec = (real64 *)ascrealloc(vec,(newcap * sizeof(real64)));      newvec = (real64 *)ascrealloc(vec,(newcap * sizeof(real64)));
893    else    else
894      newvec=newcap > 0 ?      newvec=newcap > 0 ?
895        (real64 *)ascmalloc( newcap * sizeof(real64) ) : NULL;        ASC_NEW_ARRAY(real64,newcap ) : NULL;
896    return newvec;    return newvec;
897  }  }
898    
# Line 1180  static void invert_alt(linsol_system_t s Line 1180  static void invert_alt(linsol_system_t s
1180    
1181     length = sys->rng.high - sys->rng.low + 1;     length = sys->rng.high - sys->rng.low + 1;
1182     tmp = length > 0 ?     tmp = length > 0 ?
1183        (real64 *)ascmalloc( length*sizeof(real64) ) : NULL;        ASC_NEW_ARRAY(real64,length) : NULL;
1184    
1185     sys->smallest_pivot = MAXDOUBLE;     sys->smallest_pivot = MAXDOUBLE;
1186     candidates.row.high = sys->rng.high;     candidates.row.high = sys->rng.high;

Legend:
Removed from v.707  
changed lines
  Added in v.708

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22