This project solves the linear elasticity problem using PETSc in 2d and 3d for lamé coefficients which are constant or constant by cell on a cartesian grid using the GenEO algorithm described in https://hal.archives-ouvertes.fr/hal-01170059/document.

Installation

To install this package, you need first an installation of anaconda. If you don't have anaconda on your system, you can download miniconda for Python 3 (https://conda.io/miniconda.html).

To install this project, you have to clone it

git clone https://github.com/gouarin/GenEO.git
cd GenEO

Next, we will create an environment with all the needed packages using the following command.

conda env create -f environment.yml

To activate your environment

source activate petsc-geneo

Then

python setup.py install

Execute demo file

In the directory of this project you have a demos directory with 2d and 3d examples.

It's important to be in the conda environment created previously. If it is not the case

source activate petsc-geneo

This is an example of how to test one of them

mpiexec -np 4 python demo_elasticity_2d.py -AMPCG_verbose -ksp_monitor -PCBNN_verbose

Visualize the results

If the execution of demo_elasticity_2d.py succeeded, you should have a file name 'solution_2d_asm.vts'.

To visualize this file, you have to install paraview (https://www.paraview.org/download/).

• Start paraview and select file->load state.
• Then select the file in the directory paraview of this project called visu_2d.pvsm.
• Then select the vts file.

You should see the results.