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
89
90
91
92
93
94
95
96
97
98
99
100
101
|
#include "ElectronFlux.h"
#include "ATC_Error.h"
#include <iostream>
#include <fstream>
#include <map>
#include <string>
#include <vector>
using ATC_Utility::command_line;
using ATC_Utility::str2dbl;
using std::fstream;
using std::map;
using std::string;
using std::vector;
namespace ATC {
ElectronFlux::ElectronFlux() :
maskX_(false),maskY_(false),maskZ_(false)
{}
ElectronFluxLinear::ElectronFluxLinear(
fstream &fileId, map<string,double> & parameters)
: ElectronFlux(),
electronMobility_(0),
electronDiffusivity_(0)
{
if (!fileId.is_open()) throw ATC_Error("cannot open material file");
vector<string> line;
while(fileId.good()) {
command_line(fileId, line);
if (line.size() == 0) continue;
if (line[0] == "end") return;
double value = str2dbl(line[1]);
if (line[0] == "mobility") {
electronMobility_ = value;
parameters["electron_mobility"] = electronMobility_;
}
else if (line[0] == "diffusivity") {
electronDiffusivity_ = value;
parameters["electron_diffusivity"] = electronDiffusivity_;
}
else if (line[0] == "mask_x") { maskX_ = true; }
else if (line[0] == "mask_y") { maskY_ = true; }
else if (line[0] == "mask_z") { maskZ_ = true; }
else {
throw ATC_Error( "unrecognized material function "+line[0]);
}
}
}
ElectronFluxThermopower::ElectronFluxThermopower(
fstream &fileId, map<string,double> & parameters)
: ElectronFlux(),
electronMobility_(0),
seebeckCoef_(0)
{
if (!fileId.is_open()) throw ATC_Error("cannot open material file");
vector<string> line;
while(fileId.good()) {
command_line(fileId, line);
if (line.size() == 0) continue;
if (line[0] == "end") return;
double value = str2dbl(line[1]);
if (line[0] == "mobility") {
electronMobility_ = value;
parameters["electron_mobility"] = electronMobility_;
}
else if (line[0] == "seebeck") {
seebeckCoef_ = value;
parameters["seebeck_coefficient"] = seebeckCoef_;
}
else {
throw ATC_Error( "unrecognized material function "+line[0]);
}
}
}
ElectronFluxConvection::ElectronFluxConvection(
fstream &fileId, map<string,double> & /* parameters */)
: ElectronFlux()
{
if (!fileId.is_open()) throw ATC_Error("cannot open material file");
vector<string> line;
while(fileId.good()) {
command_line(fileId, line);
if (line.size() == 0) continue;
if (line[0] == "end") return;
else if (line[0] == "mask_x") { maskX_ = true; }
else if (line[0] == "mask_y") { maskY_ = true; }
else if (line[0] == "mask_z") { maskZ_ = true; }
else {
throw ATC_Error( "unrecognized material function "+line[0]);
}
}
}
}
|
