| 977 | } |
| 978 | |
| 979 | void |
| 980 | FE_Element::addM_ForceSensitivity(int gradNumber, const Vector &vect, double fact) |
| 981 | { |
| 982 | // Get the components we need out of the vector |
| 983 | // and place in a temporary vector |
| 984 | Vector tmp(numDOF); |
| 985 | for (int i=0; i<numDOF; i++) { |
| 986 | int loc = myID(i); |
| 987 | if (loc >= 0) { |
| 988 | tmp(i) = vect(loc); |
| 989 | } |
| 990 | else { |
| 991 | tmp(i) = 0.0; |
| 992 | } |
| 993 | } |
| 994 | if (theResidual->addMatrixVector(1.0, myEle->getMassSensitivity(gradNumber),tmp,fact) < 0) { |
| 995 | opserr << "WARNING FE_Element::addM_ForceSensitivity() - "; |
| 996 | opserr << "- addMatrixVector returned error\n"; |
| 997 | } |
| 998 | } |
| 999 | |
| 1000 | void |
| 1001 | FE_Element::addD_ForceSensitivity(int gradNumber, const Vector &vect, double fact) |
no test coverage detected