This repo was forked from JesseLu's tutorial repo. I've been using it as a sort of testing area for implementation of the matrix exponential with PETSc. The python version is expm.py, and a C version is in the works.
A set of walk-throughs to gradually learn how to use PETSc via the petsc4py package.
- vec_serial.py - work with vectors on a single process.
- mat_serial.py - work with matrices on a single process.
- ksp_serial.py - solve a linear system on a single process.
- da_serial.py - work with distributed arrays on a single process.
- da_maxwell_serial.py - solving maxwell's equations.
- Implement a 2D E&M solver (that will work well with both modes).
- Remove "long way" comments in serial codes.
- Do something with the maxwell 3D code.