> BTW, there is a bug in dof_map.C, which I check out from svn, > around line 915, > >>> if (var == n_vars-1) > _var_first_local_df.push_back(next_free_dof); > > _var_first_local_df[nvar+1] record the offset of vars, like > [u_begin, v_begin, w_begin, p_begin, p_end], > but without "if (var == n_vars-1)", it becomes > [u_begin, u_end, v_begin, v_end, w_begin, w_end, p_begin, p_end]
This actually the expected behavior. Consider v_begin to be u_end in the C++ container sense - so to get all the u's: u: [u_begin, v_begin) v: [v_begin, w_begin) -Ben ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
