Changeset 1156 in Sophya for trunk/SophyaLib/TArray/matharr.cc
- Timestamp:
- Aug 29, 2000, 6:10:32 PM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/TArray/matharr.cc
r1099 r1156 29 29 throw RangeCheckError("MathArray<T>::ApplyFunctionInPlace(TArray<T> & a..) Not Allocated Array a !"); 30 30 T * pe; 31 uint_8j,k;31 sa_size_t j,k; 32 32 if (a.AvgStep() > 0) { // regularly spaced elements 33 uint_8step = a.AvgStep();34 uint_8maxx = a.Size()*step;33 sa_size_t step = a.AvgStep(); 34 sa_size_t maxx = a.Size()*step; 35 35 pe = a.Data(); 36 36 for(k=0; k<maxx; k+=step ) pe[k] = (T)(f((double)pe[k])); 37 37 } 38 38 else { // Non regular data spacing ... 39 uint_4 ka = a.MaxSizeKA();40 uint_8step = a.Step(ka);41 uint_8gpas = a.Size(ka)*step;42 uint_8naxa = a.Size()/a.Size(ka);39 int_4 ka = a.MaxSizeKA(); 40 sa_size_t step = a.Step(ka); 41 sa_size_t gpas = a.Size(ka)*step; 42 sa_size_t naxa = a.Size()/a.Size(ka); 43 43 for(j=0; j<naxa; j++) { 44 44 pe = a.DataBlock().Begin()+a.Offset(ka,j); … … 61 61 throw RangeCheckError("MathArray<T>::ApplyFunctionInPlace(TArray<T> & a..) Not Allocated Array a !"); 62 62 T * pe; 63 uint_8j,k;63 sa_size_t j,k; 64 64 if (a.AvgStep() > 0) { // regularly spaced elements 65 uint_8step = a.AvgStep();66 uint_8maxx = a.Size()*step;65 sa_size_t step = a.AvgStep(); 66 sa_size_t maxx = a.Size()*step; 67 67 pe = a.Data(); 68 68 for(k=0; k<maxx; k+=step ) pe[k] = (T)(f((float)pe[k])); 69 69 } 70 70 else { // Non regular data spacing ... 71 uint_4 ka = a.MaxSizeKA();72 uint_8step = a.Step(ka);73 uint_8gpas = a.Size(ka)*step;74 uint_8naxa = a.Size()/a.Size(ka);71 int_4 ka = a.MaxSizeKA(); 72 sa_size_t step = a.Step(ka); 73 sa_size_t gpas = a.Size(ka)*step; 74 sa_size_t naxa = a.Size()/a.Size(ka); 75 75 for(j=0; j<naxa; j++) { 76 76 pe = a.DataBlock().Begin()+a.Offset(ka,j); … … 119 119 throw RangeCheckError("MathArray<T>::MeanSigma(TArray<T> const & a..) Not Allocated Array a !"); 120 120 const T * pe; 121 uint_8j,k;121 sa_size_t j,k; 122 122 mean=0.; 123 123 sig = 0.; 124 124 double valok; 125 125 if (a.AvgStep() > 0) { // regularly spaced elements 126 uint_8step = a.AvgStep();127 uint_8maxx = a.Size()*step;126 sa_size_t step = a.AvgStep(); 127 sa_size_t maxx = a.Size()*step; 128 128 pe = a.Data(); 129 129 for(k=0; k<maxx; k+=step ) { … … 133 133 } 134 134 else { // Non regular data spacing ... 135 uint_4 ka = a.MaxSizeKA();136 uint_8step = a.Step(ka);137 uint_8gpas = a.Size(ka)*step;138 uint_8naxa = a.Size()/a.Size(ka);135 int_4 ka = a.MaxSizeKA(); 136 sa_size_t step = a.Step(ka); 137 sa_size_t gpas = a.Size(ka)*step; 138 sa_size_t naxa = a.Size()/a.Size(ka); 139 139 for(j=0; j<naxa; j++) { 140 140 pe = a.DataBlock().Begin()+a.Offset(ka,j);
Note:
See TracChangeset
for help on using the changeset viewer.