Neko  0.8.1
A portable framework for high-order spectral element flow simulations
pnpn_res_fctry.f90
Go to the documentation of this file.
1 ! Copyright (c) 2022-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 neko_config
36  use utils
41  implicit none
42  private
43 
44  public :: pnpn_prs_res_t, pnpn_vel_res_t, &
46 
47 contains
48 
49  subroutine pnpn_prs_res_factory(prs_res)
50  class(pnpn_prs_res_t), allocatable, intent(inout) :: prs_res
51 
52  if (allocated(prs_res)) then
53  deallocate(prs_res)
54  end if
55 
56 
57  if (neko_bcknd_sx .eq. 1) then
58  allocate(pnpn_prs_res_sx_t::prs_res)
59  else if (neko_bcknd_device .eq. 1) then
61  else
62  allocate(pnpn_prs_res_cpu_t::prs_res)
63  end if
64 
65  end subroutine pnpn_prs_res_factory
66 
67  subroutine pnpn_vel_res_factory(vel_res)
68  class(pnpn_vel_res_t), allocatable, intent(inout) :: vel_res
69 
70  if (allocated(vel_res)) then
71  deallocate(vel_res)
72  end if
73 
74  if (neko_bcknd_sx .eq. 1) then
75  allocate(pnpn_vel_res_sx_t::vel_res)
76  else if (neko_bcknd_device .eq. 1) then
78  else
79  allocate(pnpn_vel_res_cpu_t::vel_res)
80  end if
81 
82 
83  end subroutine pnpn_vel_res_factory
84 
85 end module pnpn_res_fctry
Build configurations.
Definition: neko_config.f90:34
integer, parameter neko_bcknd_sx
Definition: neko_config.f90:39
integer, parameter neko_bcknd_device
Definition: neko_config.f90:44
Residuals in the Pn-Pn formulation (CPU version)
Definition: pnpn_res_cpu.f90:2
Defines Pressure residual factory for the Pn-Pn formulation.
subroutine, public pnpn_prs_res_factory(prs_res)
subroutine, public pnpn_vel_res_factory(vel_res)
Residuals in the Pn-Pn formulation (SX version)
Definition: pnpn_res_sx.f90:2
Defines Pressure and velocity residuals in the Pn-Pn formulation.
Definition: pnpn_res.f90:34
Utilities.
Definition: utils.f90:35
Abstract type to compute pressure residual.
Definition: pnpn_res.f90:47
Abstract type to compute velocity residual.
Definition: pnpn_res.f90:53