Line data Source code
1 : /*
2 : * Copyright (c) 2014-2015: G-CSC, Goethe University Frankfurt
3 : * Author: Andreas Vogel
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 __H__UG__LIB_DISC__SPATIAL_DISC__ROTATING_VELOCITY__
34 : #define __H__UG__LIB_DISC__SPATIAL_DISC__ROTATING_VELOCITY__
35 :
36 : #include "lib_disc/spatial_disc/user_data/std_glob_pos_data.h"
37 :
38 : namespace ug{
39 :
40 : class RotatingVelocity2d
41 : : public StdGlobPosData<RotatingVelocity2d, MathVector<2>, 2>
42 : {
43 : public:
44 0 : RotatingVelocity2d(double _cx, double _cy, double _nu)
45 0 : : cx(_cx), cy(_cy), nu(_nu)
46 0 : {};
47 :
48 : inline void evaluate(MathVector<2>& val, const MathVector<2>& pos, number time, int si) const
49 : {
50 0 : const number x = pos[0];
51 0 : const number y = pos[1];
52 :
53 0 : val[0] = nu*(y - cx);
54 0 : val[1] = nu*(cy - x);
55 : }
56 :
57 : protected:
58 : double cx, cy, nu;
59 : };
60 :
61 : }
62 : #endif
|