search

KratosMultiphysics / Kratos

Kratos Multiphysics (A.K.A Kratos) is a framework for building parallel multi-disciplinary simulation software. Modularity, extensibility and HPC are the main objectives. Kratos has BSD license and is written in C++ with extensive Python interface.
https://kratosmultiphysics.github.io/Kratos/
Other
985 stars 242 forks source link

[Core][IO] wrong output in vtk files when Point3D elements are present in model #12415

Open aronnoordam opened 1 month ago

aronnoordam commented 1 month ago

In case Point3D elements are present in my model, all values on all integration points are inf or Nan.

My guess is that it happens because Point3D has zero integration points and in vtk_output.cpp in WriteIntegrationScalarContainerVariable(, there is a division by the number of integration points.

loumalouomega commented 1 month ago

Well, the extrapolation to nodes from GP is my implementation, I never thought about possibility of 0D elements. I can do a fix.

aronnoordam commented 1 month ago

nice thanks!

aronnoordam commented 1 month ago

I think on line 920 of vtk_ouput.cpp: aux_value /= static_cast<double>(integration_points_number); , there is a division by 0 for point3D elements