2022-05-18 07:57:06 +00:00
|
|
|
#include "driving.h"
|
2022-05-26 19:16:44 +00:00
|
|
|
#include "navier-stokes.h"
|
2022-05-18 07:57:06 +00:00
|
|
|
#include <math.h>
|
|
|
|
|
2022-05-26 19:16:44 +00:00
|
|
|
// test driving function
|
|
|
|
int g_test(
|
|
|
|
_Complex double* g,
|
|
|
|
int K1,
|
|
|
|
int K2
|
2022-05-18 07:57:06 +00:00
|
|
|
){
|
2022-05-26 19:16:44 +00:00
|
|
|
int kx,ky;
|
|
|
|
for(kx=-K1;kx<=K1;kx++){
|
|
|
|
for (ky=-K2;ky<=K2;ky++){
|
|
|
|
if(kx==2 && ky==-1){
|
2022-05-26 21:28:48 +00:00
|
|
|
g[klookup(kx,ky,2*K1+1,2*K2+1)]=0.5+sqrt(3)/2*I;
|
2022-05-26 19:16:44 +00:00
|
|
|
}
|
|
|
|
else if(kx==-2 && ky==1){
|
2022-05-26 21:28:48 +00:00
|
|
|
g[klookup(kx,ky,2*K1+1,2*K2+1)]=0.5-sqrt(3)/2*I;
|
2022-05-26 19:16:44 +00:00
|
|
|
}
|
|
|
|
else{
|
2022-05-26 21:28:48 +00:00
|
|
|
g[klookup(kx,ky,2*K1+1,2*K2+1)]=0.;
|
2022-05-26 19:16:44 +00:00
|
|
|
}
|
|
|
|
}
|
2022-05-18 07:57:06 +00:00
|
|
|
}
|
2022-05-26 19:16:44 +00:00
|
|
|
|
|
|
|
return 0;
|
2022-05-18 07:57:06 +00:00
|
|
|
}
|
|
|
|
|
2022-05-26 19:16:44 +00:00
|
|
|
int g_zero(
|
|
|
|
_Complex double* g,
|
|
|
|
int K1,
|
|
|
|
int K2
|
2022-05-19 15:51:45 +00:00
|
|
|
){
|
2022-05-26 19:16:44 +00:00
|
|
|
int i;
|
|
|
|
for(i=0;i<(2*K1+1)*(2*K2+1);i++){
|
|
|
|
g[i]=0.;
|
|
|
|
}
|
2022-05-19 15:51:45 +00:00
|
|
|
|
2022-05-26 19:16:44 +00:00
|
|
|
return 0;
|
|
|
|
}
|