forked from ignotur/NINA
-
Notifications
You must be signed in to change notification settings - Fork 0
/
without_field_decay.cpp
52 lines (40 loc) · 1.34 KB
/
without_field_decay.cpp
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
#include <cmath>
#include "stars.h"
using namespace std;
double NeutronStar::get_B (double t) {
return B;
}
double NeutronStar::get_incl (double t) {
return i_incl;
}
double NeutronStar::get_P(double t) {
double res, I;
t = t - tau;
//I = 2./5. * M*M_sol*pow(R,2);
I=1e45;
res = P * sqrt(1. + 16*t*3.2e7*pow(1e6, 6)*B*B*pi*pi/(3.*pow(light_velocity,3)*I*P*P));
//res = sqrt(P+2*B*B*t*3.2e7/1.e39);
//cout<<R<<endl;
return res;
}
double NeutronStar::get_dot_P (double t) {
double res, I;
//I = 2./5. * M*M_sol*pow(R,2);
//res = pow(B/3.2e19, 2)/get_P(t);
I=1e45;
res = 8*pi*pi*pow(1e6,6)/3./pow(light_velocity, 3)/I/get_P(t)*B*B;
//res = 1*B*B/1.e39/get_P(t);
return res;
}
parametrs_B::parametrs_B (ifstream * in) {
}
void parametrs_B::print_description (ostream * out) {
*out<<"// Модель без убывания магнитного поля //"<<endl;
*out<<"//----------------------------------------------------------//"<<endl;
}
void parametrs_B::print_parametrs (ostream * out) {
*out<<"// Так как поле не убывает параметров у модели тоже нет //"<<endl;
*out<<"//----------------------------------------------------------//"<<endl;
}
void parametrs_B::print_short (ostream * out) {
}