Fixed indices in FESeries::Fourier::set_k_vectors().

diff --git a/doc/news/changes/minor/20190108Fehling b/doc/news/changes/minor/20190108Fehling
new file mode 100644 (file)
index 0000000..804f400
--- /dev/null
+++ b/doc/news/changes/minor/20190108Fehling
@@ -0,0 +1,5 @@
+Fixed: Tensor indices of set_k_vectors() in FESeries::Fourier.
+Coefficients will now be distributed in every
+coordinate direction correctly.
+<br>
+(Marc Fehling, 2019/01/08)

diff --git a/source/fe/fe_series_fourier.cc b/source/fe/fe_series_fourier.cc
index 260a221aa28a9a1af454bacd482f6672bd8dea9b..1917836c80e86bb45dd9f4099bb2a5d8a51dff66 100644 (file)
--- a/source/fe/fe_series_fourier.cc
+++ b/source/fe/fe_series_fourier.cc
@@ -54,8 +54,8 @@ namespace
         for (unsigned int k = 0; k < N; ++k)
           {
             k_vectors(i, j, k)[0] = 2. * numbers::PI * i;
-            k_vectors(i, j, k)[0] = 2. * numbers::PI * j;
-            k_vectors(i, j, k)[0] = 2. * numbers::PI * k;
+            k_vectors(i, j, k)[1] = 2. * numbers::PI * j;
+            k_vectors(i, j, k)[2] = 2. * numbers::PI * k;
           }
   }