I accidentally replied off-list, oops ...

--- Begin Message ---
Magnus Hagdorn writes:
 > Hello world,
 
 > Do you have estimates for the memory requirements for the stiffness
 > matrix and load vector? What size problems can be reasonably run?

The stiffness matrix will always require way more storage, so let's
concentrate on that.  It depends on the number of variables Nv in your
system and the connectivity of the mesh.  If we consider e.g. a
structured 3D mesh of 8-node hexahedra, each node has 27 neighbors
(including itself).  Each sparse matrix row will in general contain
27*Nv nonzero entries.  Supposing that your structured grid contains
Nx*Ny*Nz nodes, the total number of nonzeros will be Nx*Ny*Nz*27*Nv.
Assuming each nonzero requires 8 bytes of storage (double precision)
and maybe another 8 bytes or so of sparsity information, the total
storage requirement is approximately Nx*Ny*Nz*27*Nv*16 bytes.

Concrete example: 20^3 nodes, 3 variables per node =
20*20*20*27*3*16 = 10368000 bytes = about 10 MB

If you are using tetrahedra, I believe the number of nonzeros per row
jumps dramatically due to the connectivity, so you have to take this
into account.

The DofMap and Mesh have recently been parallelized by other developers,
so I think I will let them comment a bit more on that aspect of storage.

-J

--- End Message ---
-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2008.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to