removing unnnecessary .

Author: leios

contents/quantum_systems/code/c++/energy.cpp

@@ -5,23 +5,23 @@
 
 #include <fftw3.h>
 
-void fft(std::vector<std::complex<double>> x, int n, bool inverse) {
+void fft(std::vector<std::complex<double>> x, bool inverse) {
  std::vector<std::complex<double>> y(x.size(), std::complex<double>(0.0, 0.0));
 
  fftw_plan p;
 
  fftw_complex *in = reinterpret_cast<fftw_complex*>(x.data());
  fftw_complex *out = reinterpret_cast<fftw_complex*>(y.data());
 
- p = fftw_plan_dft_1d(n, in, out,
+ p = fftw_plan_dft_1d(x.size(), in, out,
  (inverse ? FFTW_BACKWARD : FFTW_FORWARD), FFTW_ESTIMATE);
 
 
  fftw_execute(p);
  fftw_destroy_plan(p);
 
  for (size_t i = 0; i < x.size(); ++i) {
- x[i] = y[i] / sqrt(static_cast<double>(n));
+ x[i] = y[i] / sqrt(static_cast<double>(x.size()));
  }
 }
 
@@ -31,7 +31,7 @@ double calculate_energy(std::vector<std::complex<double>> wfc,
  double dx, size_t size) {
  std::vector<std::complex<double>> wfc_k(wfc);
  std::vector<std::complex<double>> wfc_c(size);
- fft(wfc_k, size, false);
+ fft(wfc_k, false);
 
  for (size_t i = 0; i < size; ++i) {
  wfc_c[i] = conj(wfc[i]);
@@ -44,7 +44,7 @@ double calculate_energy(std::vector<std::complex<double>> wfc,
  energy_k[i] = wfc_k[i] * h_k[i];
  }
 
- fft(energy_k, size, true);
+ fft(energy_k, true);
 
  for (size_t i = 0; i < size; ++i) {
  energy_k[i] *= wfc_c[i];

contents/split-operator_method/code/c++/split_op.cpp

@@ -79,20 +79,20 @@ struct Operators {
  vector_complex wfc;
 };
 
-void fft(vector_complex &x, int n, bool inverse) {
+void fft(vector_complex &x, bool inverse) {
  std::vector<std::complex<double>> y(x.size(), std::complex<double>(0.0, 0.0));
  fftw_plan p;
 
  fftw_complex *in = reinterpret_cast<fftw_complex*>(x.data());
  fftw_complex *out = reinterpret_cast<fftw_complex*>(y.data());
- p = fftw_plan_dft_1d(n, in, out,
+ p = fftw_plan_dft_1d(x.size(), in, out,
  (inverse ? FFTW_BACKWARD : FFTW_FORWARD), FFTW_ESTIMATE);
 
  fftw_execute(p);
  fftw_destroy_plan(p);
 
- for (size_t i = 0; i < n; ++i) {
- x[i] = y[i] / sqrt(static_cast<double>(n));
+ for (size_t i = 0; i < x.size(); ++i) {
+ x[i] = y[i] / sqrt(static_cast<double>(x.size()));
  }
 }
 
@@ -104,13 +104,13 @@ void split_op(Params &par, Operators &opr) {
  opr.wfc[j] *= opr.pe[j];
  }
 
- fft(opr.wfc, opr.size, false);
+ fft(opr.wfc, false);
 
  for (size_t j = 0; j < opr.size; ++j) {
  opr.wfc[j] *= opr.ke[j];
  }
 
- fft(opr.wfc, opr.size, true);
+ fft(opr.wfc, true);
 
  for (size_t j = 0; j < opr.size; ++j) {
  opr.wfc[j] *= opr.pe[j];
@@ -159,7 +159,7 @@ double calculate_energy(Params &par, Operators &opr) {
  vector_complex wfc_r(opr.wfc);
  vector_complex wfc_k(opr.wfc);
  vector_complex wfc_c(opr.size);
- fft(wfc_k, opr.size, false);
+ fft(wfc_k, false);
 
  for (size_t i = 0; i < opr.size; ++i) {
  wfc_c[i] = conj(wfc_r[i]);
@@ -172,7 +172,7 @@ double calculate_energy(Params &par, Operators &opr) {
  energy_k[i] = wfc_k[i] * pow(complex(par.k[i], 0.0), 2);
  }
 
- fft(energy_k, opr.size, true);
+ fft(energy_k, true);
 
  for (size_t i = 0; i < opr.size; ++i) {
  energy_k[i] *= 0.5 * wfc_c[i];