* ( &                                 (ele_val_at_quad(vel_prod,ele)) &                                  ) & !detwei                                  ) !shape_rhs             rhs_addto_disip = shape_rhs(shape_psi, detwei