-
Notifications
You must be signed in to change notification settings - Fork 6
/
systemmodel.cpp
53 lines (43 loc) · 1.36 KB
/
systemmodel.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
/*
* Copyright (c) 2012, Robert Rueger <[email protected]>
*
* This file is part of SSMC.
*
* SSMC is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* SSMC is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with SSMC. If not, see <http://www.gnu.org/licenses/>.
*/
#include "systemmodel.hpp"
// ----- SYSTEM MODEL TEMPLATE (DEFINES INTERFACE) -----
SystemModel::SystemModel( const bool& periodic, const double& B,
const double& T, const string& cwd )
: periodic( periodic ), B( B ), T( T ), time( 0 ), cwd( cwd )
{
// initialize the model's random number generator (Mersenne Twister)
rng = gsl_rng_alloc( gsl_rng_mt19937 );
gsl_rng_set( rng, rand() );
}
SystemModel::~SystemModel()
{
// free rng's memory
gsl_rng_free( rng );
}
void SystemModel::set_T( const double& newT )
{
if ( newT >= 0 ) {
T = newT;
}
}
void SystemModel::set_B( const double& newB )
{
B = newB;
}