Index: trunk/psLib/src/dataManip/psMatrixVectorArithmetic.c =================================================================== --- trunk/psLib/src/dataManip/psMatrixVectorArithmetic.c (revision 1647) +++ trunk/psLib/src/dataManip/psMatrixVectorArithmetic.c (revision 1648) @@ -30,6 +30,6 @@ * @author Ross Harman, MHPCC * - * @version $Revision: 1.26 $ $Name: not supported by cvs2svn $ - * @date $Date: 2004-08-27 23:29:36 $ + * @version $Revision: 1.27 $ $Name: not supported by cvs2svn $ + * @date $Date: 2004-08-27 23:33:31 $ * * Copyright 2004 Maui High Performance Computing Center, University of Hawaii @@ -114,7 +114,7 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - o = &((psScalar* )OUT)->data.TYPE; \ - i1 = &((psScalar* )IN1)->data.TYPE; \ - i2 = &((psScalar* )IN2)->data.TYPE; \ + o = &((psScalar* )OUT)->data.TYPE; \ + i1 = &((psScalar* )IN1)->data.TYPE; \ + i2 = &((psScalar* )IN2)->data.TYPE; \ *o = OP; \ } @@ -127,8 +127,8 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - npt = ((psVector* )IN1)->n; \ - o = ((psVector* )OUT)->data.TYPE; \ - i1 = &((psScalar* )IN1)->data.TYPE; \ - i2 = ((psVector* )IN2)->data.TYPE; \ + npt = ((psVector* )IN1)->n; \ + o = ((psVector* )OUT)->data.TYPE; \ + i1 = &((psScalar* )IN1)->data.TYPE; \ + i2 = ((psVector* )IN2)->data.TYPE; \ for (i=0; i < npt; i++, o++, i2++) { \ *o = OP; \ @@ -145,10 +145,10 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - numRows = ((psImage* )IN1)->numRows; \ - numCols = ((psImage* )IN1)->numCols; \ + numRows = ((psImage* )IN1)->numRows; \ + numCols = ((psImage* )IN1)->numCols; \ for(j = 0; j < numCols; j++) { \ - o = ((psImage* )OUT)->data.TYPE[j]; \ - i1 = &((psScalar* )IN1)->data.TYPE; \ - i2 = ((psImage* )IN2)->data.TYPE[j]; \ + o = ((psImage* )OUT)->data.TYPE[j]; \ + i1 = &((psScalar* )IN1)->data.TYPE; \ + i2 = ((psImage* )IN2)->data.TYPE[j]; \ for(i = 0; i < numRows; i++, o++, i2++) { \ *o = OP; \ @@ -165,8 +165,8 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - n1 = ((psVector* )IN1)->n; \ - o = ((psVector* )OUT)->data.TYPE; \ - i1 = ((psVector* )IN1)->data.TYPE; \ - i2 = &((psScalar* )IN2)->data.TYPE; \ + n1 = ((psVector* )IN1)->n; \ + o = ((psVector* )OUT)->data.TYPE; \ + i1 = ((psVector* )IN1)->data.TYPE; \ + i2 = &((psScalar* )IN2)->data.TYPE; \ for (i=0; i < n1; i++, o++, i1++) { \ *o = OP; \ @@ -182,13 +182,13 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - n1 = ((psVector* )IN1)->n; \ - n2 = ((psVector* )IN2)->n; \ + n1 = ((psVector* )IN1)->n; \ + n2 = ((psVector* )IN2)->n; \ if(n1 != n2) { \ psError(__func__, ": Inconsistent element count: %d vs %d", n1, n2); \ return OUT; \ } \ - o = ((psVector* )OUT)->data.TYPE; \ - i1 = ((psVector* )IN1)->data.TYPE; \ - i2 = ((psVector* )IN2)->data.TYPE; \ + o = ((psVector* )OUT)->data.TYPE; \ + i1 = ((psVector* )IN1)->data.TYPE; \ + i2 = ((psVector* )IN2)->data.TYPE; \ for (i=0; i < n1; i++, o++, i1++, i2++) { \ *o = OP; \ @@ -207,8 +207,8 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - n1 = ((psVector* )IN1)->n; \ - dim1 = ((psVector* )IN1)->type.dimen; \ - numRows2 = ((psImage* )IN2)->numRows; \ - numCols2 = ((psImage* )IN2)->numCols; \ + n1 = ((psVector* )IN1)->n; \ + dim1 = ((psVector* )IN1)->type.dimen; \ + numRows2 = ((psImage* )IN2)->numRows; \ + numCols2 = ((psImage* )IN2)->numCols; \ \ if(dim1 == PS_DIMEN_VECTOR) { /* Regular vectors */ \ @@ -218,8 +218,8 @@ } \ \ - i1 = ((psVector* )IN1)->data.TYPE; \ + i1 = ((psVector* )IN1)->data.TYPE; \ for(j = 0; j < numRows2; j++, i1++) { \ - o = ((psImage* )OUT)->data.TYPE[j]; \ - i2 = ((psImage* )IN2)->data.TYPE[j]; \ + o = ((psImage* )OUT)->data.TYPE[j]; \ + i2 = ((psImage* )IN2)->data.TYPE[j]; \ for(i = 0; i < numCols2; i++, o++, i2++) { \ *o = OP; \ @@ -233,7 +233,7 @@ \ for(j = 0; j < numRows2; j++) { \ - o = ((psImage* )OUT)->data.TYPE[j]; \ - i1 = ((psVector* )IN1)->data.TYPE; \ - i2 = ((psImage* )IN2)->data.TYPE[j]; \ + o = ((psImage* )OUT)->data.TYPE[j]; \ + i1 = ((psVector* )IN1)->data.TYPE; \ + i2 = ((psImage* )IN2)->data.TYPE[j]; \ for(i = 0; i < numCols2; i++, o++, i1++, i2++) { \ *o = OP; \ @@ -253,10 +253,10 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - numRows = ((psImage* )IN1)->numRows; \ - numCols = ((psImage* )IN1)->numCols; \ + numRows = ((psImage* )IN1)->numRows; \ + numCols = ((psImage* )IN1)->numCols; \ for(j = 0; j < numRows; j++) { \ - o = ((psImage* )OUT)->data.TYPE[j]; \ - i1 = ((psImage* )IN1)->data.TYPE[j]; \ - i2 = &((psScalar* )IN2)->data.TYPE; \ + o = ((psImage* )OUT)->data.TYPE[j]; \ + i1 = ((psImage* )IN1)->data.TYPE[j]; \ + i2 = &((psScalar* )IN2)->data.TYPE; \ for(i = 0; i < numCols; i++, o++, i1++) { \ *o = OP; \ @@ -276,8 +276,8 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - n2 = ((psVector* )IN2)->n; \ - dim2 = ((psVector* )IN2)->type.dimen; \ - numRows1 = ((psImage* )IN1)->numRows; \ - numCols1 = ((psImage* )IN1)->numCols; \ + n2 = ((psVector* )IN2)->n; \ + dim2 = ((psVector* )IN2)->type.dimen; \ + numRows1 = ((psImage* )IN1)->numRows; \ + numCols1 = ((psImage* )IN1)->numCols; \ \ if(dim2 == PS_DIMEN_VECTOR) { /* Regular vectors */ \ @@ -287,8 +287,8 @@ } \ \ - i2 = ((psVector* )IN2)->data.TYPE; \ + i2 = ((psVector* )IN2)->data.TYPE; \ for(j = 0; j < numRows1; j++, i1++) { \ - o = ((psImage* )OUT)->data.TYPE[j]; \ - i1 = ((psImage* )IN1)->data.TYPE[j]; \ + o = ((psImage* )OUT)->data.TYPE[j]; \ + i1 = ((psImage* )IN1)->data.TYPE[j]; \ for(i = 0; i < numCols1; i++, o++, i1++) { \ *o = OP; \ @@ -302,7 +302,7 @@ \ for(j = 0; j < numRows1; j++) { \ - o = ((psImage* )OUT)->data.TYPE[j]; \ - i1 = ((psVector* )IN2)->data.TYPE; \ - i2 = ((psImage* )IN1)->data.TYPE[j]; \ + o = ((psImage* )OUT)->data.TYPE[j]; \ + i1 = ((psVector* )IN2)->data.TYPE; \ + i2 = ((psImage* )IN1)->data.TYPE[j]; \ for(i = 0; i < numCols1; i++, o++, i2++, i1++) { \ *o = OP; \ @@ -323,8 +323,8 @@ ps##TYPE *i1 = NULL; \ ps##TYPE *i2 = NULL; \ - numRows1 = ((psImage* )IN1)->numRows; \ - numCols1 = ((psImage* )IN1)->numCols; \ - numRows2 = ((psImage* )IN2)->numRows; \ - numCols2 = ((psImage* )IN2)->numCols; \ + numRows1 = ((psImage* )IN1)->numRows; \ + numCols1 = ((psImage* )IN1)->numCols; \ + numRows2 = ((psImage* )IN2)->numRows; \ + numCols2 = ((psImage* )IN2)->numCols; \ if(numRows1!=numRows2 || numCols1!=numCols2) { \ psError(__func__, ": Inconsistent element count: numRows: %d vs %d numCols: %d vs %d", numRows1, \ @@ -333,7 +333,7 @@ } \ for(j = 0; j < numRows1; j++) { \ - o = ((psImage* )OUT)->data.TYPE[j]; \ - i1 = ((psImage* )IN1)->data.TYPE[j]; \ - i2 = ((psImage* )IN2)->data.TYPE[j]; \ + o = ((psImage* )OUT)->data.TYPE[j]; \ + i1 = ((psImage* )IN1)->data.TYPE[j]; \ + i2 = ((psImage* )IN2)->data.TYPE[j]; \ for(i = 0; i < numCols1; i++, o++, i1++, i2++) { \ *o = OP; \