Method addMatrixTranspose

SRC/matrix/Matrix.cpp:720–793

int
Matrix::addMatrixTranspose(double factThis, const Matrix &other, double factOther)
{
    if (factThis == 1.0 && factOther == 0.0)
      return 0;

#ifdef _G3DEBUG
    if ((other.numRows != numCols) || (other.numCols != numRows)) {
      opserr << "Matrix::addMatrixTranspose(): incompatible matrices\n";
      return -1;
    }
#endif

    if (factThis == 1.0) {

      // want: this += other^T * factOther
      if (factOther == 1.0) {
    double *dataPtr = data;
    for (int j=0; j<numCols; j++) {
      for (int i=0; i<numRows; i++)
	    *dataPtr++ += (other.data)[j+i*numCols];
    }
      } else {
	double *dataPtr = data;
    for (int j=0; j<numCols; j++) {
      for (int i=0; i<numRows; i++)
	    *dataPtr++ += (other.data)[j+i*numCols] * factOther;
    }
      }
    } 

    else if (factThis == 0.0) {

      // want: this = other^T * factOther
      if (factOther == 1.0) {
	double *dataPtr = data;
    for (int j=0; j<numCols; j++) {
      for (int i=0; i<numRows; i++)
	    *dataPtr++ = (other.data)[j+i*numCols];
    }
      } else {
	double *dataPtr = data;
    for (int j=0; j<numCols; j++) {
      for (int i=0; i<numRows; i++)
	    *dataPtr++ = (other.data)[j+i*numCols] * factOther;
    }
      }
    } 

    else {

      // want: this = this * thisFact + other^T * factOther
      if (factOther == 1.0) {
	double *dataPtr = data;
    for (int j=0; j<numCols; j++) {
      for (int i=0; i<numRows; i++) {
        double value = *dataPtr * factThis + (other.data)[j+i*numCols];
	    *dataPtr++ = value;