Loading [MathJax]/extensions/tex2jax.js
ug4
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules Pages
ivector_vec_functions.h
Go to the documentation of this file.
1/*
2 * Copyright (c) 2013-2015: G-CSC, Goethe University Frankfurt
3 * Author: Martin Rupp
4 *
5 * This file is part of UG4.
6 *
7 * UG4 is free software: you can redistribute it and/or modify it under the
8 * terms of the GNU Lesser General Public License version 3 (as published by the
9 * Free Software Foundation) with the following additional attribution
10 * requirements (according to LGPL/GPL v3 §7):
11 *
12 * (1) The following notice must be displayed in the Appropriate Legal Notices
13 * of covered and combined works: "Based on UG4 (www.ug4.org/license)".
14 *
15 * (2) The following notice must be displayed at a prominent place in the
16 * terminal output of covered works: "Based on UG4 (www.ug4.org/license)".
17 *
18 * (3) The following bibliography is recommended for citation and must be
19 * preserved in all covered files:
20 * "Reiter, S., Vogel, A., Heppner, I., Rupp, M., and Wittum, G. A massively
21 * parallel geometric multigrid solver on hierarchically distributed grids.
22 * Computing and visualization in science 16, 4 (2013), 151-164"
23 * "Vogel, A., Reiter, S., Rupp, M., Nägel, A., and Wittum, G. UG4 -- a novel
24 * flexible software system for simulating pde based models on high performance
25 * computers. Computing and visualization in science 16, 4 (2013), 165-179"
26 *
27 * This program is distributed in the hope that it will be useful,
28 * but WITHOUT ANY WARRANTY; without even the implied warranty of
29 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
30 * GNU Lesser General Public License for more details.
31 */
32
33#ifndef IVECTOR_VEC_FUNCTIONS_H_
34#define IVECTOR_VEC_FUNCTIONS_H_
35
36
38// T, IVector functions
39template<typename T>
40inline void VecAdd(double a1, T &v1, double a2, const IVector &v2)
41{
42 VecAdd(a1, v1, a2, v2.downcast<T>());
43}
44
45template<typename T>
46inline void VecAdd(double a1, T &v1, double a2, const IVector &v2, double a3, const IVector &v3)
47{
48 VecAdd(a1, v1, a2, v2.downcast<T>(), a3, v3.downcast<T>());
49}
50
51template<typename T>
52inline void VecAdd(double a1, T &v1, double a2, const IVector &v2, double a3, const IVector &v3, double a4, const IVector &v4)
53{
54 VecAdd(a1, v1, a2, v2.downcast<T>(), a3, v3.downcast<T>(), a4, v4.downcast<T>());
55}
56
57template<typename T>
58inline double VecProd(T &v1, IVector &v2)
59{
60 return VecProd(v1, v2.downcast<T>());
61}
62
64// IVector, IVector functions
65void VecAdd(double a1, IVector &v1, double a2, const IVector &v2)
66{
67 v1.vec_add(a1, a2, v2);
68}
69
70void VecAdd(double a1, IVector &v1, double a2, const IVector &v2, double a3, const IVector &v3)
71{
72 v1.vec_add(a1, a2, v2, a3, v3);
73}
74
75void VecAdd(double a1, IVector &v1, double a2, const IVector &v2, double a3, const IVector &v3, double a4, const IVector &v4)
76{
77 v1.vec_add(a1, a2, v2, a3, v3, a4, v4);
78}
79
80double VecProd(const IVector &v1, const IVector &v2)
81{
82 return v1.vec_prod(v2);
83}
84
85double VecNorm2(const IVector &v1)
86{
87 return v1.norm2();
88}
89
90double VecNorm(const IVector &v1)
91{
92 return v1.norm();
93}
94
95
96#endif /* IVECTOR_VEC_FUNCTIONS_H_ */
double VecNorm(const IVector &v1)
Definition ivector_vec_functions.h:90
double VecProd(T &v1, IVector &v2)
Definition ivector_vec_functions.h:58
double VecNorm2(const IVector &v1)
Definition ivector_vec_functions.h:85
void VecAdd(double a1, T &v1, double a2, const IVector &v2)
Definition ivector_vec_functions.h:40