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 and
system_rhs. I am getting the two following solutions from the same
tangent_matrix and system_rhs (data files containing the tangent_matrix and
system_rhs attached).
Matlab solution deal.ii solution
0 0
-0.003033422503555 -0.003033560000000
0 0
-0.002439900952733 -0.002439990000000
0 0
-0.002737307578376 -0.002737430000000
0.001007226333711 0.001007240000000
-0.001069740617723 -0.001069780000000
0.001663965890031 0.001663970000000
-0.002605547185476 -0.002605610000000
0.001113032138953 0.001113040000000
-0.001469379092482 -0.001469430000000
0.000916401048149 0.000916380000000
0 0
0.001674445763014 0.001674440000000
0 0
0.001041269520982 0.001041260000000
-0.000000000000000 -0.000000000000000
0 0
0.003033422503555 0.003033560000000
0 0
0.002439900952733 0.002439990000000
0 0
0.002737307578376 0.002737430000000
0.001007226333711 0.001007240000000
0.001069740617723 0.001069780000000
0.001663965890031 0.001663970000000
0.002605547185476 0.002605610000000
0.001113032138953 0.001113040000000
0.001469379092482 0.001469430000000
I am trying to figure out what could be the reason behind it and how to get
a more accurate solution.
Thank you in advance.
Krish
--
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/ac75c294-7b25-4351-82e5-dd7ad3fd410bn%40googlegroups.com.
0 0 0 -2.1684e-19 0 0 0 0 0.020944 -0.020944 0 0 0 0 0.020944 0 0 0 0 0 0 0 0 0
0 0 0.020944 0.020944 0 0
Component (0,0)
3.30322e+01 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 2.90888e+00 0.00000e+00 1.97804e-01 0.00000e+00 -2.72271e+00
6.51590e-01 -1.21009e+00 1.86168e-01 3.72337e-01 -7.44674e-01 -1.21009e+00
9.30842e-01 0.00000e+00 1.86168e-01 0.00000e+00 -1.48935e+00 5.18945e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 0.00000e+00 3.30322e+01 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 1.97804e-01 0.00000e+00 1.72206e+00 0.00000e+00 -1.32645e+00
-9.30842e-02 2.32711e-02 -5.92481e-01 -5.81776e-01 3.72337e-01 1.86168e-01
4.65421e-02 0.00000e+00 5.75493e-01 0.00000e+00 -3.72337e-01 3.02524e-01
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 3.30322e+01 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 -2.72271e+00 0.00000e+00 -1.32645e+00 0.00000e+00 6.14356e+00
5.58505e-01 4.18879e-01 -5.58505e-01 -9.77384e-01 -7.44674e-01 -3.16486e+00
-8.37758e-01 0.00000e+00 -8.37758e-01 0.00000e+00 1.86168e+00 -2.69944e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 6.51590e-01 0.00000e+00 -9.30842e-02 0.00000e+00 5.58505e-01
5.21272e+01 -1.97339e+01 3.72337e+00 -2.60636e+00 -4.46804e+01 2.15955e+01
1.10770e+01 0.00000e+00 2.11767e+00 0.00000e+00 -1.49400e+01 8.00524e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 -1.21009e+00 0.00000e+00 2.32711e-02 0.00000e+00 4.18879e-01
-1.97339e+01 5.32442e+01 -1.48935e+00 4.00262e+00 1.63828e+01 -4.52389e+01
-3.53720e+00 0.00000e+00 -9.30842e-01 0.00000e+00 7.07440e+00 -2.55516e+01
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 1.86168e-01 0.00000e+00 -5.92481e-01 0.00000e+00 -5.58505e-01
3.72337e+00 -1.48935e+00 2.13759e+01 -5.75447e+00 -2.23402e+01 8.19141e+00
1.07047e+00 0.00000e+00 2.34619e+00 0.00000e+00 -5.16617e+00 1.86168e-01
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 3.72337e-01 0.00000e+00 -5.81776e-01 0.00000e+00 -9.77384e-01
-2.60636e+00 4.00262e+00 -5.75447e+00 1.91754e+01 1.19148e+01 -2.28987e+01
1.86168e-01 0.00000e+00 -2.94518e+00 0.00000e+00 3.72337e-01 9.77384e-01
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 -7.44674e-01 0.00000e+00 3.72337e-01 0.00000e+00 -7.44674e-01
-4.46804e+01 1.63828e+01 -2.23402e+01 1.19148e+01 8.93609e+01 -2.68083e+01
-1.00531e+01 0.00000e+00 -8.65683e+00 0.00000e+00 1.17286e+01 -6.70206e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 -1.21009e+00 0.00000e+00 1.86168e-01 0.00000e+00 -3.16486e+00
2.15955e+01 -4.52389e+01 8.19141e+00 -2.28987e+01 -2.68083e+01 8.48928e+01
1.86168e+00 0.00000e+00 4.09571e+00 0.00000e+00 -8.93609e+00 5.02655e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 9.30842e-01 0.00000e+00 4.65421e-02 0.00000e+00 -8.37758e-01
1.10770e+01 -3.53720e+00 1.07047e+00 1.86168e-01 -1.00531e+01 1.86168e+00
3.48600e+01 0.00000e+00 3.32776e+00 0.00000e+00 -3.92350e+01 0.00000e+00
0.00000e+00 -9.30842e-01 0.00000e+00 -4.65421e-02 0.00000e+00 8.37758e-01
1.10770e+01 3.53720e+00 1.07047e+00 -1.86168e-01 -1.00531e+01 -1.86168e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 3.30322e+01 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 1.86168e-01 0.00000e+00 5.75493e-01 0.00000e+00 -8.37758e-01
2.11767e+00 -9.30842e-01 2.34619e+00 -2.94518e+00 -8.65683e+00 4.09571e+00
3.32776e+00 0.00000e+00 2.36732e+01 0.00000e+00 -1.96873e+01 0.00000e+00
0.00000e+00 -1.86168e-01 0.00000e+00 -5.75493e-01 0.00000e+00 8.37758e-01
2.11767e+00 9.30842e-01 2.34619e+00 2.94518e+00 -8.65683e+00 -4.09571e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 3.30322e+01 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 -1.48935e+00 0.00000e+00 -3.72337e-01 0.00000e+00 1.86168e+00
-1.49400e+01 7.07440e+00 -5.16617e+00 3.72337e-01 1.17286e+01 -8.93609e+00
-3.92350e+01 0.00000e+00 -1.96873e+01 0.00000e+00 7.70737e+01 0.00000e+00
0.00000e+00 1.48935e+00 0.00000e+00 3.72337e-01 0.00000e+00 -1.86168e+00
-1.49400e+01 -7.07440e+00 -5.16617e+00 -3.72337e-01 1.17286e+01 8.93609e+00
0.00000e+00 5.18945e+00 0.00000e+00 3.02524e-01 0.00000e+00 -2.69944e+00
8.00524e+00 -2.55516e+01 1.86168e-01 9.77384e-01 -6.70206e+00 5.02655e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 8.93609e+01
0.00000e+00 5.18945e+00 0.00000e+00 3.02524e-01 0.00000e+00 -2.69944e+00
-8.00524e+00 -2.55516e+01 -1.86168e-01 9.77384e-01 6.70206e+00 5.02655e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
3.30322e+01 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0 0 0 0 0 0
0 0 0 0 0 0
-9.30842e-01 0.00000e+00 -1.86168e-01 0.00000e+00 1.48935e+00 5.18945e+00
0.00000e+00 2.90888e+00 0.00000e+00 1.97804e-01 0.00000e+00 -2.72271e+00
-6.51590e-01 -1.21009e+00 -1.86168e-01 3.72337e-01 7.44674e-01 -1.21009e+00
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 3.30322e+01 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0 0 0 0 0 0
0 0 0 0 0 0
-4.65421e-02 0.00000e+00 -5.75493e-01 0.00000e+00 3.72337e-01 3.02524e-01
0.00000e+00 1.97804e-01 0.00000e+00 1.72206e+00 0.00000e+00 -1.32645e+00
9.30842e-02 2.32711e-02 5.92481e-01 -5.81776e-01 -3.72337e-01 1.86168e-01
0 0 0 0 0 0
0 0 0 0 0 0
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 3.30322e+01 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0 0 0 0 0 0
0 0 0 0 0 0
8.37758e-01 0.00000e+00 8.37758e-01 0.00000e+00 -1.86168e+00 -2.69944e+00
0.00000e+00 -2.72271e+00 0.00000e+00 -1.32645e+00 0.00000e+00 6.14356e+00
-5.58505e-01 4.18879e-01 5.58505e-01 -9.77384e-01 7.44674e-01 -3.16486e+00
0 0 0 0 0 0
0 0 0 0 0 0
1.10770e+01 0.00000e+00 2.11767e+00 0.00000e+00 -1.49400e+01 -8.00524e+00
0.00000e+00 -6.51590e-01 0.00000e+00 9.30842e-02 0.00000e+00 -5.58505e-01
5.21272e+01 1.97339e+01 3.72337e+00 2.60636e+00 -4.46804e+01 -2.15955e+01
0 0 0 0 0 0
0 0 0 0 0 0
3.53720e+00 0.00000e+00 9.30842e-01 0.00000e+00 -7.07440e+00 -2.55516e+01
0.00000e+00 -1.21009e+00 0.00000e+00 2.32711e-02 0.00000e+00 4.18879e-01
1.97339e+01 5.32442e+01 1.48935e+00 4.00262e+00 -1.63828e+01 -4.52389e+01
0 0 0 0 0 0
0 0 0 0 0 0
1.07047e+00 0.00000e+00 2.34619e+00 0.00000e+00 -5.16617e+00 -1.86168e-01
0.00000e+00 -1.86168e-01 0.00000e+00 5.92481e-01 0.00000e+00 5.58505e-01
3.72337e+00 1.48935e+00 2.13759e+01 5.75447e+00 -2.23402e+01 -8.19141e+00
0 0 0 0 0 0
0 0 0 0 0 0
-1.86168e-01 0.00000e+00 2.94518e+00 0.00000e+00 -3.72337e-01 9.77384e-01
0.00000e+00 3.72337e-01 0.00000e+00 -5.81776e-01 0.00000e+00 -9.77384e-01
2.60636e+00 4.00262e+00 5.75447e+00 1.91754e+01 -1.19148e+01 -2.28987e+01
0 0 0 0 0 0
0 0 0 0 0 0
-1.00531e+01 0.00000e+00 -8.65683e+00 0.00000e+00 1.17286e+01 6.70206e+00
0.00000e+00 7.44674e-01 0.00000e+00 -3.72337e-01 0.00000e+00 7.44674e-01
-4.46804e+01 -1.63828e+01 -2.23402e+01 -1.19148e+01 8.93609e+01 2.68083e+01
0 0 0 0 0 0
0 0 0 0 0 0
-1.86168e+00 0.00000e+00 -4.09571e+00 0.00000e+00 8.93609e+00 5.02655e+00
0.00000e+00 -1.21009e+00 0.00000e+00 1.86168e-01 0.00000e+00 -3.16486e+00
-2.15955e+01 -4.52389e+01 -8.19141e+00 -2.28987e+01 2.68083e+01 8.48928e+01
