]> https://gitweb.dealii.org/ - dealii.git/commitdiff
more sophisticated permutation
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Wed, 30 Apr 2003 12:55:39 +0000 (12:55 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Wed, 30 Apr 2003 12:55:39 +0000 (12:55 +0000)
git-svn-id: https://svn.dealii.org/trunk@7499 0785d39b-7218-0410-832d-ea1e28bc413d

tests/lac/solver.cc

index 5edaf5d21ab7446ea5a4c0fafa2055f5cf60f405..6e831397301786fa132c87bc5e020f4b6baa3ce0 100644 (file)
@@ -104,8 +104,21 @@ int main()
 
       std::vector<unsigned int> permutation(dim);
       std::vector<unsigned int> inverse_permutation(dim);
-      for (unsigned int i=0;i<dim;++i)
-       permutation[i] = dim-i-1;
+
+                                      // Create a permutation: Blocks
+                                      // backwards and every second
+                                      // block backwards
+      unsigned int k = 0;
+      for (unsigned int i=0;i<size-1;++i)
+       for (unsigned int j=0;j<size-1;++j)
+         {
+           if (i % 2)
+             permutation[k++] = (size-i-2) * (size-1) + j;
+           else
+             permutation[k++] = (size-i-2) * (size-1) + size-j-2;
+         }
+      
+
       for (unsigned int i=0;i<dim;++i)
        inverse_permutation[permutation[i]] = i;
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.