fork download

copy

#include <iostream>
#include <cmath>
using namespace std;
 
const int Nx = 199;
const int Ny = 199;
const int Nz = 199;
double kx_e[Nx];
double ky_e[Ny];
double kz_e[Nz];
 
double E_p(double alpha, double beta, double mu, double kx, double ky, double kz){
return (alpha + 9/4*beta)*(pow(kx,2)+pow(ky,2)+pow(kz,2))-mu;
}
 
double E_m(double alpha, double beta, double mu, double kx, double ky, double kz){
return (alpha + 1/4*beta)*(pow(kx,2)+pow(ky,2)+pow(kz,2))-mu;
}
 
double c_xz(double beta, double kx, double kz){
return sqrt(3)*beta*kx*kz;
}
 
double c_yz(double beta, double ky, double kz){
return sqrt(3)*beta*ky*kz;
}
 
double Pfaff(double kx, double ky, double kz, double Delta_0){
return pow(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu),2)+4*pow(Delta_0,2)*(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu)+pow(c_xz(kx,kz,beta),2)+pow(c_yz(ky,kz,beta),2));
}
 
void write_data(double x, double y, double z)
{
cout << x << '\t' << y << '\t' << z << endl;
}
 
 
int main(){
double ke_xlimits[2] = {-a, a};
double ke_ylimits[2] = {-a, a};
double ke_zlimits[2] = {-a, a};
 
kx_e[0] = ke_xlimits[0];
ky_e[0] = ke_ylimits[0];
kz_e[0] = ke_zlimits[0];
 
 
for(int i = 0; i < Nx; i++){
kx_e[i] = kx_e[0] + i*(ke_xlimits[1]-ke_xlimits[0])/(Nx-1);
}
for(int j = 0; j < Ny; j++){
ky_e[j] = ky_e[0] + j*(ke_ylimits[1]-ke_ylimits[0])/(Ny-1);
}
for(int l = 1; l < Nz; l++){
kz_e[l] = kz_e[0] + (ke_zlimits[1]-ke_zlimits[0])/(Nz-1);
}
if(Pfaff(kx_e[i], ky_e[j], kz_e[l], Delta_0) == 0) {
double x = kx_e[i];
double y = ky_e[j];
double z = kz_e[l];
write_data(x, y, z);
continue;
}
}
}
}
return 0;
} 
 

prog.cpp: In function ‘double Pfaff(double, double, double, double)’:
prog.cpp:29:28: error: ‘alpha’ was not declared in this scope
 return pow(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu),2)+4*pow(Delta_0,2)*(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu)+pow(c_xz(kx,kz,beta),2)+pow(c_yz(ky,kz,beta),2));
                            ^~~~~
prog.cpp:29:28: note: suggested alternative: ‘isalpha’
 return pow(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu),2)+4*pow(Delta_0,2)*(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu)+pow(c_xz(kx,kz,beta),2)+pow(c_yz(ky,kz,beta),2));
                            ^~~~~
                            isalpha
prog.cpp:29:35: error: ‘beta’ was not declared in this scope
 return pow(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu),2)+4*pow(Delta_0,2)*(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu)+pow(c_xz(kx,kz,beta),2)+pow(c_yz(ky,kz,beta),2));
                                   ^~~~
prog.cpp:29:35: note: suggested alternative: ‘getw’
 return pow(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu),2)+4*pow(Delta_0,2)*(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu)+pow(c_xz(kx,kz,beta),2)+pow(c_yz(ky,kz,beta),2));
                                   ^~~~
                                   getw
prog.cpp:29:41: error: ‘mu’ was not declared in this scope
 return pow(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu),2)+4*pow(Delta_0,2)*(E_p(kx, ky, kz, alpha, beta, mu)*E_m(kx, ky, kz, alpha, beta, mu)+pow(c_xz(kx,kz,beta),2)+pow(c_yz(ky,kz,beta),2));
                                         ^~
prog.cpp: In function ‘int main()’:
prog.cpp:39:26: error: ‘a’ was not declared in this scope
 double ke_xlimits[2] = {-a, a};
                          ^
prog.cpp:57:15: error: ‘i’ was not declared in this scope
 if(Pfaff(kx_e[i], ky_e[j], kz_e[l], Delta_0) == 0) {
               ^
prog.cpp:57:24: error: ‘j’ was not declared in this scope
 if(Pfaff(kx_e[i], ky_e[j], kz_e[l], Delta_0) == 0) {
                        ^
prog.cpp:57:24: note: suggested alternative: ‘jn’
 if(Pfaff(kx_e[i], ky_e[j], kz_e[l], Delta_0) == 0) {
                        ^
                        jn
prog.cpp:57:33: error: ‘l’ was not declared in this scope
 if(Pfaff(kx_e[i], ky_e[j], kz_e[l], Delta_0) == 0) {
                                 ^
prog.cpp:57:37: error: ‘Delta_0’ was not declared in this scope
 if(Pfaff(kx_e[i], ky_e[j], kz_e[l], Delta_0) == 0) {
                                     ^~~~~~~
prog.cpp:62:1: error: continue statement not within a loop
 continue;
 ^~~~~~~~
prog.cpp: At global scope:
prog.cpp:65:1: error: expected declaration before ‘}’ token
 }
 ^

Standard output is empty