Hi all
I am looking at step-43 using different method. I do a semi discretization
in space and end up with an ode for the saturation and use different
technique for time integration for the saturation.
For this I will need to assemble the matrix H after solving for the
velocity in the darcy system. Individual entries of H are given by:
\mathbf{H}_{ij} & = - \Delta t^{(n)}_c \Big( F\left(S^{(n-1)}\right)
\mathbf{v}_i,\nabla\phi_j\Big)_{\Omega}
I am struggling with sparsity pattern that I have to define for this
matrix. I try this:
assemble_matrix_H.clear();
DynamicSparsityPattern dsp (n_s, n_u);
DoFTools::make_sparsity_pattern (saturation_dof_handler, darcy_dof_handler,
dsp,
saturation_constraints, false);
ssemble_matrix_H.reinit(dsp);
dealii is complaining and is saying
/home/franckm/apps/candi/deal.II-toolchain/deal.II-v8.4.0/include/deal.II/dofs/dof_tools.h:563:3:
note: template<class DoFHandlerType, class SparsityPatternType> void
dealii::DoFTools::make_sparsity_pattern(const DoFHandlerType&, const
DoFHandlerType&, SparsityPatternType&)
make_sparsity_pattern (const DoFHandlerType &dof_row,
^
/home/franckm/apps/candi/deal.II-toolchain/deal.II-v8.4.0/include/deal.II/dofs/dof_tools.h:563:3:
note: template argument deduction/substitution failed:
/home/franckm/candi-examples/step-43b/step-43b.cc:918:68: note: candidate
expects 3 arguments, 5 provided
saturation_constraints, false);
^
make[3]: *** [CMakeFiles/step-43b.dir/step-43b.cc.o] Error 1
How one define a sparsity pattern for assembling H in the above scenario?
cheers
--
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.
