Neko  0.9.0
A portable framework for high-order spectral element flow simulations
pc_identity_device.f90
Go to the documentation of this file.
1 ! Copyright (c) 2021-2023, The Neko Authors
2 ! All rights reserved.
3 !
4 ! Redistribution and use in source and binary forms, with or without
5 ! modification, are permitted provided that the following conditions
6 ! are met:
7 !
8 ! * Redistributions of source code must retain the above copyright
9 ! notice, this list of conditions and the following disclaimer.
10 !
11 ! * Redistributions in binary form must reproduce the above
12 ! copyright notice, this list of conditions and the following
13 ! disclaimer in the documentation and/or other materials provided
14 ! with the distribution.
15 !
16 ! * Neither the name of the authors nor the names of its
17 ! contributors may be used to endorse or promote products derived
18 ! from this software without specific prior written permission.
19 !
20 ! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
21 ! "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
22 ! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
23 ! FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
24 ! COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
25 ! INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
26 ! BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
27 ! LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
28 ! CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
29 ! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
30 ! ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
31 ! POSSIBILITY OF SUCH DAMAGE.
32 !
35  use device
36  use device_math
37  use precon, only : pc_t
38  use num_types, only : rp
39  use, intrinsic :: iso_c_binding, only : c_ptr
40  implicit none
41  private
42 
44  type, public, extends(pc_t) :: device_ident_t
45  contains
46  procedure, pass(this) :: solve => device_ident_solve
47  procedure, pass(this) :: update => device_ident_update
48  end type device_ident_t
49 
50 contains
51 
53  subroutine device_ident_solve(this, z, r, n)
54  integer, intent(in) :: n
55  class(device_ident_t), intent(inout) :: this
56  real(kind=rp), dimension(n), intent(inout) :: z
57  real(kind=rp), dimension(n), intent(inout) :: r
58  type(c_ptr) :: z_d, r_d
59 
60  z_d = device_get_ptr(z)
61  r_d = device_get_ptr(r)
62 
63  call device_copy(z_d, r_d, n)
64 
65  end subroutine device_ident_solve
66 
68  subroutine device_ident_update(this)
69  class(device_ident_t), intent(inout) :: this
70  end subroutine device_ident_update
71 
72 end module device_identity
Return the device pointer for an associated Fortran array.
Definition: device.F90:81
Identity Krylov preconditioner for accelerators.
subroutine device_ident_update(this)
Mandatory update routine (NOP)
subroutine device_ident_solve(this, z, r, n)
The (default) naive preconditioner .
subroutine, public device_copy(a_d, b_d, n)
Copy a vector .
Definition: device_math.F90:76
Device abstraction, common interface for various accelerators.
Definition: device.F90:34
integer, parameter, public rp
Global precision used in computations.
Definition: num_types.f90:12
Krylov preconditioner.
Definition: precon.f90:34
Defines a canonical Krylov preconditioner for accelerators.
Defines a canonical Krylov preconditioner.
Definition: precon.f90:40