Thank you for your answer, I could make it work.
Regards,
Hermes
El mar, 7 sept 2021 a las 4:28, Wolfgang Bangerth ()
escribió:
> On 9/6/21 3:27 AM, Hermes Sampedro wrote:
> > Thank you very much. I have been trying and I only get the real part of
> the
> > solution and it seems not to be right
On 9/6/21 3:27 AM, Hermes Sampedro wrote:
Thank you very much. I have been trying and I only get the real part of the
solution and it seems not to be right so far. Would be possible to make it
work for complex numbers?
Hermes,
do you think you could come up with a small, self-contained
Thank you very much. I have been trying and I only get the real part of the
solution and it seems not to be right so far. Would be possible to make it
work for complex numbers?
Thank you
H.
El vie, 27 ago 2021 a las 9:33, 'peterrum' via deal.II User Group (<
dealii@googlegroups.com>) escribió:
Hi Hermes,
I have extended the documentation a bit.
See https://github.com/dealii/dealii/pull/12714.
Maybe you could give us a short feedback if `VectorTools::point_values()`
is working for complex numbers. If this is not the case, we can extend the
function for this use case.
PM
On
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;
>
>
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
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
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
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
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
Thank you very much, I could make it work.
Regards,
Hermes
El mar, 10 ago 2021 a las 19:46, Jean-Paul Pelteret ()
escribió:
> Another thing: You probably also need to initialise with the right number
> of components. So something like
> Vector vecSol (fe.n_components());
>
>
> On 10. Aug 2021,
Another thing: You probably also need to initialise with the right number of
components. So something like
Vector vecSol (fe.n_components());
> On 10. Aug 2021, at 19:40, Jean-Paul Pelteret wrote:
>
> Hi Hermes,
>
> You don’t say what errors you’re seeing, but my guess is that it now doesn’t
Hi Hermes,
You don’t say what errors you’re seeing, but my guess is that it now doesn’t
compile. This variant of the function (the one that Daniel linked to) returns
void, so you should call it before outputting the result:
Vector vecSol;
VectorTools::point_value(dof_handler,
Thank you for your answer. I am not sure if I fully understand your
suggestion. Do you mean something like that:
Vector<*double*> vecSol;
std::cout << "Solution at (0.2,0.2): "<<
VectorTools::point_value(dof_handler, solution,Point<2>(0.2, 0.2),vecSol)<<
std::endl;
I still get some errors.
Hermes,
Use another overload. The one returning the solution as a parameter should
work:
https://www.dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#acd358e9b110ccbf4a7f76796d206b9c7
Best,
Daniel
Am Di., 10. Aug. 2021 um 09:41 Uhr schrieb Hermes Sampedro <
hermesampe...@gmail.com>:
Dear all,
It is explained in Step-3 how to evaluate the solution in a point. I am
trying to do the same for Step-29, to evaluate the real and imaginary parts
separately in a single point:
*std::cout << "Solution at (0.2,0.2): "<<
VectorTools::point_value(dof_handler, solution,Point<2>(0.2,
16 matches
Mail list logo