| 981 | } // setConstraintForce |
| 982 | |
| 983 | void |
| 984 | CIBMethod::getConstraintForce(Vec* L, const double data_time) |
| 985 | { |
| 986 | #if !defined(NDEBUG) |
| 987 | TBOX_ASSERT(IBTK::rel_equal_eps(data_time, d_current_time) || IBTK::rel_equal_eps(data_time, d_new_time)); |
| 988 | #else |
| 989 | NULL_USE(data_time); |
| 990 | #endif |
| 991 | const int struct_ln = getStructuresLevelNumber(); |
| 992 | Pointer<LData> ptr_lagmultpr = d_l_data_manager->getLData("lambda", struct_ln); |
| 993 | Vec lambda = ptr_lagmultpr->getVec(); |
| 994 | *L = lambda; |
| 995 | |
| 996 | return; |
| 997 | } // getConstraintForce |
| 998 | |
| 999 | void |
| 1000 | CIBMethod::getFreeRigidVelocities(Vec* U, const double data_time) |
no test coverage detected