lammps-sys 0.6.0

Generates bindings to LAMMPS' C interface (with optional builds from source)
Documentation 
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

#include "ElectronHeatFlux.h"
#include "ATC_Error.h"

#include <iostream>
#include <fstream>
#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 {

ElectronHeatFlux::ElectronHeatFlux(ElectronHeatCapacity * electronHeatCapacity)
  :
  electronHeatCapacity_(electronHeatCapacity)
{
  // do nothing
}

ElectronHeatFluxLinear::ElectronHeatFluxLinear(fstream &fileId, map<string,double> & parameters,
                                               ElectronHeatCapacity * electronHeatCapacity) 
  : ElectronHeatFlux(electronHeatCapacity),
  conductivity_(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;
    else if (line[0] == "conductivity") {
      conductivity_ = str2dbl(line[1]);
      parameters["electron_thermal_conductivity"] = conductivity_;
    }
    else {
      throw ATC_Error( "unrecognized material function "+line[0]);
    }
  }
}

ElectronHeatFluxPowerLaw::ElectronHeatFluxPowerLaw(fstream &fileId, map<string,double> & parameters,
                                                   ElectronHeatCapacity * electronHeatCapacity) 
  : ElectronHeatFlux(electronHeatCapacity),
  conductivity_(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;
    else if (line[0] == "conductivity") {
      conductivity_ = str2dbl(line[1]);
      parameters["electron_thermal_conductivity"] = conductivity_;
    }
    else {
      throw ATC_Error( "unrecognized material function "+line[0]);
    }
  }
}

ElectronHeatFluxThermopower::ElectronHeatFluxThermopower(
  fstream &fileId, map<string,double> & parameters,
  /*const*/ ElectronFlux * electronFlux,
  ElectronHeatCapacity * electronHeatCapacity) 
  : ElectronHeatFlux(electronHeatCapacity),
    conductivity_(0),
    seebeckCoef_(0),
    electronFlux_(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;
    double value = str2dbl(line[1]);
    if (line[0] == "conductivity") {
      conductivity_ = value;
      parameters["electron_thermal_conductivity"] = conductivity_;
    }
    else {
      throw ATC_Error( "unrecognized material function "+line[0]);
    }
    
    seebeckCoef_ = parameters["seebeck_coefficient"];
  }
}

}