************************************************************ CalculiX Version 2.16, Copyright(C) 1998-2019 Guido Dhondt CalculiX comes with ABSOLUTELY NO WARRANTY. This is free software, and you are welcome to redistribute it under certain conditions, see gpl.htm ************************************************************ You are using an executable made on Mo 25. Nov 18:56:47 CET 2019 Decascading the MPC's Determining the structure of the matrix: number of equations 1400 number of nonzero lower triangular matrix elements 20932 Starting FSI analysis via preCICE using the geometrically non-linear CalculiX solver... Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. Using up to 1 cpu(s) for the symmetric stiffness/mass contributions. Factoring the system of equations using the symmetric spooles solver Using 1 cpu for spooles. Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. Setting up preCICE participant Calculix, using config file: config.yml ---[precice]  This is preCICE version 2.2.0 ---[precice]  Revision info: no-info [Git failed/Not a repository] ---[precice]  Configuration: Release (Debug and Trace log unavailable) ---[precice]  Configuring preCICE with configuration "./precice-config.xml" ---[precice]  I am participant "Calculix" Set ID Found Read data 'Forces0' found with ID # '2'. Write data 'Displacements0' found with ID # '3'. ---[precice]  Setting up master communication to coupling partner/s ---[precice]  Masters are connected ---[precice]  Setting up preliminary slaves communication to coupling partner/s ---[precice]  Prepare partition for mesh Calculix_Mesh ---[precice]  Gather mesh Calculix_Mesh ---[precice]  Send global mesh Calculix_Mesh ---[precice]  Setting up slaves communication to coupling partner/s ---[precice]  Slaves are connected ---[precice]  iteration: 1 of 50, time-window: 1, time: 0 of 10, time-window-size: 0.001, max-timestep-length: 0.001, ongoing: yes, time-window-complete: no, write-iteration-checkpoint Initializing coupling data Adapter writing coupling data... ---[precice]  initializeData is skipped since no data has to be initialized. Adapter reading coupling data... Reading FORCES coupling data with ID '2'. Adjusting time step for transient step precice_dt dtheta = 0.000010, dtheta = 0.000010, solver_dt = 0.001000 Adapter reading coupling data... Reading FORCES coupling data with ID '2'. Adapter writing checkpoint... increment 1 attempt 1 increment size= 1.000000e-03 sum of previous increments=0.000000e+00 actual step time=1.000000e-03 actual total time=1.000000e-03 iteration 1 Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. Using up to 1 cpu(s) for the symmetric stiffness/mass contributions. Factoring the system of equations using the symmetric spooles solver Using 1 cpu for spooles. Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. average force= 0.017446 time avg. forc= 0.017446 largest residual force= 0.000062 in node 283 and dof 1 largest increment of disp= 5.178901e-06 largest correction to disp= 5.178901e-06 in node 283 and dof 1 no convergence iteration 2 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions. Factoring the system of equations using the symmetric spooles solver Using 1 cpu for spooles. Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. average force= 0.017452 time avg. forc= 0.017452 largest residual force= 0.000000 in node 567 and dof 1 largest increment of disp= 5.178901e-06 largest correction to disp= 1.616256e-09 in node 66 and dof 2 convergence Adapter writing coupling data... Writing DISPLACEMENTS coupling data with ID '3'. Adapter calling advance()... ---[precice]  relative convergence measure: relative two-norm diff = -nan, limit = 0.0001, normalization = -nan, conv = false ---[precice]  relative convergence measure: relative two-norm diff = 1, limit = 0.0001, normalization = 0.953757, conv = false ---[precice]  iteration: 2 of 50, time-window: 1, time: 0 of 10, time-window-size: 0.001, max-timestep-length: 0.001, ongoing: yes, time-window-complete: no, read-iteration-checkpoint Adapter reading checkpoint... Adjusting time step for transient step precice_dt dtheta = 0.000010, dtheta = 0.000010, solver_dt = 0.001000 Adapter reading coupling data... Reading FORCES coupling data with ID '2'. increment 1 attempt 2 increment size= 1.000000e-03 sum of previous increments=0.000000e+00 actual step time=1.000000e-03 actual total time=1.000000e-03 iteration 1 Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. Using up to 1 cpu(s) for the symmetric stiffness/mass contributions. Factoring the system of equations using the symmetric spooles solver Using 1 cpu for spooles. Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. average force= 0.017446 time avg. forc= 0.017446 largest residual force= 0.000062 in node 283 and dof 1 largest increment of disp= 5.178901e-06 largest correction to disp= 5.178901e-06 in node 283 and dof 1 no convergence iteration 2 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions. Factoring the system of equations using the symmetric spooles solver Using 1 cpu for spooles. Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the energy calculation. average force= 0.017452 time avg. forc= 0.017452 largest residual force= 0.000000 in node 567 and dof 1 largest increment of disp= 5.178901e-06 largest correction to disp= 1.616256e-09 in node 66 and dof 2 convergence Adapter writing coupling data... Writing DISPLACEMENTS coupling data with ID '3'. Adapter calling advance()... ---[precice]  relative convergence measure: relative two-norm diff = -nan, limit = 0.0001, normalization = 3.49057e-05, conv = false ---[precice]  relative convergence measure: relative two-norm diff = 0, limit = 0.0001, normalization = 0.953757, conv = true ---[precice] ERROR:  The quasi-Newton update contains NaN values. This means that the quasi-Newton acceleration failed to converge. When writing your own adapter this could indicate that you give wrong information to preCICE, such as identical data in succeeding iterations. Or you do not properly save and reload checkpoints. If you give the correct data this could also mean that the coupled problem is too hard to solve. Try to use a QR filter or increase its threshold (larger epsilon).