enable Double64 for mpr, change nepsType and tol_rel for mpr

Author: AndreaNeumayr

src/Basics/constantsAndFunctions.jl

@@ -9,7 +9,7 @@
 # Epsilon and sign
 const neps = sqrt( eps() )
 
-nepsType(::Type{T}) where {T} = sqrt( eps(T) ) # mpr and bounding box calculation use this
+nepsType(::Type{T}) where {T} = 100.0 * eps(T) # mpr and bounding box calculation use this
 
 function sign_eps(value::T) where {T}
     seps::T = 100.0*nepsType(T)

src/Composition/scene.jl

@@ -421,7 +421,7 @@ mutable struct Scene <: Modia3D.AbstractScene
     function Scene(;gravityField          = UniformGravityField(),
             useOptimizedStructure         = true,
             enableContactDetection        = true,
-            mpr_tol_rel                   = 1.0e-7,
+            mpr_tol_rel                   = 1.0e-20,
             elasticContactReductionFactor = 1.0,
             gap                           = 0.001,
             enableVisualization           = true,

src/Modia3D.jl

@@ -96,8 +96,8 @@ convertAndStripUnit(TargetType, requiredUnit, value) =
             convert(TargetType, ustrip.( uconvert.(requiredUnit, value))) : convert(TargetType, value)
 
 # MPRFloatType is used to change betweeen Double64 and Float64 for mpr calculations
-# using DoubleFloats
-const MPRFloatType = Float64
+using DoubleFloats
+const MPRFloatType = Double64
 
 # Include sub-modules
 include(joinpath("Frames"           , "_module.jl"))