3DSpinWave/main.cpp at master · paavalipopov/3DSpinWave · GitHub

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
#include <iostream>
#include <string>
#include <Eigen>
#include "1DPeriodicCrystal.h"

using namespace Eigen;
using namespace std;

int main(int argc, char* argv[])
{
	int N = 6;
	bool debug = true;
	int kStep = 400;
	int omegaStep = 800;

	double l1 = 10e-4;
	double l2 = 10e-4;
	double d = 1e-4;
	double gamma = 1.76e7;
	double H1 = 350;
	double H2 = 420;
	double M4pi = 1750;

	double omegaH1 = gamma * H1;
	double omegaH2 = gamma * H2;
	double omegaM = gamma * M4pi;

//	double omegaStart = 1.0e10;
//	double omegaStart = omegaH1;
	double omegaStart = omegaH2;
//	double omegaStart = 1;
	double omegaEnd = 1.2e10;
//	double omegaEnd = omegaH2;
//	double omegaEnd = sqrt(omegaH1*(omegaH1 + omegaM));
//	double omegaEnd = sqrt(omegaH2*(omegaH2 + omegaM));
//	double omegaEnd = 2;
//	double ka = 0.4;
	double ka = 1.0;

	double b = 2.0 * M_PI / (l1+l2-2.0*d);

	double kStart= -2.0 * b;
	double kEnd = 2.0 * b;

	omegaEnd = ka * (omegaEnd - omegaStart) + omegaStart;


//	if(argc != 4) {
//		cout << "Not enough arguments" << endl;
//		return -1;
//	}
//SpinWaveProblem1D(int N, int baseSplit, int kStpes, int omegaSteps, double H2, double l1, double l2, double omegaStart, double omegaEnd, double kStart, double kEnd);
//	cin.get();
//	SpinWaveProblem1D A = SpinWaveProblem1D(std::stoi(argv[1]), std::stoi(argv[2]), std::stoi(argv[3]));

	SpinWaveProblem1D D = SpinWaveProblem1D(N, 0, kStep, omegaStep, H1, H2, l1, l2, d, omegaStart, omegaEnd, kStart, kEnd, debug);
}