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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
|
#ifndef ATC_COUPLING_MASS_H
#define ATC_COUPLING_MASS_H
/** owned field/s: MASS_DENSITY */
// ATC headers
#include "ATC_Coupling.h"
// Other headers
#include <map>
#include <string>
namespace ATC {
// Forward declarations
class FE_Engine;
class SpeciesTimeIntegrator;
class ChargeRegulator;
class ConcentrationRegulator;
/**
* @class ATC_CouplingMass
* @brief A class for atom-continuum transfers & control for species transport
*/
//--------------------------------------------------------
//--------------------------------------------------------
// Class ATC_CouplingMass
//--------------------------------------------------------
//--------------------------------------------------------
class ATC_CouplingMass : public ATC_Coupling {
public:
// constructor
ATC_CouplingMass(std::string groupName,
double **& perAtomArray,
LAMMPS_NS::Fix * thisFix,
std::string matParamFile,
ExtrinsicModelType extrinsic = NO_MODEL);
// destructor
virtual ~ATC_CouplingMass();
/** parser/modifier */
virtual bool modify(int narg, char **arg);
/** pre time integration */
virtual void initialize();
/** prior to exchange */
virtual void pre_exchange();
virtual void reset_atoms() { resetNlocal_=true;}
/** compute vector for output */
virtual double compute_vector(int n);
/** output routines */
virtual void output();
protected:
// functions
//---------------------------------------------------------------
/** initialization routines */
//---------------------------------------------------------------
/** constructs all data which is updated with time integration, i.e. fields */
//virtual void construct_time_integration_data();
/** set up data which is dependency managed */
virtual void construct_transfers();
/** sets the position/velocity of the ghost atoms */
virtual void set_ghost_atoms(){};
/** physics specific filter initialization */
void init_filter();
bool resetNlocal_;
// i.e. we only need the correct shape function matrix for restriction
};
};
#endif
|
