Q_rhs_vec is not set up properly in IBFEMethod::computeLagrangianFluidSource

[sharpgao]

IBAMR/src/IB/IBFEMethod.cpp

Line 888 in ac4614e

Q_dof_map.constrain_element_vector(Q_rhs_e, Q_dof_indices);

In IBFEMethod::computeLagrangianFluidSource(double data_time), the following two lines 888-889 need to move out of the for loop on qp,
// Apply constraints (e.g., enforce periodic boundary conditions) // and add the elemental contributions to the global vector. Q_dof_map.constrain_element_vector(Q_rhs_e, Q_dof_indices); Q_rhs_vec->add_vector(Q_rhs_e, Q_dof_indices);