2019-10-08 01:35:13 +08:00
|
|
|
--- a/src/struct_ls/pfmg3_setup_rap.c Thu Jun 09 15:56:18 2016
|
|
|
|
|
+++ b/src/struct_ls/pfmg3_setup_rap.c Thu Dec 22 17:45:20 2016
|
2017-05-16 17:00:05 +08:00
|
|
|
@@ -19,13 +19,28 @@
|
|
|
|
|
* allow for coarsening to be done in the x- and y-directions also.
|
|
|
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
|
|
|
|
|
|
-#define MapIndex(in_index, cdir, out_index) \
|
|
|
|
|
- hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 2); \
|
|
|
|
|
- cdir = (cdir + 1) % 3; \
|
|
|
|
|
- hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 0); \
|
|
|
|
|
- cdir = (cdir + 1) % 3; \
|
|
|
|
|
- hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 1); \
|
|
|
|
|
- cdir = (cdir + 1) % 3;
|
|
|
|
|
+//#define MapIndex(in_index, cdir, out_index) \
|
|
|
|
|
+// hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 2); \
|
|
|
|
|
+// cdir = (cdir + 1) % 3; \
|
|
|
|
|
+// hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 0); \
|
|
|
|
|
+// cdir = (cdir + 1) % 3; \
|
|
|
|
|
+// hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 1); \
|
|
|
|
|
+// cdir = (cdir + 1) % 3;
|
|
|
|
|
+
|
|
|
|
|
+template<typename T, typename U, typename V>
|
|
|
|
|
+void MapIndex(T& in_index, U& cdir, V& out_index)
|
|
|
|
|
+{
|
|
|
|
|
+ if(cdir > 3) cdir = cdir - (cdir / 3) * 3;
|
|
|
|
|
+ hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 2);
|
|
|
|
|
+ if(cdir == 3) cdir = 1;
|
|
|
|
|
+ else ++cdir;
|
|
|
|
|
+ hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 0);
|
|
|
|
|
+ if(cdir == 3) cdir = 1;
|
|
|
|
|
+ else ++cdir;
|
|
|
|
|
+ hypre_IndexD(out_index, cdir) = hypre_IndexD(in_index, 1);
|
|
|
|
|
+ if(cdir == 3) cdir = 1;
|
|
|
|
|
+ else ++cdir;
|
|
|
|
|
+}
|
|
|
|
|
|
|
|
|
|
/*--------------------------------------------------------------------------
|
|
|
|
|
* Sets up new coarse grid operator stucture.
|
