Took a brief hiatus. and managed to pass a huddle following Marc's
suggestion.
Now I have set to construct convergence tables for the polymer flow problem
I am solving. It involves solving a Stokes Problem for velocity and
pressure (u,p), (size dim+1) and the Constitutive Equations for the stress
tensor tau, (size dim*dim). The solution has 7 components for dim=2.
To calculate the errors, I need the exact solution [(u,p), tau], ( size dim
+1 + dim*dim). To pick the correct component I use the
ComponentSelectFunction:
const ComponentSelectFunction<dim*dim+dim+1>
velocity_mask(std::make_pair(0, dim), dim*dim+dim+1);
const ComponentSelectFunction<dim*dim+dim+1>
pressure_mask (dim, dim+3);
const ComponentSelectFunction<dim*dim+dim+1>
stress_mask(std::pair(dim+3, dim*dim+dim+1), dim*dim+dim+1);
Now the errors that I am getting are:
/home/jurombo/binaires/dealii/nnf/nnewton/oldroyd_b_stab/oldroyd_bn.cc:951:43:
error:
no matching fu
nction for call to ‘integrate_difference(dealii::MappingQ<2, 2>&,
dealii::DoFHandler<2, 2>&, dealii:
:BlockVector<double>&, Viscoelastic::ExactSolution<7>&,
dealii::Vector<double>&, dealii::QGauss<7>&,
dealii::VectorTools::NormType, const dealii::ComponentSelectFunction<7,
double>*)’
951 | VectorTools::integrate_difference (mapping, dof_handler,
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
952 | solution,
| ~~~~~~~~~
953 | exact_solution,
| ~~~~~~~~~~~~~~~
954 | cellwise_errors,
| ~~~~~~~~~~~~~~~~
955 | quadrature,
| ~~~~~~~~~~~
956 | VectorTools::H1_seminorm,
| ~~~~~~~~~~~~~~~~~~~~~~~~~
957 | &velocity_mask);
| ~~~~~~~~~~~~~~~
/usr/local/share/dealii/include/deal.II/numerics/vector_tools_integrate_difference.h:143:3:
note:
ca
ndidate: ‘template<int dim, class InVector, class OutVector, int spacedim>
void dealii::VectorTools:
:integrate_difference(const dealii::Mapping<dim, spacedim>&, const
dealii::DoFHandler<dim, spacedim>
&, const InVector&, const dealii::Function<spacedim, typename
InVector::value_type>&, OutVector&, co
nst dealii::Quadrature<dim>&, const NormType&, const
dealii::Function<spacedim, double>*, double)’
143 | integrate_difference(
| ^~~~~~~~~~~~~~~~~~~~
/usr/local/share/dealii/include/deal.II/numerics/vector_tools_integrate_difference.h:143:3:
note:
template argument deduction/substitution failed:
/home/jurombo/binaires/dealii/nnf/nnewton/oldroyd_b_stab/oldroyd_bn.cc:951:43:
note:
deduced confl
icting values for non-type parameter ‘spacedim’ (‘2’ and ‘7’)
951 | VectorTools::integrate_difference (mapping, dof_handler,
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
952 | solution,
| ~~~~~~~~~
953 | exact_solution,
| ~~~~~~~~~~~~~~~
954 | cellwise_errors,
| ~~~~~~~~~~~~~~~~
955 | quadrature,
| ~~~~~~~~~~~
956 | VectorTools::H1_seminorm,
| ~~~~~~~~~~~~~~~~~~~~~~~~~
957 | &velocity_mask);
| ~~~~~~~~~~~~~~~
/usr/local/share/dealii/include/deal.II/numerics/vector_tools_integrate_difference.h:160:3:
note:
ca
ndidate: ‘template<int dim, class InVector, class OutVector, int spacedim>
void dealii::VectorTools:
:integrate_difference(const dealii::DoFHandler<dim, spacedim>&, const
InVector&, const dealii::Funct
ion<spacedim, typename InVector::value_type>&, OutVector&, const
dealii::Quadrature<dim>&, const Nor
mType&, const dealii::Function<spacedim, double>*, double)’
The issues seem to be emanating from the form of the *exact_solution *and
how the components are selected.
Ant suggestions on resolving this?
----------------------------------------------------------------------------
On Monday, 6 February 2023 at 03:11:19 UTC+2 wrote:
Hello Jack,
the error message is quite descriptive and tells you that something is off
with the number of components of your Function object:
*The violated condition was: exact_solution.n_components ==
n_components*
Check line 200 of your code:
ExactSolution () : Function<dim>(dim+1) {}
You call the constructor for Function with 3 components (since dim=2 in
your case). However, your implementation returns 7 components. Make sure
that you call the constructor with the right number of components.
Marc
On Saturday, February 4, 2023 at 8:54:47 AM UTC-7 jack urombo wrote:
i am using the implimenation by Jaekwang Kim to model polymer flows. The
major change that I want to do is to output convergence tables like in
Step-7 and Step-22, which means I have to have an ExactSolution.
The code is encountering an error because the dimension of the solutions do
not match at a some point when generating output:
*An error occurred in line <820> of file
</home/jurombo/Downloads/dealii-9.3.2/include/deal.II/numerics/vector_tools_integrate_difference.templates.h>
in function void
dealii::VectorTools::internal::do_integrate_difference(const
dealii::hp::MappingCollection<dim, spacedim>&, const
dealii::DoFHandler<dim, spacedim>&, const InVector&, const
dealii::Function<spacedim, typename InVector::value_type>&, OutVector&,
const dealii::hp::QCollection<dim>&, const dealii::VectorTools::NormType&,
const dealii::Function<spacedim>*, double) [with int dim = 2; int spacedim
= 2; InVector = dealii::BlockVector<double>; OutVector =
dealii::Vector<double>; typename InVector::value_type = double] The
violated condition was: exact_solution.n_components == n_components
Additional information: Dimension 3 not equal to 7. Stacktrace:
----------- #0 /usr/local/share/dealii/lib/libdeal_II.g.so.9.3.2: #1
/usr/local/share/dealii/lib/libdeal_II.g.so.9.3.2: void
dealii::VectorTools::integrate_difference<2, dealii::BlockVector<double>,
dealii::Vector<double>, 2>(dealii::Mapping<2, 2> const&,
dealii::DoFHandler<2, 2> const&, dealii::BlockVector<double> const&,
dealii::Function<2, dealii::BlockVector<double>::value_type> const&,
dealii::Vector<double>&, dealii::Quadrature<2> const&,
dealii::VectorTools::NormType const&, dealii::Function<2, double> const*,
double) #2 ./oldroyd_bn:
Viscoelastic::StokesProblem<2>::output_results(unsigned int) #3
./oldroyd_bn: Viscoelastic::StokesProblem<2>::run() #4 ./oldroyd_bn: main
--------------------------------------------------------*
I have checked the ExactSolution in the code to try and figure wher the
error is occurring without success.
Can anyone point me the direction to look.
--
The information in this message is confidential and legally
privileged.
It is intended solely for the addressee(s). Access to this message
by
anyone else is unauthorized. If received in error, please accept our
apologies
and notify the sender immediately. You must also delete the
original message
from your machine. If you are not the intended recipient,
any use, disclosure,
copying, distribution or action taken in reliance of
it, is prohibited and may
be unlawful. The information, attachments,
opinions or advice contained in this
email are not the views or opinions of
Harare Institute of Technology, its
subsidiaries or affiliates. Although
this email and any attachments are
believed to be free of any virus or
other defects which might affect any
computer or IT system into which they
are received, no responsibility is
accepted by Harare Institute of
Technology and/or its subsidiaries for any loss
or damage arising in any
way from the receipt or use thereof.
--
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/c77b68f3-c243-49f5-977e-15d5cd183937n%40googlegroups.com.
