From: David K. <dknez@MIT.EDU> - 2009-04-30 17:50:41
|
Hi Roy, This patch should take care of those changes. - Dave Roy Stogner wrote: > > On Tue, 28 Apr 2009, David Knezevic wrote: > >> So, I've added the multiple SCALARs stuff (see the patch), and that >> works fine. To do this I added a std::map member to DofMap to store >> the Order of each SCALAR. > > Do we need that? We've already committed the sin of duplicating data > once, so we might as well take advantage of that: we should be able to > grab the variable order data from DofMap::_variables. > >> Also, I added the SCALAR code from DofMap::dof_indices to >> DofMap::old_dof_indices, with n_old_dofs() in place of n_dofs(). > > Thanks; that will be important. > >> However, I realized that a system with a SCALAR variable breaks on >> calling equation_systems.reinit() after doing (say) a uniform >> refinement, since vector projection calls fe_scalar_shape_2D.C. I >> haven't looked at this in detail yet, but I guess we should just map >> each SCALAR value in the old vector to the corresponding entry in the >> projected vector... > > Right. And if that's too much work for now, just stick a > libmesh_not_implemented() inside a loop testing for SCALARs in > System::project_vector(), I think. Likewise, let's switch those > libmesh_error() calls in DofMap to libmesh_not_implemented() > --- > Roy |