Dear all,
I am trying to extend my implementation with parallel computing. I would
like to get the solution in a single point when having a distributed
system. For example, in step-40 I have the *locally_relevant_solution* at
the output_result(...) function.
In a non distributed implementation I u
Dear Hermes,
Did you try VectorTools::point_values,
https://www.dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#a1ad6eceb6cbeaa505baf7f938289bbde
?
As an alternative, you can also check if the point is in the locally
owned cell of the respective cell, and run the evaluation only
Thank you for pointing that function.
I am having some problems when creating the mapping and I am not sure if
the function call is right. Could you please help me with that?
std::vector<*double*> pointSol;
Mapping mapping;
dealii::Utilities::MPI::RemotePointEvaluation remot
Hi Hermes,
You don't need to do anything regarding the setup (it is done within the
function). Just take at a look
at:
https://github.com/dealii/dealii/blob/44b6aadb35aca2333e4bfb6c9ce29c613f4dc9e9/tests/remote_point_evaluation/vector_tools_evaluate_at_points_01.cc#L214-L216
I'll extend the do
Thank you.
Doing something similar:
*const* MappingQ1 mapping;
Utilities::MPI::RemotePointEvaluation evaluation_cache;
VectorTools::point_values(mapping, dof_handler,
locally_relevant_solution, Point<2>(0.2, 0.2), evaluation_cache);
I get the following error:
error: no mat
You need to pass a vector of points.
PM
On Tuesday, 24 August 2021 at 15:26:05 UTC+2 hermes...@gmail.com wrote:
> Thank you.
>
> Doing something similar:
>
> *const* MappingQ1 mapping;
>
> Utilities::MPI::RemotePointEvaluation evaluation_cache;
>
> VectorTools::point_values(map
Hello everyone!
This is deal.II newsletter #179.
It automatically reports recently merged features and discussions about the
deal.II finite element library.
## Below you find a list of recently proposed or merged features:
#12703: Require GINKGO 1.4.0 (proposed by gfcas)
https://github.com/de
Ok, I see - it is that EPS_TARGET_MAGNITUDE flag that was missing. It works
now! Thanks a lot, Simon.
Anton.
On Tuesday, August 24, 2021 at 12:14:20 AM UTC-6 simon...@gmail.com wrote:
> Hi,
>
> I think you need to set both which eigen pair and the target value. For
> example, if you look for th
