Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Including PETSc4py as linear solver #362

Open
wants to merge 29 commits into
base: main
Choose a base branch
from
Open

Including PETSc4py as linear solver #362

wants to merge 29 commits into from

Conversation

enricofacca
Copy link
Collaborator

The linear algebra operations required for solving the optimal transport problem via finite volume are solved via PETSc and its python interfaces.

  • The linalg module has a new class to access the PETSc utilities.
  • Lagrange multiplier for fixing the pressure is removed by passing the kernel of the jacobian to KSP solver.
  • Reduction from the the full jacobian to the reduced system is done via fieldsplit preconditioners.
  • AMG is removed since hypre is available via PETSc and seems way faster.
  • Direct sparse solver based on MUMPS is available.

@enricofacca
including petsc4py in requirement
f3974fe
@enricofacca
formatting according to black
c317a3a
@enricofacca
Merge branch 'main' into dev_facca
8f96ecc
@enricofacca
using only petsc ksp sovlers
f333276
@enricofacca
passing ises and kernel to linear solver
8feb034
@enricofacca
newton works, bregman not
86bfab6
@enricofacca
passing kernel info
c220d18
@enricofacca
able to define a user defined preconditioner for the Schur complement
771780c
@enricofacca
passing contrxt to ksp and pc
5f2e8f4
@enricofacca
cleaing linear algebra and added sinkhorn
7a31345
@enricofacca
test passed
3ec454f
@enricofacca
working in geomloss inclusinon
9699818
@enricofacca
passing field ises
65c17af
@enricofacca
main test for wasserstein
6ae2229
@enricofacca
both dmk and pdhg works but only on squared grids and there is probab…
5142b1c 
…ly some eroor in pdhg with to scaling
@enricofacca
fixed issue for using LU factorization. Fixed memory accumualation ad…
46b2c89 
…ding a method to free it.
@enricofacca
remove amg. kappa for bregmann and newton
c226454
@enricofacca
test with zoom
3f3be72
@enricofacca
orginal unit test
4067855
@enricofacca
restoring old test
bc93ce8
@enricofacca
test for geodetic wasserstein distance
70b5a9a
@enricofacca
removed print from dmk
1971a04
@enricofacca
using proper image refinement
0cf473c
@enricofacca
proper refinement and exact 2d wasserstein distance
cd55e15
@enricofacca
other solvers are removed for the moment
dbb75fa
@enricofacca
before pull request
c3b4cec
@enricofacca
removed kappa for pull requests
8540272
@enricofacca
testing all full and reduced formulation in newton
61d9d96
@enricofacca
removed reference to lagrange multipliers. Add a lift to jacobian mat…
1d1a447 
…rix if a direct solver is used
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Reviewers
No reviews
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

Successfully merging this pull request may close these issues.
1 participant
@enricofacca