From 480ce57afa138d0daac4ae750a63ac2ee9e56bd4 Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Thu, 26 May 2022 17:28:48 -0400 Subject: Fix denomination: enstrophy -> alpha --- src/driving.c | 6 +++--- src/main.c | 10 +++++----- src/navier-stokes.c | 4 ++-- src/navier-stokes.h | 4 ++-- 4 files changed, 12 insertions(+), 12 deletions(-) diff --git a/src/driving.c b/src/driving.c index 70c2564..fa26b4b 100644 --- a/src/driving.c +++ b/src/driving.c @@ -12,13 +12,13 @@ int g_test( for(kx=-K1;kx<=K1;kx++){ for (ky=-K2;ky<=K2;ky++){ if(kx==2 && ky==-1){ - g[klookup(kx,ky,K1,K2)]=0.5+sqrt(3)/2*I; + g[klookup(kx,ky,2*K1+1,2*K2+1)]=0.5+sqrt(3)/2*I; } else if(kx==-2 && ky==1){ - g[klookup(kx,ky,K1,K2)]=0.5-sqrt(3)/2*I; + g[klookup(kx,ky,2*K1+1,2*K2+1)]=0.5-sqrt(3)/2*I; } else{ - g[klookup(kx,ky,K1,K2)]=0.; + g[klookup(kx,ky,2*K1+1,2*K2+1)]=0.; } } } diff --git a/src/main.c b/src/main.c index f074720..2fc2fcd 100644 --- a/src/main.c +++ b/src/main.c @@ -38,7 +38,7 @@ _Complex double* set_driving(unsigned int driving, nstrophy_parameters parameter _Complex double* set_init(unsigned int init, nstrophy_parameters parameters); #define COMMAND_UK 1 -#define COMMAND_ENSTROPHY 2 +#define COMMAND_ALPHA 2 #define COMMAND_QUIET 3 #define COMMAND_ENERGY 4 @@ -89,8 +89,8 @@ int main ( else if (command==COMMAND_ENERGY){ energy(parameters.K1, parameters.K2, parameters.N1, parameters.N2, parameters.nsteps, parameters.nu, parameters.delta, parameters.L, u0, g, parameters.print_freq, nthreads); } - else if(command==COMMAND_ENSTROPHY){ - enstrophy(parameters.K1, parameters.K2, parameters.N1, parameters.N2, parameters.nsteps, parameters.nu, parameters.delta, parameters.L, u0, g, parameters.print_freq, nthreads); + else if(command==COMMAND_ALPHA){ + alpha(parameters.K1, parameters.K2, parameters.N1, parameters.N2, parameters.nsteps, parameters.nu, parameters.delta, parameters.L, u0, g, parameters.print_freq, nthreads); } else if(command==COMMAND_QUIET){ quiet(parameters.K1, parameters.K2, parameters.N1, parameters.N2, parameters.nsteps, parameters.nu, parameters.delta, parameters.L, u0, g, nthreads); @@ -206,8 +206,8 @@ int read_args( if(strcmp(argv[i], "uk")==0){ *command=COMMAND_UK; } - else if(strcmp(argv[i], "enstrophy")==0){ - *command=COMMAND_ENSTROPHY; + else if(strcmp(argv[i], "alpha")==0){ + *command=COMMAND_ALPHA; } else if(strcmp(argv[i], "quiet")==0){ *command=COMMAND_QUIET; diff --git a/src/navier-stokes.c b/src/navier-stokes.c index fadc87e..24b4508 100644 --- a/src/navier-stokes.c +++ b/src/navier-stokes.c @@ -120,8 +120,8 @@ int energy( return(0); } -// compute enstrophy as a function of time in the I-NS equation -int enstrophy( +// compute alpha as a function of time in the I-NS equation +int alpha( int K1, int K2, int N1, diff --git a/src/navier-stokes.h b/src/navier-stokes.h index e64eeda..07b7550 100644 --- a/src/navier-stokes.h +++ b/src/navier-stokes.h @@ -18,8 +18,8 @@ int uk( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double d // compute the energy as a function of time int energy( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double delta, double L, _Complex double* u0, _Complex double* g, unsigned int print_freq, unsigned int nthreads); -// compute enstrophy -int enstrophy( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double delta, double L, _Complex double* u0, _Complex double* g, unsigned int print_freq, unsigned int nthreads); +// compute alpha +int alpha( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double delta, double L, _Complex double* u0, _Complex double* g, unsigned int print_freq, unsigned int nthreads); // compute solution as a function of time, but do not print anything (useful for debugging) int quiet( int K1, int K2, int N1, int N2, unsigned int nsteps, double nu, double delta, double L, _Complex double* u0, _Complex double* g, unsigned int nthreads); -- cgit v1.2.3-54-g00ecf