[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Getfem-users] Segmentation fault (core dumped)

From: Zhenghuai Guo
Subject: [Getfem-users] Segmentation fault (core dumped)
Date: Sun, 5 May 2019 14:26:54 +0000

Hi, Konstantinos and Yves  and Andriy,

I use gmsh to make meshes. In case 1, I generate triangle elements. In case 2, 
I use the recombine function in 
gmsh to generate quadrangle elements - it is actually mixed quadrangle and 
triangle elements, but majority is quadrangle elements. 
This reduces the mesh size.

I run the exact same code. In the code I only add linear_elasticity_brick:
  getfem::add_isotropic_linearized_elasticity_brick(model, mim_1, "u_1", 
"lambda_1", "mu_1");
  getfem::add_isotropic_linearized_elasticity_brick(model, mim_2, "u_2", 
"lambda_2", "mu_2");

and Dirichlet_condition

  getfem::add_Dirichlet_condition_with_multipliers(model, mim_1, "u_1", mf_u_1, 
total_bounday_1,"Dirichlet_value" ); 
  getfem::add_Dirichlet_condition_with_multipliers(model, mim_2, "u_2", mf_u_2, 
total_bounday_2,"Dirichlet_value" );

no source term, no mass term

The code runs ok for meshes of the case 1 and case 2.

But, if I add the contact_brick  as below, only the meshes of the case 1 is ok. 
 ind = getfem::add_integral_large_sliding_contact_brick_raytracing(model, "1", 
release_dist, "0", "1", false); ///
mim_contact_1, total_bounday_1, true, true, "u_1", "pre_mflambda_1"); 
mim_contact_2, total_bounday_2, true, true, "u_2", "pre_mflambda_2"); ///

When I use meshes of the case 2, it gives an error of
Trace 2 in, line 2698: Global generic assembly RHS
Segmentation fault (core dumped)

Could you please advise about this? 

Thank you very much
Zhenghuai Guo

reply via email to

[Prev in Thread] Current Thread [Next in Thread]