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
|
#ifndef DYNMAT_H
#define DYNMAT_H
#include "stdio.h"
#include "stdlib.h"
#include "string.h"
#include "memory.h"
#include "interpolate.h"
class DynMat {
public:
DynMat(int, char**);
~DynMat();
int nx, ny, nz, nucell;
int sysdim, fftdim;
double eml2f, eml2fc;
char *funit;
void getDMq(double *);
void getDMq(double *, double *);
void writeDMq(double *);
void writeDMq(double *, const double, FILE *fp);
int geteigen(double *, int);
void reset_interp_method();
doublecomplex **DM_q;
int flag_latinfo;
int npt, fftdim2;
double Tmeasure, basevec[9], ibasevec[9];
double *M_inv_sqrt;
double **basis;
int *attyp;
private:
int flag_skip, flag_reset_gamma;
Interpolate *interpolate;
Memory *memory;
int nasr;
void EnforceASR();
char *binfile, *dmfile;
double boltz, q[3];
doublecomplex **DM_all;
void car2dir(); // to convert basis from cartisian coordinate into factional.
void real2rec();
void GaussJordan(int, double *);
void help();
void ShowInfo();
void ShowVersion();
void Define_Conversion_Factor();
};
#endif
|
