| 134 | } |
| 135 | |
| 136 | void |
| 137 | Pressure_Constraint::setPressure(double p) |
| 138 | { |
| 139 | if (pval != 0) { |
| 140 | pval[0] = p; |
| 141 | return; |
| 142 | } |
| 143 | |
| 144 | Node* pnode = this->getPressureNode(); |
| 145 | if (pnode == 0) return; |
| 146 | const Vector& vel = pnode->getVel(); |
| 147 | Vector newvel(vel); |
| 148 | newvel.Zero(); |
| 149 | newvel(0) = p; |
| 150 | pnode->setTrialVel(newvel); |
| 151 | pnode->commitState(); |
| 152 | } |
| 153 | |
| 154 | void |
| 155 | Pressure_Constraint::setPdot(double pdot) |
nothing calls this directly
no test coverage detected