Skip to content
Snippets Groups Projects
Commit 45b1006a authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

remove unused data

parent 1334fe2c
No related branches found
No related tags found
No related merge requests found
......@@ -143,50 +143,28 @@ fluid_solver_base<R>::fluid_solver_base( \
this->kx = new double[this->cd->sizes[2]]; \
this->ky = new double[this->cd->subsizes[0]]; \
this->kz = new double[this->cd->sizes[1]]; \
this->knullx = new bool[this->cd->sizes[2]]; \
this->knully = new bool[this->cd->subsizes[0]]; \
this->knullz = new bool[this->cd->sizes[1]]; \
this->nonzerokx = int(this->rd->sizes[2] / 3); \
this->kMx = this->dkx*(this->nonzerokx-1); \
this->nonzeroky = int(this->rd->sizes[1] / 3); \
this->kMy = this->dky*(this->nonzeroky-1); \
this->nonzeroky = 2*this->nonzeroky - 1; \
this->nonzerokz = int(this->rd->sizes[0] / 3); \
this->kMz = this->dkz*(this->nonzerokz-1); \
this->nonzerokz = 2*this->nonzerokz - 1; \
this->kMx = this->dkx*(int(this->rd->sizes[2] / 3)-1); \
this->kMy = this->dky*(int(this->rd->sizes[1] / 3)-1); \
this->kMz = this->dkz*(int(this->rd->sizes[0] / 3)-1); \
int i, ii; \
for (i = 0; i<this->cd->sizes[2]; i++) \
{ \
this->kx[i] = i*this->dkx; \
if (i < this->nonzerokx) \
this->knullx[i] = false; \
else \
this->knullx[i] = true; \
} \
for (i = 0; i<this->cd->subsizes[0]; i++) \
{ \
int tval = (this->nonzeroky+1)/2; \
ii = i + this->cd->starts[0]; \
if (ii <= this->rd->sizes[1]/2) \
this->ky[i] = this->dky*ii; \
else \
this->ky[i] = this->dky*(ii - this->rd->sizes[1]); \
if (ii < tval || (this->rd->sizes[1] - ii) < tval) \
this->knully[i] = false; \
else \
this->knully[i] = true; \
} \
for (i = 0; i<this->cd->sizes[1]; i++) \
{ \
int tval = (this->nonzerokz+1)/2; \
if (i <= this->rd->sizes[0]/2) \
this->kz[i] = this->dkz*i; \
else \
this->kz[i] = this->dkz*(i - this->rd->sizes[0]); \
if (i < tval || (this->rd->sizes[0] - i) < tval) \
this->knullz[i] = false; \
else \
this->knullz[i] = true; \
} \
this->kM = this->kMx; \
if (this->kM < this->kMy) this->kM = this->kMy; \
......@@ -257,9 +235,6 @@ fluid_solver_base<R>::~fluid_solver_base() \
delete[] this->kx;\
delete[] this->ky;\
delete[] this->kz;\
delete[] this->knullx;\
delete[] this->knully;\
delete[] this->knullz;\
\
delete this->cd; \
delete this->rd; \
......
......@@ -55,8 +55,8 @@ class fluid_solver_base
/* mode and dealiasing information */
double kMx, kMy, kMz, kM, kM2;
double *kx, *ky, *kz;
bool *knullx, *knully, *knullz;
int nonzerokx, nonzeroky, nonzerokz;
//double *Fourier_filter;
//int Fourier_filter_size;
double *kshell;
int64_t *nshell;
int nshells;
......@@ -74,6 +74,7 @@ class fluid_solver_base
~fluid_solver_base();
void low_pass_Fourier(cnumber *a, int howmany, double kmax);
//void dealias(cnumber *a, int howmany, double kmax);
void force_divfree(cnumber *a);
void symmetrize(cnumber *a, int howmany);
void clean_up_real_space(rnumber *a, int howmany);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment