thanks!! Full error message is as follow
*/Users/jaekwangjk/repos/trial/sphere-deal.ii/mystokes.cc:1258:7: **error: *
*no*
* matching function for call to 'integrate_difference'*
VectorTools::integrate_difference (dof_handler,
* ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
*/Users/jaekwangjk/repos/trial/sphere-deal.ii/mystokes.cc:1220:7: note: *in
instantiation of member function
'MyStokes::StokesProblem<2>::process_solution' requested here
process_solution ();
* ^*
*/Users/jaekwangjk/repos/trial/sphere-deal.ii/mystokes.cc:1298:21: note: *in
instantiation of member function 'MyStokes::StokesProblem<2>::run'
requested here
flow_problem.run ();
* ^*
*/Users/jaekwangjk/Programs/dealii/include/deal.II/numerics/vector_tools.h:1940:8:
note: *
candidate template ignored: could not match 'Function' against
'Solution'
void integrate_difference (const DoFHandler<dim,spacedim> &dof,
* ^*
*/Users/jaekwangjk/Programs/dealii/include/deal.II/numerics/vector_tools.h:1968:8:
note: *
candidate template ignored: could not match 'dealii::hp::DoFHandler'
against 'dealii::DoFHandler'
void integrate_difference (const hp::DoFHandler<dim,spacedim> &dof,
* ^*
*/Users/jaekwangjk/Programs/dealii/include/deal.II/numerics/vector_tools.h:1925:8:
note: *
candidate function template not viable: requires at least 7
arguments, but
6 were provided
void integrate_difference (const Mapping<dim,spacedim> &mapping,
* ^*
*/Users/jaekwangjk/Programs/dealii/include/deal.II/numerics/vector_tools.h:1953:8:
note: *
candidate function template not viable: requires at least 7
arguments, but
6 were provided
void integrate_difference (const hp::MappingCollection<dim,spacedim>
&mapping,
I will check the mapping as you advised for my other questions.!!!!!!
2016년 9월 20일 화요일 오전 9시 51분 13초 UTC-5, Martin Kronbichler 님의 말:
>
> Jaekwang,
>
> can you post the full error message? To me the call looks reasonable (this
> is a copy from step-7) and the function definition looks good as well.
>
> Ah, and in case you want to use this in the context of your other question
> regarding high order methods: Put a "mapping" as first argument to this
> call and make sure to use QGauss<dim>(fe.degree+2) or so, which increases
> the order of quadrature as the element order gets higher.
>
> Best,
> Martin
>
> On 09/20/2016 04:46 PM, JAEKWANG KIM wrote:
>
>
>
> I wanted to evaluate Vectortool::integrate_difference to estimate my
> error.
>
>
>
> but I couldn't figure out what I need to enter for the third component. (I
> marked it as red)
>
>
>
>
> VectorTools::integrate_difference (dof_handler,
>
> solution,
>
> *Solution<dim+**1**>()*,
>
> difference_per_cell,
>
> QGauss<dim>(3),
>
> VectorTools::L2_norm);
>
>
> I have solution vector over my domain,(by the way I am solving vector
> valued problems...)
>
> I want to compare this with continuous exact function, "May be
> Solution<dim+1>" , I have generated it as follow.
>
> However, I got a error message that
>
>
> */Users/jaekwangjk/repos/trial/sphere-deal.ii/mystokes.cc:1258:7: **error:
> **no*
>
> * matching function for call to 'integrate_difference'*
>
>
> What's might be wrong with my code......
>
>
> template <int dim>
>
> class Solution : public Function<dim>
>
> {
>
> public:
>
>
>
> Solution () : Function<dim>(dim+1) {}
>
>
>
> //const unsigned int b_id;
>
> virtual double value (const Point<dim> &p,
>
> const unsigned int component = 0) const;
>
>
>
> virtual void vector_value (const Point<dim> &p,
>
> Vector<double> &value) const;
>
>
>
> };
>
>
>
>
>
> template <int dim>
>
> double
>
> Solution<dim>::value (const Point<dim> &p,
>
> const unsigned int component) const
>
> {
>
> Assert (component < this->n_components,
>
> ExcIndexRange (component, 0, this->n_components));
>
>
>
> double U=-1.; //inflow velocity
>
> double obj_rad=0.2; //Object Radius
>
> double r_sph = sqrt( pow(p[0],2) + pow(p[1],2) );
>
> double phi = (-1) * acos(p[1]/r_sph);
>
>
>
> double q_r_sph=U*cos(phi)*( 1+
> 0.5*pow(obj_rad,3)/pow(r_sph,3)-1.5*obj_rad/r_sph
> );
>
> double q_phi=(-1)* U*sin(phi)*(1 -0.25 *pow(obj_rad,3)/pow(r_sph,3
> )-0.75*obj_rad/r_sph);
>
>
>
> double q_x=q_r_sph*sin(phi)+q_phi*cos(phi);
>
> double q_y=(-1)*q_r_sph*cos(phi)+q_phi*sin(phi);
>
>
>
>
>
> if (component == 0)
>
> return q_x;
>
> if (component == 1)
>
> return q_y;
>
> else
>
> return 0; //
>
>
>
>
>
> }
>
>
> template <int dim>
>
> void
>
> Solution<dim>::vector_value (const Point<dim> &p,
>
> Vector<double> &values) const
>
> {
>
> for (unsigned int c=0; c<this->n_components; ++c)
>
> values(c) = Solution<dim>::value (p, c);
>
> }
>
> --
> 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+un...@googlegroups.com <javascript:>.
> For more options, visit https://groups.google.com/d/optout.
>
>
>
--
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.
For more options, visit https://groups.google.com/d/optout.
