Este proyecto resuelve el problema de elasticidad lineal utilizando PETSc en 2d y 3d para coeficientes lamé que son constantes o constantes por celda en una cuadrícula cartesiana utilizando el algoritmo GenEO descrito en https://hal.archives-ouvertes.fr/hal-01170059/document .
Para instalar este paquete, primero necesita una instalación de anaconda. Si no tiene anaconda en su sistema, puede descargar miniconda para Python 3 (https://conda.io/miniconda.html).
Para instalar este proyecto, debes clonarlo.
git clone https://github.com/gouarin/GenEO.git
cd GenEO
A continuación, crearemos un entorno con todos los paquetes necesarios usando el siguiente comando.
conda env create -f environment.yml
Para activar tu entorno
source activate petsc-geneo
Entonces
python setup.py install
En el directorio de este proyecto tienes un directorio demos
con ejemplos 2d y 3d.
Es importante estar en el entorno conda creado anteriormente. Si no es el caso
source activate petsc-geneo
Este es un ejemplo de cómo probar uno de ellos.
mpiexec -np 4 python demo_elasticity_2d.py -AMPCG_verbose -ksp_monitor -PCBNN_verbose
Si la ejecución de demo_elasticity_2d.py
tuvo éxito, debería tener un nombre de archivo 'solution_2d_asm.vts'
.
Para visualizar este archivo, debe instalar paraview (https://www.paraview.org/download/).
paraview
y seleccione archivo->estado de carga.paraview
de este proyecto llamado visu_2d.pvsm
.vts
.Deberías ver los resultados.