Hello Simon,
Thank you very much for the suggestion. The Matlab and deal.ii solutions
matched with each other after increasing the number of digits written to
the std::ofstream object.
Thanks,
Krish
On Saturday, 29 July 2023 at 20:13:48 UTC+5:30 simon...@gmail.com wrote:
> Hi.
>
> From your
Hi.
From your attached files it looks like you have saved the matrix and rhs with
quite few digits. Do you get the same result if you increase the number of
digits you write to your std::ofstream object?
Compare:
https://en.cppreference.com/w/cpp/io/manip/setprecision
Best,
Simon
On
Hello,
I am solving a problem that has 30 degrees of freedom. The newton update
solution obtained from deal.ii (using both Direct and CG solver) seems to
be less accurate (differs after 4 significant digits) from the newton
update solution obtained from Matlab using the same tangent_matrix
