Zachary Ferguson
Zachary Ferguson
I don't think so, #13 was a compilation error if you did not specify the right rounding policy. This error is related to how the trigonometric functions are rounded using...
For transcendental functions, you need to specify the rounding policy. For example: ``` typedef boost::numeric::interval< double, boost::numeric::interval_lib::policies< boost::numeric::interval_lib::save_state< boost::numeric::interval_lib::rounded_transc_std >, boost::numeric::interval_lib::checking_base > > Interval; ``` Read more about the rounding...
Can you post the JSON file you are using? Are you doing something like ``` "time": { "dt": 0.01, "tend": 1, "quasistatic": true } ``` I believe you might need...
I have never tested contact with Spline elements. I know contact works with higher-order Lagrange elements (both tets and hexes). If you could share an example JSON and mesh with...
Copying this discussion from Slack for reference: > Lumping for higher-order elements cannot be done by summing rows (this works only for P1): the resulting lumped matrix may be not...
I have a fix for this I will push it today.
Strange... it is working for me. Can you run it in debug mode? The easiest way to do this is to run the sim in release until it crashes and...
Is this run in debug mode (i.e. run CMake with `-DCMAKE_BUILD_TYPE=Debug`)?
Can you use `gdb` or `lldb` to find where the segfault comes from?
Ah, okay. I am using CHOLMOD instead of MKL Pardiso. I can try to reproduce it on a Linux machine with MKL support. In the meantime try to switch to...