SLES snes |
MAT A |
Vec x,F |
integer n,its,ierr |
call MatCreat(MPI_COMM_WORLD,n,n,J,ierr) |
call VecCreate(MPI_COMM_WORLD,n,x,ierr) |
call VecDuplicate(x,F,ierr) |
call SNESCreate(MPI_COMM_WORLD,SNES_NONLINEAR_EQUATIONS,snes,ierr) |
call SNESSetFunction(snes,F,EvaluateFunction,PETSC_NULL,ierr) |
call SNESSetJacobian(snes,J,EvaluateJacobian,PETSC_NULL,ierr) |
call SNESSetFromOptions(sles,ierr) |
call SNESSolve(sles,b,x,its,ierr) |
call SNESDestroy(snes,ierr) |