Hello Everyone!
First of all i'd like to introduce myself, my name is Jorge and i am working
on a surgical simulation project. I have been looking for an open source FEM
library that i could use for simulating elastic deformations. After an
evaluation process i have decided to go for this lib!
Based on elastostatic.cc test, i have performed a test using a simple cube
surface (triangle elements) mesh (8 points, 12 convexes). I have succesfully
imported the mesh and generated the mesh_fem objects (m_u and m_rhs)
defining its qdim to 3 (i want to simulate the deformation vector) and
"FEM_PK(2,1)". The integration method has been set "IM_TRIANGLE(6)".
Trying to solve the equation with a gravitational force, Dirichlet Boundary
conditions for each point:
for (int i = 0; i < (numberOfDOFs/N); ++i) {
gmm::copy(f, gmm::sub_vector(F, gmm::sub_interval(i*N, N)));
}
// Volumic source term brick.
int numberOfDOFs = this->RHS->nb_dof();
plain_vector F(numberOfDOFs);
this->Model->add_initialized_fem_data("VolumicData", *(this->RHS), F);
getfem::add_source_term_brick(*(this->Model),*(this->IntegrationMethod),
"u", "VolumicData");
// Dirichlet condition
plain_vector F2(numberOfDOFs);
this->Model->add_initialized_fem_data("DirichletData", *(this->RHS), F2);
getfem::add_Dirichlet_condition_with_multipliers(*(this->Model),
*(this->IntegrationMethod), "u", *(this->FEM), 0, "DirichletData");
I get the following error:
Trace 2 in getfem_models.cc, line 1658: Source term assembly for Dirichlet
condition
terminate called after throwing an instance of 'gmm::gmm_error'
what(): Error in getfem_superlu.cc, line 214 void
gmm::SuperLU_solve(const gmm::csc_matrix<T, 0>&, T*, T*, double&, int) [with
T = double]:
SuperLU solve failed: info=25
Aborted
I'd really appreciate any hint on this.
Thanks.
Jorge Ballesteros
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
