Changeset 81
 Timestamp:
 Feb 16, 2010 4:57:12 PM (9 years ago)
 Location:
 trunk/ADOLC/src
 Files:

 3 edited
Legend:
 Unmodified
 Added
 Removed

trunk/ADOLC/src/convolut.c
r42 r81 78 78 /**/ 79 79 /* olvo 980616 nf */ 80 /* Increments truncated convolution of a and b to c and sets a to zero*/80 /* Increments truncated convolution of a and b to c */ 81 81 void inconv1( int dim, revreal *a, revreal *b, revreal *c ) { 82 82 revreal tmpVal; … … 84 84 for (i=dim1; i>=0; i) { 85 85 tmpVal = a[i]*b[0]; 86 a[i] = 0;87 86 for (j=1; j<=i; j++) 88 87 tmpVal += a[ij]*b[j]; … … 129 128 for (i=dim1; i>=0; i) { 130 129 tmpVal = a[i]*b[0]; 130 for (j=1; j<=i; j++) 131 tmpVal += a[ij]*b[j]; 132 c[i] = tmpVal; 133 } 134 } 135 136 /**/ 137 /* Decrements truncated convolution of a and b to c and sets a to zero */ 138 void deconvZeroR( int dim, revreal *a, revreal *b, revreal *c ) { 139 double tmpVal; 140 int i,j; 141 for (i=dim1; i>=0; i) { 142 tmpVal = a[i]*b[0]; 143 a[i] = 0; 131 144 for (j=1; j<=i; j++) 132 145 tmpVal += a[ij]*b[j]; 
trunk/ADOLC/src/convolut.h
r42 r81 53 53 void deconv0( int dim, double* a, revreal *b, double* c ); 54 54 void deconv1( int dim, revreal* a, revreal *b, revreal* c ); 55 void deconvZeroR( int dim, revreal *a, revreal *b, revreal *c ); 55 56 56 57 
trunk/ADOLC/src/ho_rev.c
r51 r81 1092 1092 AARG2_INC_O; 1093 1093 1094 deconv 0(k,Ares,Targ1,Aarg2);1095 deconv 0(k,Ares,Targ2,Aarg1);1094 deconv1(k,Ares,Targ1,Aarg2); 1095 deconv1(k,Ares,Targ2,Aarg1); 1096 1096 1097 1097 HOV_INC(Ares, k) … … 1228 1228 Tres, dp_Atemp); 1229 1229 VEC_COMPUTED_END 1230 deconv (k,Ares,dp_Atemp,Aarg);1230 deconvZeroR(k,Ares,dp_Atemp,Aarg); 1231 1231 1232 1232 HOV_INC(Ares, k) … … 2060 2060 } /* endswitch */ 2061 2061 2062 /* Get the next operation */ 2062 /* Get the next operation */ 2063 2063 operation=get_op_r(); 2064 2064 #if defined(ADOLC_DEBUG)
Note: See TracChangeset
for help on using the changeset viewer.