Dear All,
I have gone through a lot of tutorials. I would like to understand the basic relation between dof_handler and FESystem. I am trying to solve *two* 4th order equations on the same mesh (or triangulation). And each equation is split into two second order equations. FESystem<dim> fe_one, fe_two; BlockSparseMatrix<double> system_matrix_one, system_matrix_two; BlockVector<double> solution_one, solution_two; BlockVector<double> old_solution_one, old_solution_two; BlockVector<double> system_rhs_one, system_rhs_two; Question (1) Do I need one or two dof_handler? DoFHandler<dim> dof_handler; Where, fe_one(FE_Q<dim>(degree),1, FE_Q<dim>(degree),1) fe_two(FE_Q<dim>(degree),1, FE_Q<dim>(degree),1) Some additional information, both the fe_one and fe_two systems are *independently solved* in a semi-implicitly manner. Meaning system_matrix_one is a function of old_solution_one and old_solution_two, likewise for system_matrix_two. Question (2) I am assuming one block spartisty_pattern is enough and use the same for both the system_matrix? BlockSparsityPattern sparsity_pattern; If there would be other multiple declarations please let me know. I look forward to your response. Best regards, Karthi. -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/465e55a9-9d17-4bcf-ac51-209217e5467an%40googlegroups.com.
