FinePreconditioner_complete.F90

Go to the documentation of this file.

! DOUG - Domain decomposition On Unstructured Grids
! Copyright (C) 1998-2006 Faculty of Computer Science, University of Tartu and
! Department of Mathematics, University of Bath
!
! This library is free software; you can redistribute it and/or
! modify it under the terms of the GNU Lesser General Public
! License as published by the Free Software Foundation; either
! version 2.1 of the License, or (at your option) any later version.
!
! This library is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
! Lesser General Public License for more details.
!
! You should have received a copy of the GNU Lesser General Public
! License along with this library; if not, write to the Free Software
! Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
! or contact the authors (University of Tartu, Faculty of Computer Science, Chair
! of Distributed Systems, Liivi 2, 50409 Tartu, Estonia, http://dougdevel.org,
! mailto:info(at)dougdevel.org)

module FinePreconditioner_complete_mod
  use Preconditioner_base_mod
  use subsolvers

  implicit none

contains
  subroutine FinePreconditioner_complete_Init(FP)
    type(FinePreconditioner),intent(inout) :: FP
    
    FP%type = FINE_PRECONDITIONER_TYPE_COMPLETE
    allocate(FP%complete)
    FP%complete%factored = .false.
    FP%complete%nsubsolves = 0
    FP%complete%subsolve_ids => NULL()
  end subroutine FinePreconditioner_complete_Init

  subroutine FinePreconditioner_complete_Apply(FP, sol, rhs)
    type(FinePreconditioner),intent(inout) :: FP
    real(kind=rk),dimension(:),pointer :: sol !< solution
    real(kind=rk),dimension(:),pointer :: rhs !< right hand side
    
    if (.not.FP%complete%factored) then!{
      if (sctls%verbose>4) write(stream,*) "Factorizing 1. level"
      call Factorise_subdomains(FP%domains, FP%distr%A, FP%distr%A_ghost, &
           FP%complete%subsolve_ids)
      FP%complete%factored=.true.
    end if

    ! solve
    if (sctls%verbose>4) write(stream,*) "Solving 1. level"
    call solve_subdomains(sol,FP%domains,FP%complete%subsolve_ids,rhs)

  end subroutine FinePreconditioner_complete_Apply
end module FinePreconditioner_complete_mod