Repository: ross-rotordynamics/ross
Pull request #834 (2021-08-07 01:11:45)

Correction of 6 dof Gyroscopic matrix and TEST
Hello (@raphaeltimbo )! I have made changes to the gyroscopic matrix and Vergilio helped me checking the tests. I updated the gyroscopic matrix  in the tests and it seems they passed. I hope everthing know works!

marcelosamora: Hello Raphael! Does it still have errors right? As you can see in the picture attached I runned the tests in shaft elements, rotor assembly and test results and the 3 passed. I'm going to review again.
![TESTE ROSS](https://user-images.githubusercontent.com/86723282/128901280-2da0b396-5a28-41b7-b67f-b7b95a2be210.png)

raphaeltimbo: Hi @marcelosamora, sorry for the long time to respond.
Taking as an example the [Tests / build (3.6, ubuntu-latest) (pull_request) ](https://github.com/ross-rotordynamics/ross/pull/834/checks?check_run_id=3284980775), we can see that the tests that are not passing are related to the rubbing calculation (marked as F):

```
ross/tests/test_results.py .........                                     [ 77%]
ross/tests/test_rotor_assembly.py ..............s............s..s...     [ 85%]
ross/tests/test_rubbing.py ..FF                                          [ 86%]
ross/tests/test_shaft_element.py .......................                 [ 92%]
```

Since we know the code in the PR is right, values in the expected array that we have in the test are wrong.
To fix that, you could run the test in debug mode and get the value for `rub.forces_rub[rub.posRUB * 6, :]`, and then replace the expected values that we have in the test with these values.
Let me know if you need any help.
marcelosamora: Hi @raphaeltimbo! No worries! I created a rotor as it is in the test rubbing and Iprinted the rub.forces_rub[rub.posRUB * 6, :]
as you said to me. I updated the forces data but I think there is something wrong with the codes. I tried to check what I did but all the test are not passing as you can see in the picture attached. I  think maybe is something wrong in the ross.
![4](https://user-images.githubusercontent.com/86723282/129496117-740530ab-4015-457d-9b19-9a5ac1e5263f.png)

raphaeltimbo: Hi @marcelosamora !

I am not aware how the test plugin in vscode works, but I think this could be related to some path issue.
Another way that you could try to run the tests is to go to the 'Terminal' within vscode and run the command `pytest`.
To run the command you should be in the tests parent directory (the ross directory with the `__init__.py` file in it).
marcelosamora: Hi @raphaeltimbo! ![Pytest](https://user-images.githubusercontent.com/86723282/129659536-a2dfc651-c2a8-42d4-8402-0f753df91d77.png)
marcelosamora: Hi @raphaeltimbo! I think that now worked! As you can see in the pic the pytest for the test_rubbing is without errors.
![Pytest](https://user-images.githubusercontent.com/86723282/129659595-76aa939f-8f77-434f-b33d-d569c77b592f.png)

codecov[bot]: # [Codecov](https://codecov.io/gh/ross-rotordynamics/ross/pull/834?src=pr&el=h1&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics) Report
> Merging [#834](https://codecov.io/gh/ross-rotordynamics/ross/pull/834?src=pr&el=desc&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics) (347c0c1) into [master](https://codecov.io/gh/ross-rotordynamics/ross/commit/4e60013a4e82797ccd18f9c615439d50d90d3a25?el=desc&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics) (4e60013) will **decrease** coverage by `0.03%`.
> The diff coverage is `100.00%`.

[![Impacted file tree graph](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/graphs/tree.svg?width=650&height=150&src=pr&token=SZjY7TrJnS&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics)](https://codecov.io/gh/ross-rotordynamics/ross/pull/834?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics)

```diff
@@            Coverage Diff             @@
##           master     #834      +/-   ##
==========================================
- Coverage   93.41%   93.37%   -0.04%     
==========================================
  Files          32       31       -1     
  Lines        6541     6503      -38     
==========================================
- Hits         6110     6072      -38     
  Misses        431      431              
```


| [Impacted Files](https://codecov.io/gh/ross-rotordynamics/ross/pull/834?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics) | Coverage Δ | |
|---|---|---|
| [ross/shaft\_element.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9zaGFmdF9lbGVtZW50LnB5) | `88.55% <100.00%> (+0.05%)` | :arrow_up: |
| [ross/element.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9lbGVtZW50LnB5) | `84.78% <0.00%> (-2.01%)` | :arrow_down: |
| [ross/materials.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9tYXRlcmlhbHMucHk=) | `88.75% <0.00%> (-0.67%)` | :arrow_down: |
| [ross/defects/abs\_defect.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9kZWZlY3RzL2Fic19kZWZlY3QucHk=) | `95.65% <0.00%> (-0.10%)` | :arrow_down: |
| [ross/stochastic/st\_results.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9zdG9jaGFzdGljL3N0X3Jlc3VsdHMucHk=) | `85.71% <0.00%> (-0.06%)` | :arrow_down: |
| [ross/defects/rubbing.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9kZWZlY3RzL3J1YmJpbmcucHk=) | `93.67% <0.00%> (-0.04%)` | :arrow_down: |
| [ross/stochastic/st\_rotor\_assembly.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9zdG9jaGFzdGljL3N0X3JvdG9yX2Fzc2VtYmx5LnB5) | `93.27% <0.00%> (-0.04%)` | :arrow_down: |
| [ross/bearing\_seal\_element.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9iZWFyaW5nX3NlYWxfZWxlbWVudC5weQ==) | `94.01% <0.00%> (-0.03%)` | :arrow_down: |
| [ross/results.py](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics#diff-cm9zcy9yZXN1bHRzLnB5) | `94.27% <0.00%> (-0.03%)` | :arrow_down: |
| ... and [3 more](https://codecov.io/gh/ross-rotordynamics/ross/pull/834/diff?src=pr&el=tree-more&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics) | |

------

[Continue to review full report at Codecov](https://codecov.io/gh/ross-rotordynamics/ross/pull/834?src=pr&el=continue&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics).
> **Legend** - [Click here to learn more](https://docs.codecov.io/docs/codecov-delta?utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics)
> `Δ = absolute <relative> (impact)`, `ø = not affected`, `? = missing data`
> Powered by [Codecov](https://codecov.io/gh/ross-rotordynamics/ross/pull/834?src=pr&el=footer&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics). Last update [4e60013...347c0c1](https://codecov.io/gh/ross-rotordynamics/ross/pull/834?src=pr&el=lastupdated&utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics). Read the [comment docs](https://docs.codecov.io/docs/pull-request-comments?utm_medium=referral&utm_source=github&utm_content=comment&utm_campaign=pr+comments&utm_term=ross-rotordynamics).

raphaeltimbo: Hi @marcelosamora ! All tests are passing. The marks showing failures are related to lint:
![image](https://user-images.githubusercontent.com/18506378/130154989-bd9f84f6-eb3c-4103-9ee3-6e9f141cff3c.png)

We use [black](https://ross.readthedocs.io/en/latest/contributing.html#code-style-black) to format our code.
But don't worry about that for now, I will make a commit to format the code using black and will merge the PR.

Thank you so much for fixing this!

Patch:
diff --git a/ross/shaft_element.py b/ross/shaft_element.py
index 4c4a5c30..d999733e 100644
--- a/ross/shaft_element.py
+++ b/ross/shaft_element.py
@@ -35,8 +35,7 @@ class ShaftElement(Element):
 
     Parameters
     ----------
-    L : float, pint.Quantity
-        Element length (m).
+
     idl : float, pint.Quantity
         Inner diameter of the element at the left position (m).
     odl : float, pint.Quantity
@@ -1893,22 +1892,38 @@ class ShaftElement6DoF(ShaftElement):
                     - ((self.idr / 2) ** 4 + (self.idl / 2) ** 4) / 2
                 )
             )
+            tempS = np.pi * (
+                ((self.odr / 2) ** 2 + (self.odl / 2) ** 2) / 2
+                - ((self.idr / 2) ** 2 + (self.idl / 2) ** 2) / 2
+            )
 
             # fmt: off
             # Gyroscopic effect matrix
-            G = (self.material.rho * tempI / (15 * L)) * np.array([
-                [  0, -36, 0,  -3*L,      0, 0,    0,   36, 0,  -3*L,      0, 0],
-                [ 36,   0, 0,     0,   -3*L, 0,  -36,    0, 0,     0,   -3*L, 0],
-                [  0,   0, 0,     0,      0, 0,    0,    0, 0,     0,      0, 0],
-                [3*L,   0, 0,     0,-4*L**2, 0, -3*L,    0, 0,     0,   L**2, 0],
-                [  0, 3*L, 0,4*L**2,      0, 0,    0, -3*L, 0, -L**2,      0, 0],
-                [  0,   0, 0,     0,      0, 0,    0,    0, 0,     0,      0, 0],
-                [  0,  36, 0,   3*L,      0, 0,    0,  -36, 0,   3*L,      0, 0],
-                [-36,   0, 0,     0,    3*L, 0,   36,    0, 0,     0,    3*L, 0],
-                [  0,   0, 0,     0,      0, 0,    0,    0, 0,     0,      0, 0],
-                [3*L,   0, 0,     0,   L**2, 0, -3*L,    0, 0,     0,-4*L**2, 0],
-                [  0, 3*L, 0, -L**2,      0, 0,    0, -3*L, 0,4*L**2,      0, 0],
-                [  0,   0, 0,     0,      0, 0,    0,    0, 0,     0,      0, 0],
+            phi=(
+            12
+            * self.material.E
+            * tempI
+            / (self.material.G_s * self.kappa * tempS * L ** 2)
+            )
+
+            g1 = 36
+            g2 = (3 - 15*phi) * L
+            g3 = (4 + 5*phi + 10*phi**2) * L**2
+            g4 = (-1 - 5*phi + 5*phi**2) * L**2
+
+            G = (self.material.rho * tempI / (15*((1+phi)**2) * L)) * np.array([
+                    [   0,  g1, 0,   -g2,   0, 0,     0, -g1, 0,   -g2,   0, 0],
+                    [ -g1,   0, 0,     0, -g2, 0,    g1,   0, 0,     0, -g2, 0],
+                    [   0,   0, 0,     0,   0, 0,     0,   0, 0,     0,   0, 0],
+                    [  g2,   0, 0,     0,  g3, 0,   -g2,   0, 0,     0,  g4, 0],
+                    [   0,  g2, 0,   -g3,   0, 0,     0, -g2, 0,   -g4,   0, 0],
+                    [   0,   0, 0,     0,   0, 0,     0,   0, 0,     0,   0, 0],
+                    [   0, -g1, 0,    g2,   0, 0,     0,  g1, 0,    g2,   0, 0],
+                    [  g1,   0, 0,     0,  g2, 0,   -g1,   0, 0,     0,  g2, 0],
+                    [   0,   0, 0,     0,   0, 0,     0,   0, 0,     0,   0, 0],
+                    [  g2,   0, 0,     0,  g4, 0,   -g2,   0, 0,     0,  g3, 0],
+                    [   0,  g2, 0,   -g4,   0, 0,     0, -g2, 0,   -g3,   0, 0],
+                    [   0,   0, 0,     0,   0, 0,     0,   0, 0,     0,   0, 0],
             ])
             # fmt: on
         else:


Test patch:
diff --git a/ross/tests/test_rubbing.py b/ross/tests/test_rubbing.py
index f902e32c..d0bc6c89 100644
--- a/ross/tests/test_rubbing.py
+++ b/ross/tests/test_rubbing.py
@@ -146,132 +146,132 @@ def test_rub_parameters_units(rub_units):
 def test_rub_forces(rub):
     assert rub.forces_rub[rub.posRUB * 6, :] == pytest.approx(
         # fmt: off
-    np.array([  0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         1.33959978,   2.38449456,   2.49659676,   1.81196092,
-         0.59693967,  -1.62826881,  -4.24183226,  -6.00328692,
-        -6.20469077,  -4.73542742,  -2.72934246,  -2.58572889,
-        -6.09380802, -11.89929423, -16.28700512, -16.40580808,
-       -12.28949661,  -6.4715516 ,  -1.9219398 ,   0.03597126,
-         0.1483467 ,   0.06399221,   0.18405941,   0.44370982,
-         0.59170921,   0.55028863,   0.61661081,   1.15972902,
-         1.99145991,   2.2827491 ,   1.21484512,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,  -0.16254596,
-        -1.15328767,  -1.31402774,  -0.57077859,   0.17589894,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         1.35177422,   4.63052083,   8.26595768,  10.77229714,
-        11.33408396,   9.97185414,   7.43725544,   4.8091569 ,
-         2.90279915,   1.9280518 ,   1.60783411,   1.51278525,
-         1.30371098,   0.82428403,   0.11364267,  -1.19150903,
-        -2.28684891,  -2.8371608 ,  -2.70293006,  -1.7191864 ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-        -0.36389397,  -0.02296953,  -0.58251754,  -1.26899237,
-        -1.37833025,  -0.83533185,   0.        ,  -1.54209802,
-        -3.90213044,  -5.84909057,  -5.54764546,  -2.82245046,
-         0.        ,   0.        ,  -0.9813435 ,  -5.02620683,
-        -7.74948602,  -7.28498629,  -4.21292376,  -0.68113248,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.49651344,   0.94005508,   1.20771384,   1.01602262,
-         0.58504824,   0.42249176,   0.74712263,   1.21263464,
-         1.24036134,   0.53226162,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.81085807,   2.24845489,
-         2.79783564,   2.36749471,   1.48390821,   0.91606906,
-         1.10872277,   1.79762142,   2.27408108,   2.02409761,
-         1.09722613,   0.03554649,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,  -3.20746741,
-        -5.67109902,  -5.84109266,  -3.27138013,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,  -0.94605595,
-        -2.89024152,  -2.82355527,  -0.91081202,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.29208245,   0.31276202,   0.55407755,
-         0.65238154,   0.38599931,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.30206142,
-         1.24645948,   2.17118657,   2.79095166,   2.92291644,
-         2.4997066 ,   1.57624218,   0.3448333 ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-        -0.25292217,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,  -1.16895344,  -3.13003066,
-        -4.31801297,  -3.79661585,  -1.72270664,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.22784948,   0.25813886,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,  -1.0783812 ,  -1.3305057 ,
-         0.04528174,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ])
+    np.array([ 0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        1.34106430e+00,  2.38612745e+00,  2.49817920e+00,  1.81347173e+00,
+        5.98445473e-01, -1.62547004e+00, -4.23917250e+00, -6.00118806e+00,
+       -6.20373067e+00, -4.73612273e+00, -2.73184699e+00, -2.58971885e+00,
+       -6.09874812e+00, -1.19048343e+01, -1.62930120e+01, -1.64120039e+01,
+       -1.22951018e+01, -6.47561269e+00, -1.92419235e+00,  3.53374623e-02,
+        1.47810258e-01,  6.34772882e-02,  1.83806641e-01,  4.43954908e-01,
+        5.92533172e-01,  5.51738633e-01,  6.18865825e-01,  1.16291085e+00,
+        1.99527847e+00,  2.28657481e+00,  1.21810729e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.55680614e-01,
+       -1.15087399e+00, -1.31570879e+00, -5.73686616e-01,  1.75075759e-01,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        1.35857558e+00,  4.63547122e+00,  8.26625410e+00,  1.07680988e+01,
+        1.13286366e+01,  9.96961037e+00,  7.44128958e+00,  4.81928761e+00,
+        2.91576005e+00,  1.93968467e+00,  1.61560462e+00,  1.51666598e+00,
+        1.30530274e+00,  8.25200469e-01,  1.14260074e-01, -1.19244337e+00,
+       -2.29090277e+00, -2.84242534e+00, -2.70568998e+00, -1.71797718e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+       -3.62954221e-01, -2.15664186e-02, -5.82573814e-01, -1.27145991e+00,
+       -1.38207772e+00, -8.36333007e-01,  0.00000000e+00, -1.52638197e+00,
+       -3.88459875e+00, -5.84147533e+00, -5.55533169e+00, -2.84156510e+00,
+        0.00000000e+00,  0.00000000e+00, -9.98460245e-01, -5.03924706e+00,
+       -7.75632048e+00, -7.28589841e+00, -4.21050757e+00, -6.78335856e-01,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        4.97669123e-01,  9.43108661e-01,  1.21217301e+00,  1.02059004e+00,
+        5.88846364e-01,  4.25532247e-01,  7.49102363e-01,  1.21248611e+00,
+        1.23719949e+00,  5.26714086e-01,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  8.30579340e-01,  2.26273248e+00,
+        2.80172735e+00,  2.36060364e+00,  1.47121281e+00,  9.05484679e-01,
+        1.10607178e+00,  1.80305085e+00,  2.28352354e+00,  2.03304519e+00,
+        1.10298146e+00,  3.77694650e-02,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -3.22380450e+00,
+       -5.68021237e+00, -5.83462167e+00, -3.25195544e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -9.90196556e-01,
+       -2.90885128e+00, -2.81169132e+00, -8.81117993e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  2.69024403e-01,  3.11412852e-01,  5.58798236e-01,
+        6.40371175e-01,  3.49164726e-01,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  2.36971006e-01,
+        1.21346945e+00,  2.18249083e+00,  2.84208082e+00,  2.99156571e+00,
+        2.55421246e+00,  1.59111859e+00,  3.14246711e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+       -9.16343507e-02,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00, -1.16296507e+00, -3.23286551e+00,
+       -4.51354090e+00, -3.98257895e+00, -1.81166144e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  2.12219345e-01,  2.21171647e-01,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00, -9.56292801e-01, -1.33271793e+00,
+       -1.03821500e-02,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00])
         # fmt: on
     )
 
@@ -283,14 +283,14 @@ def test_rub_forces(rub):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-       -2.04052219e+00, -4.37686754e+00, -6.57561714e+00, -8.28614638e+00,
-       -9.37724617e+00, -9.76009440e+00, -9.25692987e+00, -7.68859694e+00,
-       -5.19404219e+00, -2.47936767e+00, -6.07948333e-01, -1.39115180e-01,
-       -2.90615182e-01,  5.97112915e-01,  3.77091904e+00,  7.17958236e+00,
-        8.50566481e+00,  6.92594573e+00,  3.76719699e+00,  1.38672233e+00,
-        1.14491558e+00,  2.47008298e+00,  3.68895320e+00,  3.55463037e+00,
-        2.15169896e+00,  7.16172655e-01,  4.58130638e-01,  1.31722645e+00,
-        2.08017814e+00,  1.73698314e+00,  3.66421390e-01,  0.00000000e+00,
+       -2.04165157e+00, -4.37810280e+00, -6.57666995e+00, -8.28694141e+00,
+       -9.37791126e+00, -9.76081035e+00, -9.25783043e+00, -7.68973505e+00,
+       -5.19539491e+00, -2.48084518e+00, -6.09414243e-01, -1.40464287e-01,
+       -2.91868675e-01,  5.94864394e-01,  3.76903307e+00,  7.17867399e+00,
+        8.50615176e+00,  6.92755001e+00,  3.76919148e+00,  1.38886654e+00,
+        1.14793265e+00,  2.47476310e+00,  3.69518317e+00,  3.56151341e+00,
+        2.15845164e+00,  7.22745269e-01,  4.64776971e-01,  1.32361970e+00,
+        2.08542187e+00,  1.74043906e+00,  3.68327123e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -301,14 +301,14 @@ def test_rub_forces(rub):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.41608952e+00,
-        3.91872725e+00,  6.08751103e+00,  6.04879589e+00,  3.34005262e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.40526270e+00,
+        3.91183574e+00,  6.08879987e+00,  6.05701802e+00,  3.35127909e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        1.32257678e+00,  4.72982449e+00,  5.65624377e+00,  3.74610313e+00,
-        2.79267497e-01, -1.61571947e+00, -2.61760973e+00, -2.76295228e+00,
-       -2.46754501e+00, -2.27912450e+00, -2.49834149e+00, -3.08123496e+00,
-       -3.77081696e+00, -4.26182369e+00, -4.32814157e+00, -3.92465771e+00,
-       -3.18784204e+00, -2.30089569e+00, -1.36585055e+00, -4.14179880e-01,
+        1.33684009e+00,  4.74117308e+00,  5.66393572e+00,  3.75144548e+00,
+        2.83459537e-01, -1.61445732e+00, -2.61883874e+00, -2.76818696e+00,
+       -2.47679180e+00, -2.29014373e+00, -2.50794047e+00, -3.08754938e+00,
+       -3.77455164e+00, -4.26560645e+00, -4.33448168e+00, -3.93383052e+00,
+       -3.19736872e+00, -2.30738448e+00, -1.36778475e+00, -4.13019001e-01,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -321,15 +321,15 @@ def test_rub_forces(rub):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-       -2.27746914e-01, -2.70148968e+00, -3.76415772e+00, -3.18400597e+00,
-       -1.61071020e+00, -1.55937047e-01,  0.00000000e+00, -3.48135704e-01,
-       -1.20173785e+00, -1.41164098e+00, -8.22789836e-01,  2.22972608e-02,
-        0.00000000e+00,  0.00000000e+00,  1.11826215e+00,  3.42190067e+00,
-        6.43880255e+00,  7.73939551e+00,  6.02642177e+00,  2.15551882e+00,
+       -2.30731635e-01, -2.71055349e+00, -3.77777456e+00, -3.19820518e+00,
+       -1.62052825e+00, -1.57786580e-01,  0.00000000e+00, -3.39194075e-01,
+       -1.19486567e+00, -1.40982868e+00, -8.24731765e-01,  1.77245380e-02,
+        0.00000000e+00,  0.00000000e+00,  1.12393272e+00,  3.42883960e+00,
+        6.44392419e+00,  7.74026022e+00,  6.02300041e+00,  2.15021680e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        1.30262708e+00,  2.86790215e+00,  2.69095371e+00,  1.30596104e+00,
-        6.27355622e-02, -1.46409693e-01, -1.61596825e-03,  2.02368230e-01,
-        7.76880971e-02, -1.13943983e-01,  0.00000000e+00,  0.00000000e+00,
+        1.31206244e+00,  2.88224433e+00,  2.70442915e+00,  1.31540273e+00,
+        6.83579157e-02, -1.44657955e-01, -9.74069771e-04,  2.01609711e-01,
+        7.58795510e-02, -1.14712672e-01,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -342,12 +342,12 @@ def test_rub_forces(rub):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00, -1.20894403e-01, -9.67454370e-02,
-       -3.75553448e-01, -6.69880393e-01, -7.81524379e-01, -8.20393160e-01,
-       -1.19541101e+00, -2.17352648e+00, -3.38963807e+00, -3.97901667e+00,
-       -3.24152409e+00, -1.21249778e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.71167570e+00,
-       -2.63258535e+00, -2.11381589e+00, -7.29071480e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00, -1.16918168e-01, -9.48748823e-02,
+       -3.74701248e-01, -6.68192840e-01, -7.77720853e-01, -8.15718667e-01,
+       -1.19404470e+00, -2.17926682e+00, -3.40215785e+00, -3.99418641e+00,
+       -3.25449483e+00, -1.22032570e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.72304245e+00,
+       -2.63776333e+00, -2.11110151e+00, -7.22241992e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -360,15 +360,15 @@ def test_rub_forces(rub):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.55809454e-01,
-       -1.33550998e-01,  1.30622371e-01,  5.14772296e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.37988828e-01,
+       -1.36497601e-01,  1.32278945e-01,  5.15571508e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  5.87400070e-01,  2.52352092e+00,  3.18829790e+00,
-        2.12154931e+00,  1.96463435e-01,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.71590151e-01,
-        8.49398637e-02, -3.92130275e-03, -3.10907319e-01, -7.32520442e-01,
-       -1.04949045e+00, -1.03758926e+00, -5.57304261e-01,  0.00000000e+00,
+        0.00000000e+00,  8.42548433e-01,  2.72091640e+00,  3.24542027e+00,
+        2.03399132e+00,  3.77699186e-02,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.99799699e-01,
+        5.92260640e-02, -1.02272668e-02, -3.19939917e-01, -7.49813663e-01,
+       -1.06945353e+00, -1.04483902e+00, -5.36963585e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -382,15 +382,15 @@ def test_rub_forces(rub):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-       -2.22986515e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+       -3.23192576e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  1.36141783e-01, -1.62089253e-01,
-       -6.56719131e-02,  3.67462355e-01,  6.11527272e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  1.42596962e-01, -1.71901114e-01,
+       -7.45893810e-02,  3.71940362e-01,  6.24286373e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00, -7.23117550e-02, -4.63871628e-02,
+        0.00000000e+00,  0.00000000e+00, -9.16154067e-02, -8.72890060e-02,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -399,8 +399,8 @@ def test_rub_forces(rub):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  8.90041966e-01,  1.14524591e+00,
-        4.21896347e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  8.36067934e-01,  1.13854180e+00,
+        4.79721529e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00])
         # fmt: on
@@ -410,132 +410,132 @@ def test_rub_forces(rub):
 def test_rub_forces_units(rub_units):
     assert rub_units.forces_rub[rub_units.posRUB * 6, :] == pytest.approx(
         # fmt: off
-    np.array([  0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         1.33959978,   2.38449456,   2.49659676,   1.81196092,
-         0.59693967,  -1.62826881,  -4.24183226,  -6.00328692,
-        -6.20469077,  -4.73542742,  -2.72934246,  -2.58572889,
-        -6.09380802, -11.89929423, -16.28700512, -16.40580808,
-       -12.28949661,  -6.4715516 ,  -1.9219398 ,   0.03597126,
-         0.1483467 ,   0.06399221,   0.18405941,   0.44370982,
-         0.59170921,   0.55028863,   0.61661081,   1.15972902,
-         1.99145991,   2.2827491 ,   1.21484512,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,  -0.16254596,
-        -1.15328767,  -1.31402774,  -0.57077859,   0.17589894,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         1.35177422,   4.63052083,   8.26595768,  10.77229714,
-        11.33408396,   9.97185414,   7.43725544,   4.8091569 ,
-         2.90279915,   1.9280518 ,   1.60783411,   1.51278525,
-         1.30371098,   0.82428403,   0.11364267,  -1.19150903,
-        -2.28684891,  -2.8371608 ,  -2.70293006,  -1.7191864 ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-        -0.36389397,  -0.02296953,  -0.58251754,  -1.26899237,
-        -1.37833025,  -0.83533185,   0.        ,  -1.54209802,
-        -3.90213044,  -5.84909057,  -5.54764546,  -2.82245046,
-         0.        ,   0.        ,  -0.9813435 ,  -5.02620683,
-        -7.74948602,  -7.28498629,  -4.21292376,  -0.68113248,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.49651344,   0.94005508,   1.20771384,   1.01602262,
-         0.58504824,   0.42249176,   0.74712263,   1.21263464,
-         1.24036134,   0.53226162,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.81085807,   2.24845489,
-         2.79783564,   2.36749471,   1.48390821,   0.91606906,
-         1.10872277,   1.79762142,   2.27408108,   2.02409761,
-         1.09722613,   0.03554649,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,  -3.20746741,
-        -5.67109902,  -5.84109266,  -3.27138013,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,  -0.94605595,
-        -2.89024152,  -2.82355527,  -0.91081202,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.29208245,   0.31276202,   0.55407755,
-         0.65238154,   0.38599931,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.30206142,
-         1.24645948,   2.17118657,   2.79095166,   2.92291644,
-         2.4997066 ,   1.57624218,   0.3448333 ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-        -0.25292217,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,  -1.16895344,  -3.13003066,
-        -4.31801297,  -3.79661585,  -1.72270664,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.22784948,   0.25813886,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,  -1.0783812 ,  -1.3305057 ,
-         0.04528174,   0.        ,   0.        ,   0.        ,
-         0.        ,   0.        ,   0.        ,   0.        ,
-         0.        ])
+    np.array([ 0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        1.34106430e+00,  2.38612745e+00,  2.49817920e+00,  1.81347173e+00,
+        5.98445473e-01, -1.62547004e+00, -4.23917250e+00, -6.00118806e+00,
+       -6.20373067e+00, -4.73612273e+00, -2.73184699e+00, -2.58971885e+00,
+       -6.09874812e+00, -1.19048343e+01, -1.62930120e+01, -1.64120039e+01,
+       -1.22951018e+01, -6.47561269e+00, -1.92419235e+00,  3.53374623e-02,
+        1.47810258e-01,  6.34772882e-02,  1.83806641e-01,  4.43954908e-01,
+        5.92533172e-01,  5.51738633e-01,  6.18865825e-01,  1.16291085e+00,
+        1.99527847e+00,  2.28657481e+00,  1.21810729e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.55680614e-01,
+       -1.15087399e+00, -1.31570879e+00, -5.73686616e-01,  1.75075759e-01,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        1.35857558e+00,  4.63547122e+00,  8.26625410e+00,  1.07680988e+01,
+        1.13286366e+01,  9.96961037e+00,  7.44128958e+00,  4.81928761e+00,
+        2.91576005e+00,  1.93968467e+00,  1.61560462e+00,  1.51666598e+00,
+        1.30530274e+00,  8.25200469e-01,  1.14260074e-01, -1.19244337e+00,
+       -2.29090277e+00, -2.84242534e+00, -2.70568998e+00, -1.71797718e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+       -3.62954221e-01, -2.15664186e-02, -5.82573814e-01, -1.27145991e+00,
+       -1.38207772e+00, -8.36333007e-01,  0.00000000e+00, -1.52638197e+00,
+       -3.88459875e+00, -5.84147533e+00, -5.55533169e+00, -2.84156510e+00,
+        0.00000000e+00,  0.00000000e+00, -9.98460245e-01, -5.03924706e+00,
+       -7.75632048e+00, -7.28589841e+00, -4.21050757e+00, -6.78335856e-01,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        4.97669123e-01,  9.43108661e-01,  1.21217301e+00,  1.02059004e+00,
+        5.88846364e-01,  4.25532247e-01,  7.49102363e-01,  1.21248611e+00,
+        1.23719949e+00,  5.26714086e-01,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  8.30579340e-01,  2.26273248e+00,
+        2.80172735e+00,  2.36060364e+00,  1.47121281e+00,  9.05484679e-01,
+        1.10607178e+00,  1.80305085e+00,  2.28352354e+00,  2.03304519e+00,
+        1.10298146e+00,  3.77694650e-02,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -3.22380450e+00,
+       -5.68021237e+00, -5.83462167e+00, -3.25195544e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -9.90196556e-01,
+       -2.90885128e+00, -2.81169132e+00, -8.81117993e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  2.69024403e-01,  3.11412852e-01,  5.58798236e-01,
+        6.40371175e-01,  3.49164726e-01,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  2.36971006e-01,
+        1.21346945e+00,  2.18249083e+00,  2.84208082e+00,  2.99156571e+00,
+        2.55421246e+00,  1.59111859e+00,  3.14246711e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+       -9.16343507e-02,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00, -1.16296507e+00, -3.23286551e+00,
+       -4.51354090e+00, -3.98257895e+00, -1.81166144e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  2.12219345e-01,  2.21171647e-01,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00, -9.56292801e-01, -1.33271793e+00,
+       -1.03821500e-02,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00])
         # fmt: on
     )
 
@@ -547,14 +547,14 @@ def test_rub_forces_units(rub_units):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-       -2.04052219e+00, -4.37686754e+00, -6.57561714e+00, -8.28614638e+00,
-       -9.37724617e+00, -9.76009440e+00, -9.25692987e+00, -7.68859694e+00,
-       -5.19404219e+00, -2.47936767e+00, -6.07948333e-01, -1.39115180e-01,
-       -2.90615182e-01,  5.97112915e-01,  3.77091904e+00,  7.17958236e+00,
-        8.50566481e+00,  6.92594573e+00,  3.76719699e+00,  1.38672233e+00,
-        1.14491558e+00,  2.47008298e+00,  3.68895320e+00,  3.55463037e+00,
-        2.15169896e+00,  7.16172655e-01,  4.58130638e-01,  1.31722645e+00,
-        2.08017814e+00,  1.73698314e+00,  3.66421390e-01,  0.00000000e+00,
+       -2.04165157e+00, -4.37810280e+00, -6.57666995e+00, -8.28694141e+00,
+       -9.37791126e+00, -9.76081035e+00, -9.25783043e+00, -7.68973505e+00,
+       -5.19539491e+00, -2.48084518e+00, -6.09414243e-01, -1.40464287e-01,
+       -2.91868675e-01,  5.94864394e-01,  3.76903307e+00,  7.17867399e+00,
+        8.50615176e+00,  6.92755001e+00,  3.76919148e+00,  1.38886654e+00,
+        1.14793265e+00,  2.47476310e+00,  3.69518317e+00,  3.56151341e+00,
+        2.15845164e+00,  7.22745269e-01,  4.64776971e-01,  1.32361970e+00,
+        2.08542187e+00,  1.74043906e+00,  3.68327123e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -565,14 +565,14 @@ def test_rub_forces_units(rub_units):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.41608952e+00,
-        3.91872725e+00,  6.08751103e+00,  6.04879589e+00,  3.34005262e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.40526270e+00,
+        3.91183574e+00,  6.08879987e+00,  6.05701802e+00,  3.35127909e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        1.32257678e+00,  4.72982449e+00,  5.65624377e+00,  3.74610313e+00,
-        2.79267497e-01, -1.61571947e+00, -2.61760973e+00, -2.76295228e+00,
-       -2.46754501e+00, -2.27912450e+00, -2.49834149e+00, -3.08123496e+00,
-       -3.77081696e+00, -4.26182369e+00, -4.32814157e+00, -3.92465771e+00,
-       -3.18784204e+00, -2.30089569e+00, -1.36585055e+00, -4.14179880e-01,
+        1.33684009e+00,  4.74117308e+00,  5.66393572e+00,  3.75144548e+00,
+        2.83459537e-01, -1.61445732e+00, -2.61883874e+00, -2.76818696e+00,
+       -2.47679180e+00, -2.29014373e+00, -2.50794047e+00, -3.08754938e+00,
+       -3.77455164e+00, -4.26560645e+00, -4.33448168e+00, -3.93383052e+00,
+       -3.19736872e+00, -2.30738448e+00, -1.36778475e+00, -4.13019001e-01,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -585,15 +585,15 @@ def test_rub_forces_units(rub_units):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-       -2.27746914e-01, -2.70148968e+00, -3.76415772e+00, -3.18400597e+00,
-       -1.61071020e+00, -1.55937047e-01,  0.00000000e+00, -3.48135704e-01,
-       -1.20173785e+00, -1.41164098e+00, -8.22789836e-01,  2.22972608e-02,
-        0.00000000e+00,  0.00000000e+00,  1.11826215e+00,  3.42190067e+00,
-        6.43880255e+00,  7.73939551e+00,  6.02642177e+00,  2.15551882e+00,
+       -2.30731635e-01, -2.71055349e+00, -3.77777456e+00, -3.19820518e+00,
+       -1.62052825e+00, -1.57786580e-01,  0.00000000e+00, -3.39194075e-01,
+       -1.19486567e+00, -1.40982868e+00, -8.24731765e-01,  1.77245380e-02,
+        0.00000000e+00,  0.00000000e+00,  1.12393272e+00,  3.42883960e+00,
+        6.44392419e+00,  7.74026022e+00,  6.02300041e+00,  2.15021680e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        1.30262708e+00,  2.86790215e+00,  2.69095371e+00,  1.30596104e+00,
-        6.27355622e-02, -1.46409693e-01, -1.61596825e-03,  2.02368230e-01,
-        7.76880971e-02, -1.13943983e-01,  0.00000000e+00,  0.00000000e+00,
+        1.31206244e+00,  2.88224433e+00,  2.70442915e+00,  1.31540273e+00,
+        6.83579157e-02, -1.44657955e-01, -9.74069771e-04,  2.01609711e-01,
+        7.58795510e-02, -1.14712672e-01,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -606,12 +606,12 @@ def test_rub_forces_units(rub_units):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00, -1.20894403e-01, -9.67454370e-02,
-       -3.75553448e-01, -6.69880393e-01, -7.81524379e-01, -8.20393160e-01,
-       -1.19541101e+00, -2.17352648e+00, -3.38963807e+00, -3.97901667e+00,
-       -3.24152409e+00, -1.21249778e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.71167570e+00,
-       -2.63258535e+00, -2.11381589e+00, -7.29071480e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00, -1.16918168e-01, -9.48748823e-02,
+       -3.74701248e-01, -6.68192840e-01, -7.77720853e-01, -8.15718667e-01,
+       -1.19404470e+00, -2.17926682e+00, -3.40215785e+00, -3.99418641e+00,
+       -3.25449483e+00, -1.22032570e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.72304245e+00,
+       -2.63776333e+00, -2.11110151e+00, -7.22241992e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -624,15 +624,15 @@ def test_rub_forces_units(rub_units):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.55809454e-01,
-       -1.33550998e-01,  1.30622371e-01,  5.14772296e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  1.37988828e-01,
+       -1.36497601e-01,  1.32278945e-01,  5.15571508e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  5.87400070e-01,  2.52352092e+00,  3.18829790e+00,
-        2.12154931e+00,  1.96463435e-01,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.71590151e-01,
-        8.49398637e-02, -3.92130275e-03, -3.10907319e-01, -7.32520442e-01,
-       -1.04949045e+00, -1.03758926e+00, -5.57304261e-01,  0.00000000e+00,
+        0.00000000e+00,  8.42548433e-01,  2.72091640e+00,  3.24542027e+00,
+        2.03399132e+00,  3.77699186e-02,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  0.00000000e+00, -1.99799699e-01,
+        5.92260640e-02, -1.02272668e-02, -3.19939917e-01, -7.49813663e-01,
+       -1.06945353e+00, -1.04483902e+00, -5.36963585e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -646,15 +646,15 @@ def test_rub_forces_units(rub_units):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-       -2.22986515e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+       -3.23192576e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  1.36141783e-01, -1.62089253e-01,
-       -6.56719131e-02,  3.67462355e-01,  6.11527272e-01,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  1.42596962e-01, -1.71901114e-01,
+       -7.45893810e-02,  3.71940362e-01,  6.24286373e-01,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00, -7.23117550e-02, -4.63871628e-02,
+        0.00000000e+00,  0.00000000e+00, -9.16154067e-02, -8.72890060e-02,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
@@ -663,8 +663,8 @@ def test_rub_forces_units(rub_units):
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
-        0.00000000e+00,  0.00000000e+00,  8.90041966e-01,  1.14524591e+00,
-        4.21896347e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
+        0.00000000e+00,  0.00000000e+00,  8.36067934e-01,  1.13854180e+00,
+        4.79721529e-01,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00,  0.00000000e+00,  0.00000000e+00,  0.00000000e+00,
         0.00000000e+00])
         # fmt: on
diff --git a/ross/tests/test_shaft_element.py b/ross/tests/test_shaft_element.py
index 7a239ee9..8f0fa2d3 100644
--- a/ross/tests/test_shaft_element.py
+++ b/ross/tests/test_shaft_element.py
@@ -382,19 +382,20 @@ def test_s6_M_matrix(s6_eb):
 
 def test_s6_G_matrix(s6_eb):
     # fmt: off
-    G_mat_6DoF = 0.01 * np.array(
-        [[           0,  -9.7239e+10,            0,  -8.1032e+08,            0,            0,            0,   9.7239e+10,            0,  -8.1032e+08,            0,            0],
-         [  9.7239e+10,            0,            0,            0,  -8.1032e+08,            0,  -9.7239e+10,            0,            0,            0,  -8.1032e+08,            0],
+    G_mat_6DoF = np.array(
+        [[           0,   3.3495e-02,            0,   2.3779e+02,            0,            0,            0,   -3.3495e-02,           0,   2.3779e+02,            0,            0],
+         [ -3.3495e-02,            0,            0,            0,   2.3779e+02,            0,   3.3495e-02,            0,            0,            0,   2.3779e+02,            0],
          [           0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0],
-         [  8.1032e+08,            0,            0,            0,  -1.0804e+08,            0,  -8.1032e+08,            0,            0,            0,   2.7011e+07,            0],
-         [           0,   8.1032e+08,            0,   1.0804e+08,            0,            0,            0,  -8.1032e+08,            0,  -2.7011e+07,            0,            0],
+         [-2.37798e+02,            0,            0,            0,   2.7010e+06,            0,   2.3779e+02,            0,            0,            0,   1.3505e+06,            0],
+         [           0,  -2.3779e+02,            0,  -2.7010e+06,            0,            0,            0,   2.3779e+02,            0,  -1.3505e+06,            0,            0],
          [           0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0],
-         [           0,   9.7239e+10,            0,   8.1032e+08,            0,            0,            0,  -9.7239e+10,            0,   8.1032e+08,            0,            0],
-         [ -9.7239e+10,            0,            0,            0,   8.1032e+08,            0,   9.7239e+10,            0,            0,            0,   8.1032e+08,            0],
+         [           0,  -3.3495e-02,            0,  -2.3779e+02,            0,            0,            0,   3.3495e-02,            0,  -2.3779e+02,            0,            0],
+         [  3.3495e-02,            0,            0,            0,  -2.3779e+02,            0,  -3.3495e-02,            0,            0,            0,  -2.3779e+02,            0],
          [           0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0],
-         [  8.1032e+08,            0,            0,            0,   2.7011e+07,            0,  -8.1032e+08,            0,            0,            0,  -1.0804e+08,            0],
-         [           0,   8.1032e+08,            0,  -2.7011e+07,            0,            0,            0,  -8.1032e+08,            0,   1.0804e+08,            0,            0],
+         [ -2.3779e+02,            0,            0,            0,   1.3505e+06,            0,   2.3779e+02,            0,            0,            0,   2.7010e+06,            0],
+         [           0,  -2.3779e+02,            0,  -1.3505e+06,            0,            0,            0,   2.3779e+02,            0,  -2.7010e+06,            0,            0],
          [           0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0,            0]])
+
     # fmt: on
     assert_allclose(s6_eb.G(), G_mat_6DoF, rtol=1e-3)