Example 3 - deformable walls

[Lina Gundelwein]

Hello,

I want to incorporate deformable walls into my simulation.
After successfully running examples 1 and 2, I am trying to get example 3 of the simulation guide to work.
I followed the simulation guide, the preprocessing with the creation of the displacement.vtp gives following fatal error after the initial iterative parameters and preprocessor and preconditioner parameters:

• ************************************************************
  ****************** f a t a l e r r o r ******************
  ************************************************************
  
  routine cgsisrcgsr
  insufficient floating point workspace
  
  ************************************************************

When running the svsolver, I get the following warning:

• WARNING more data in restart than needed: keeping 1st 5
  
  
  ****** Error occurred in routine <restar >
  Error code :nshg r : 8591
  
  
  Life: 0.00000

and the residuals are all NAN:

• 1 1.000E+00 NaN (****) NaN NaN < 1599- 1|****> [ 0 - 0]
  1 1.000E+00 NaN (****) NaN NaN < 1- 1|****> [ 0 - 0]
  1 2.000E+00 NaN (****) NaN NaN < 1- 1|****> [ 0 - 0]
  1 3.000E+00 NaN (****) NaN NaN < 1- 1|****> [ 0 - 0]
  ...

Same problem occurs if I use the default solver (NS) instead of GMRES.
I am using OS X Yosemite, and I am running everything from command line.

To create the restart file I run

Code: Select all

svpost -all -indir "." -outdir "." -sim_units_cm -start 0 -stop 500 -incr 20 -vtu "all_results" -vtp "all_results.vtp" -ph -laststep

in the directory of example 2, then take the restart.500.0, copy it into the directory of example 3 and rename it to restart.1.0
Is this is where the mistake lies?
Can you help me out?

Best regards,
Lina

[Lina Gundelwein]

I just saw that a new Version is available.
Using that, I still get the fatal error, but the solver runs without errors!

[Hongzhi Lan]

Glad you are using the latest version. Just make download/install both SimVascular and svSolver.

The fatal error is not necessarily fatal. You need to check the displacement.vtp and see if the values are reasonable.

SimVascular Development Team