https://github.com/hpc-maths/GenEO
Revision 907bdfd332a7ec2c33ea11201cf029397cc15eab authored by Nicole Spillane on 20 May 2021, 09:54:10 UTC, committed by Nicole Spillane on 20 May 2021, 09:54:10 UTC
1 parent 7df0dea
Tip revision: 907bdfd332a7ec2c33ea11201cf029397cc15eab authored by Nicole Spillane on 20 May 2021, 09:54:10 UTC
the user can now specify directly the values of tau for the GenEO eigenvalue problems. A separate file for testing Apos has been created. The option PCNewH2addCS was removed as it is not a good choice for the H2 that is incorporated into H3.
the user can now specify directly the values of tau for the GenEO eigenvalue problems. A separate file for testing Apos has been created. The option PCNewH2addCS was removed as it is not a good choice for the H2 that is incorporated into H3.
Tip revision: 907bdfd
readme.md
[![Binder](https://mybinder.org/badge.svg)](https://mybinder.org/v2/gh/gouarin/GenEO/master)
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.
Computing file changes ...