Re: [Libmesh-devel] [Libmesh-users] block matrix implementation

2013-05-31 Thread Kirk, Benjamin (JSC-EG311)
On May 31, 2013, at 4:04 AM, Lorenzo Alessio Botti wrote: > I don't know if it is the "correct" way to do it but the memory footprint has > not changed. > I'm also interested in knowing if the memory savings take place on the > libMesh side or on the PETSc side, or on both sides. That's wha

Re: [Libmesh-devel] [Libmesh-users] block matrix implementation

2013-05-31 Thread Lorenzo Alessio Botti
I tried to modify dof_map.h as follows unsigned int block_size() const { #ifdef LIBMESH_ENABLE_BLOCKED_STORAGE return (this->has_blocked_representation() ? this->n_variables()*4 : 1); #else return 1; #endif } I don't know if it is the "correct" way to do it but the memory footprint