Alsvinn  0.5.3
The fast FVM simulator with UQ support
tecno_utils.hpp
Go to the documentation of this file.
1 /* Copyright (c) 2018 ETH Zurich, Kjetil Olsen Lye
2  * This program is free software: you can redistribute it and/or modify
3  * it under the terms of the GNU General Public License as published by
4  * the Free Software Foundation, either version 3 of the License, or
5  * (at your option) any later version.
6  *
7  * This program is distributed in the hope that it will be useful,
8  * but WITHOUT ANY WARRANTY; without even the implied warranty of
9  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
10  * GNU General Public License for more details.
11  *
12  * You should have received a copy of the GNU General Public License
13  * along with this program. If not, see <http://www.gnu.org/licenses/>.
14  */
15 
16 #pragma once
17 #include "alsfvm/types.hpp"
19 
20 namespace alsfvm {
21 namespace numflux {
22 inline __device__ __host__ real bar(real left, real right) {
23  return (left + right) / 2;
24 }
25 
26 inline __device__ __host__ real diff(real left, real right) {
27  return right - left;
28 }
29 inline __device__ __host__ real ln(real left, real right) {
30  if (diff(left, right) == 0 || diff(log(left), log(right)) == 0) {
31  return 0;
32  }
33 
34  return diff(left, right) / diff(log(left), log(right));
35 }
36 
39 inline __device__ __host__ real divLn(real left, real right, real F) {
40  auto div = ln(left, right);
41 
42  if (div == 0) {
43  return 0;
44  } else {
45  return F / div;
46  }
47 }
48 
50  if (a == 0) {
51  return 0;
52  } else {
53  return a / b;
54  }
55 }
56 }
57 }
__device__ __host__ real div(real a, real b)
Definition: tecno_utils.hpp:49
#define __host__
Definition: types.hpp:46
double real
Definition: types.hpp:65
__device__ __host__ real diff(real left, real right)
Definition: tecno_utils.hpp:26
__device__ __host__ real bar(real left, real right)
Definition: tecno_utils.hpp:22
#define __device__
Definition: types.hpp:45
Various utility functions to implement the tecno flux.
Definition: types.hpp:30
__device__ __host__ real divLn(real left, real right, real F)
Definition: tecno_utils.hpp:39
__device__ __host__ real ln(real left, real right)
Definition: tecno_utils.hpp:29