> There are no constraints on the equation. It is a first order time dependent > problem, so there are just initial conditions that are taken care of > separately.
and boundary conditions I would assume? Otherwise your system is not invertible. > The code is for 3D brick elements with 9 dofs per node. I just use > dof_hanlder.distribute_dofs() function to number the nodes. Can this result > in weird couplings ? It depends on how many entries you have per row in your sparsity pattern (even if the values are zero). You should check that to be sure. > Also, does deal.ii support superLU_dist ? I mean I know Petsc does but > searching for superLU + deal.ii > doesn't give any useful results to me. No, but you should be able to easily create a class for this or it might even be possible to just to use PreconditionLU and set "-pc_factor_mat_solver_package superlu" on the command line. -- Timo Heister http://www.math.clemson.edu/~heister/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
