| 193 | } |
| 194 | |
| 195 | double |
| 196 | Pressure_Constraint::getPdot(int last) |
| 197 | { |
| 198 | if (pval != 0) { |
| 199 | return pval[1]; |
| 200 | } |
| 201 | |
| 202 | Domain* theDomain = this->getDomain(); |
| 203 | if(theDomain == 0) { |
| 204 | opserr<<"WARNING: domain has not been set"; |
| 205 | opserr<<" -- Pressure_Constraint::getPressureNode\n"; |
| 206 | return 0; |
| 207 | } |
| 208 | |
| 209 | Node* pNode = theDomain->getNode(pTag); |
| 210 | if(pNode == 0) return 0.0; |
| 211 | const Vector& accel = pNode->getAccel(); |
| 212 | if(last == 1) { |
| 213 | if(accel.Size()==0) return 0.0; |
| 214 | return accel(0); |
| 215 | } |
| 216 | return 0.0; |
| 217 | } |
| 218 | |
| 219 | const ID& |
| 220 | Pressure_Constraint::getFluidElements() |