add 149k test

models/dolicorhynchops/124k/README.md

-# dolicorhynchops 124k
-
-Dolicorhynchops osborni FHSM VP404
-
-Intermediate model generated by RB with MeshLab.
-
-![screenshot](screenshot.png)
\ No newline at end of file

models/dolicorhynchops/149k/LTMJ.ply

+ply
+format binary_little_endian 1.0
+comment VCGLIB generated
+element vertex 1
+property float x
+property float y
+property float z
+property float nx
+property float ny
+property float nz
+element face 0
+property list uchar int vertex_indices
+end_header
+N?O²dCxtd^?@
\ No newline at end of file

models/dolicorhynchops/149k/LmuscleF.off

+NOFF
+1 0 0
+-100.1458893 -173.13895 227.3909

models/dolicorhynchops/149k/README.md

+# dolicorhynchops 149k
+
+Dolicorhynchops osborni FHSM VP404
+
+Intermediate model generated by RB with MeshLab.
+
+![screenshot](screenshot.png)
+
+
+First Results
+
+```
+User CPU : 2m15.11s
+Real CPU : 2m14.13s
+ 
+[TSC-STP]                         Number of steps : 1
+[TSC-ITE]              Number of mech. iterations : 4
+[TSC-INW]                         Internal energy : 24.5869
+[TSC-EXW]                 Work of external forces : 24.587
+[TSC-CPU]                           User CPU Time : 135.109
+[TSC-REA]                           Real CPU Time : 134.127
+[TSC-KER]                         Kernel CPU Time : 8.01562
+[TSC-MEM]                        Peak Memory [Kb] : 5.12406e+06
+[TSC-EXT]                                 LTMJ_Fx : 59.4756
+[TSC-EXT]                                 LTMJ_Fy : 45.4345
+[TSC-EXT]                                 LTMJ_Fz : -34.9242
+[TSC-EXT]                                 RTMJ_Fx : 47.4809
+[TSC-EXT]                                 RTMJ_Fy : -49.7536
+[TSC-EXT]                                 RTMJ_Fz : -34.8951
+[TSC-EXT]                                teeth_Fx : 21.3817
+[TSC-EXT]                                teeth_Fy : 0
+[TSC-EXT]                                teeth_Fz : 0
+[TSC-EXT]                          Lmuscle_lfs_Fx : -65.77
+[TSC-EXT]                          Lmuscle_lfs_Fy : -66.0663
+[TSC-EXT]                          Lmuscle_lfs_Fz : 36.1878
+[TSC-EXT]                          Rmuscle_lfs_Fx : -62.5683
+[TSC-EXT]                          Rmuscle_lfs_Fy : 70.3856
+[TSC-EXT]                          Rmuscle_lfs_Fz : 33.6316
+```
+
+![screenshot](results.png)

models/dolicorhynchops/149k/RTMJ.ply

+ply
+format binary_little_endian 1.0
+comment VCGLIB generated
+element vertex 1
+property float x
+property float y
+property float z
+property float nx
+property float ny
+property float nz
+element face 0
+property list uchar int vertex_indices
+end_header
+'|0nٍBqC#v$<@
\ No newline at end of file

models/dolicorhynchops/149k/RmuscleF.off

+NOFF
+1 0 0
+-108.716309 179.13171 233.8385

models/dolicorhynchops/149k/teeth.ply

+ply
+format binary_little_endian 1.0
+comment VCGLIB generated
+element vertex 2
+property float x
+property float y
+property float z
+property float nx
+property float ny
+property float nz
+element face 0
+property list uchar int vertex_indices
+end_header
+'*A|idz}#K?ĚF?
\ No newline at end of file

models/dolicorhynchops_10k.py

@@ -12,8 +12,8 @@ def parms(d={}):
         { 'file': f'{path}/Lmuscle.ply', 'force': 100., 'focalpt':f'{path}/LmuscleF.off'},
         { 'file': f'{path}/Rmuscle.off', 'force': 100., 'focalpt':f'{path}/RmuscleF.off'}
     ]
-    d.update(p)
-    return d
+    p.update(d)
+    return p
 
 def getMetafor(p={}):
     import mandiblemodel as model

models/dolicorhynchops_149k.py

+#! /usr/bin/env python3
+# -*- coding: utf-8 -*-
+
+def parms(d={}):
+    p = {}
+    path = 'dolicorhynchops/149k'
+    p['mandible'] = f'{path}/mandible_lfs.ply'
+    p['teeth'] = f'{path}/teeth.ply'
+    p['LTMJ'] = f'{path}/LTMJ.ply'
+    p['RTMJ'] = f'{path}/RTMJ.ply'
+    p['muscles'] = [
+        { 'file': f'{path}/Lmuscle_lfs.ply', 'force': 100., 'focalpt':f'{path}/LmuscleF.off'},
+        { 'file': f'{path}/Rmuscle_lfs.ply', 'force': 100., 'focalpt':f'{path}/RmuscleF.off'}
+    ]
+    p.update(d)
+    return p
+
+def getMetafor(p={}):
+    import mandiblemodel as model
+    return model.getMetafor(parms(p))

models/mandiblemodel.py

@@ -34,14 +34,14 @@ def parms(d={}):
     # numerical parameters
     p['tolNR'] = 1e-6        # [-] equilibrium tolerance
 
-    d.update(p)
-    return d
+    p.update(d)
+    return p
 
 
 def getMetafor(p={}):
 
     p = parms(p)
-    print(f'parameters=\n{p}')
+    print(f'parameters={p}')
 
     metafor = Metafor()
     domain = metafor.getDomain()