From f1cf33db4a482f0041f46c57fc1e633afe365d3a Mon Sep 17 00:00:00 2001
From: Jochem Snuverink <jochem.snuverink@psi.ch>
Date: Thu, 9 Jan 2020 13:34:53 +0100
Subject: [PATCH] Port IPPL unit tests to OPAL unit tests and cleanup - 2
---
ippl/test/doAllTests | 18 -
ippl/test/simple/ScaV1 | 13 -
ippl/test/simple/VecV1 | 13 -
ippl/test/simple/send | 53 --
ippl/test/simple/stest1 | 12 -
ippl/test/simple/stest2 | 12 -
ippl/test/simple/t1.cpp | 24 -
ippl/test/simple/text.correct.TestBC | 708 -------------------
ippl/test/simple/text.test.TestBC | 708 -------------------
ippl/test/simple/vtest1 | 12 -
ippl/test/simple/vtest2 | 12 -
ippl/test/sindex/compsindex.cpp | 169 -----
ippl/test/sindex/sindex.cpp | 120 ----
ippl/test/sindex/subfield.cpp | 233 ------
ippl/test/sindex/vertcell.cpp | 66 --
tests/ippl_src/Field/BCond.cpp | 9 +-
tests/ippl_src/Field/FieldDebug.cpp | 8 +-
tests/ippl_src/Index/Index.cpp | 346 +++++++++
tests/ippl_src/Meshes/CartesianCentering.cpp | 8 +-
19 files changed, 359 insertions(+), 2185 deletions(-)
delete mode 100755 ippl/test/doAllTests
delete mode 100755 ippl/test/simple/ScaV1
delete mode 100755 ippl/test/simple/VecV1
delete mode 100644 ippl/test/simple/send
delete mode 100755 ippl/test/simple/stest1
delete mode 100755 ippl/test/simple/stest2
delete mode 100644 ippl/test/simple/t1.cpp
delete mode 100755 ippl/test/simple/text.correct.TestBC
delete mode 100755 ippl/test/simple/text.test.TestBC
delete mode 100755 ippl/test/simple/vtest1
delete mode 100755 ippl/test/simple/vtest2
delete mode 100644 ippl/test/sindex/compsindex.cpp
delete mode 100644 ippl/test/sindex/sindex.cpp
delete mode 100644 ippl/test/sindex/subfield.cpp
delete mode 100644 ippl/test/sindex/vertcell.cpp
diff --git a/ippl/test/doAllTests b/ippl/test/doAllTests
deleted file mode 100755
index 3c4dd776d..000000000
--- a/ippl/test/doAllTests
+++ /dev/null
@@ -1,18 +0,0 @@
-#!/bin/bash
-#
-# doAllTests
-#
-#
-#
-#
-
-for arg in FFT parallel simple speed particle sindex
-do
- cd $arg
- echo start test in $arg
- make clean
- make CXX=g++3.2 | tee make.out
- rm *.o
- echo done in echo $arg
- cd ..
-done
diff --git a/ippl/test/simple/ScaV1 b/ippl/test/simple/ScaV1
deleted file mode 100755
index 5982276dc..000000000
--- a/ippl/test/simple/ScaV1
+++ /dev/null
@@ -1,13 +0,0 @@
-vnode = 0
- [-1][-1]=0 [0][-1]=0 [1][-1]=0 [2][-1]=0 [3][-1]=0 [4][-1]=0 [5][-1]=0 [6][-1]=0 [7][-1]=0 [8][-1]=0 [9][-1]=0
- [-1][0]=0 [0][0]=0 [1][0]=0.5 [2][0]=0.5 [3][0]=0.5 [4][0]=0.5 [5][0]=0.5 [6][0]=0.5 [7][0]=0.5 [8][0]=-3.5 [9][0]=0
- [-1][1]=0 [0][1]=0.5 [1][1]=2 [2][1]=2 [3][1]=2 [4][1]=2 [5][1]=2 [6][1]=2 [7][1]=2 [8][1]=-6.5 [9][1]=0
- [-1][2]=0 [0][2]=0.5 [1][2]=2 [2][2]=2 [3][2]=2 [4][2]=2 [5][2]=2 [6][2]=2 [7][2]=2 [8][2]=-6.5 [9][2]=0
- [-1][3]=0 [0][3]=0.5 [1][3]=2 [2][3]=2 [3][3]=2 [4][3]=2 [5][3]=2 [6][3]=2 [7][3]=2 [8][3]=-6.5 [9][3]=0
- [-1][4]=0 [0][4]=0.5 [1][4]=2 [2][4]=2 [3][4]=2 [4][4]=2 [5][4]=2 [6][4]=2 [7][4]=2 [8][4]=-6.5 [9][4]=0
- [-1][5]=0 [0][5]=0.5 [1][5]=2 [2][5]=2 [3][5]=2 [4][5]=2 [5][5]=2 [6][5]=2 [7][5]=2 [8][5]=-6.5 [9][5]=0
- [-1][6]=0 [0][6]=0.5 [1][6]=2 [2][6]=2 [3][6]=2 [4][6]=2 [5][6]=2 [6][6]=2 [7][6]=2 [8][6]=-6.5 [9][6]=0
- [-1][7]=0 [0][7]=0.5 [1][7]=2 [2][7]=2 [3][7]=2 [4][7]=2 [5][7]=2 [6][7]=2 [7][7]=2 [8][7]=-6.5 [9][7]=0
- [-1][8]=0 [0][8]=-3.5 [1][8]=-6.5 [2][8]=-6.5 [3][8]=-6.5 [4][8]=-6.5 [5][8]=-6.5 [6][8]=-6.5 [7][8]=-6.5 [8][8]=-7 [9][8]=0
- [-1][9]=0 [0][9]=0 [1][9]=0 [2][9]=0 [3][9]=0 [4][9]=0 [5][9]=0 [6][9]=0 [7][9]=0 [8][9]=0 [9][9]=0
-------------------
diff --git a/ippl/test/simple/VecV1 b/ippl/test/simple/VecV1
deleted file mode 100755
index 5c06964a3..000000000
--- a/ippl/test/simple/VecV1
+++ /dev/null
@@ -1,13 +0,0 @@
-vnode = 0
- [-1][-1]=( 0 , 0 ) [0][-1]=( 0 , 0 ) [1][-1]=( 0 , 0 ) [2][-1]=( 0 , 0 ) [3][-1]=( 0 , 0 ) [4][-1]=( 0 , 0 ) [5][-1]=( 0 , 0 ) [6][-1]=( 0 , 0 ) [7][-1]=( 0 , 0 ) [8][-1]=( 0 , 0 ) [9][-1]=( 0 , 0 )
- [-1][0]=( 0 , 0 ) [0][0]=( 0 , 0 ) [1][0]=( 0.5 , 0.5 ) [2][0]=( 0.5 , 1.5 ) [3][0]=( 0.5 , 2.5 ) [4][0]=( 0.5 , 3.5 ) [5][0]=( 0.5 , 4.5 ) [6][0]=( 0.5 , 5.5 ) [7][0]=( 0.5 , 6.5 ) [8][0]=( -3.5 , 3.5 ) [9][0]=( 0 , 0 )
- [-1][1]=( 0 , 0 ) [0][1]=( 0 , 0 ) [1][1]=( 1 , 0 ) [2][1]=( 1 , 0 ) [3][1]=( 1 , 0 ) [4][1]=( 1 , 0 ) [5][1]=( 1 , 0 ) [6][1]=( 1 , 0 ) [7][1]=( 1 , 0 ) [8][1]=( -7 , 0 ) [9][1]=( 0 , 0 )
- [-1][2]=( 0 , 0 ) [0][2]=( 0 , 0 ) [1][2]=( 1 , 0 ) [2][2]=( 1 , 0 ) [3][2]=( 1 , 0 ) [4][2]=( 1 , 0 ) [5][2]=( 1 , 0 ) [6][2]=( 1 , 0 ) [7][2]=( 1 , 0 ) [8][2]=( -7 , 0 ) [9][2]=( 0 , 0 )
- [-1][3]=( 0 , 0 ) [0][3]=( 0 , 0 ) [1][3]=( 1 , 0 ) [2][3]=( 1 , 0 ) [3][3]=( 1 , 0 ) [4][3]=( 1 , 0 ) [5][3]=( 1 , 0 ) [6][3]=( 1 , 0 ) [7][3]=( 1 , 0 ) [8][3]=( -7 , 0 ) [9][3]=( 0 , 0 )
- [-1][4]=( 0 , 0 ) [0][4]=( 0 , 0 ) [1][4]=( 1 , 0 ) [2][4]=( 1 , 0 ) [3][4]=( 1 , 0 ) [4][4]=( 1 , 0 ) [5][4]=( 1 , 0 ) [6][4]=( 1 , 0 ) [7][4]=( 1 , 0 ) [8][4]=( -7 , 0 ) [9][4]=( 0 , 0 )
- [-1][5]=( 0 , 0 ) [0][5]=( 0 , 0 ) [1][5]=( 1 , 0 ) [2][5]=( 1 , 0 ) [3][5]=( 1 , 0 ) [4][5]=( 1 , 0 ) [5][5]=( 1 , 0 ) [6][5]=( 1 , 0 ) [7][5]=( 1 , 0 ) [8][5]=( -7 , 0 ) [9][5]=( 0 , 0 )
- [-1][6]=( 0 , 0 ) [0][6]=( 0 , 0 ) [1][6]=( 1 , 0 ) [2][6]=( 1 , 0 ) [3][6]=( 1 , 0 ) [4][6]=( 1 , 0 ) [5][6]=( 1 , 0 ) [6][6]=( 1 , 0 ) [7][6]=( 1 , 0 ) [8][6]=( -7 , 0 ) [9][6]=( 0 , 0 )
- [-1][7]=( 0 , 0 ) [0][7]=( 0 , 0 ) [1][7]=( 1 , 0 ) [2][7]=( 1 , 0 ) [3][7]=( 1 , 0 ) [4][7]=( 1 , 0 ) [5][7]=( 1 , 0 ) [6][7]=( 1 , 0 ) [7][7]=( 1 , 0 ) [8][7]=( -7 , 0 ) [9][7]=( 0 , 0 )
- [-1][8]=( 0 , 0 ) [0][8]=( 0 , 0 ) [1][8]=( 0.5 , -0.5 ) [2][8]=( 0.5 , -1.5 ) [3][8]=( 0.5 , -2.5 ) [4][8]=( 0.5 , -3.5 ) [5][8]=( 0.5 , -4.5 ) [6][8]=( 0.5 , -5.5 ) [7][8]=( 0.5 , -6.5 ) [8][8]=( -3.5 , -3.5 ) [9][8]=( 0 , 0 )
- [-1][9]=( 0 , 0 ) [0][9]=( 0 , 0 ) [1][9]=( 0 , 0 ) [2][9]=( 0 , 0 ) [3][9]=( 0 , 0 ) [4][9]=( 0 , 0 ) [5][9]=( 0 , 0 ) [6][9]=( 0 , 0 ) [7][9]=( 0 , 0 ) [8][9]=( 0 , 0 ) [9][9]=( 0 , 0 )
-------------------
diff --git a/ippl/test/simple/send b/ippl/test/simple/send
deleted file mode 100644
index 7f2ed0973..000000000
--- a/ippl/test/simple/send
+++ /dev/null
@@ -1,53 +0,0 @@
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-TestParticleDebugParallel.cpp: IpplInfo::Comm->barrier();
-bool.cpp: sprintf(name,"out.%d",Ippl::Comm->myNode());
-bool.cpp: Ippl::Comm->barrier();
-conejo_balance.cpp: Ippl::Comm->barrier();
-doof2d.cpp: if( Ippl::Comm->myNode() == Parent ) {
-doof2d.cpp: Ippl::Comm->broadcast_all(mess, tag);
-doof2d.cpp: Message *mess = Ippl::Comm->receive_block(Parent, tag);
-doof2d.cpp: msg << "received a message on node " << Ippl::Comm->myNode();
-doof2d_complex.cpp: if( Ippl::Comm->myNode() == Parent ) {
-doof2d_complex.cpp: Ippl::Comm->broadcast_all(mess, tag);
-doof2d_complex.cpp: Message *mess = Ippl::Comm->receive_block(Parent, tag);
-doof2d_complex.cpp: msg << "received a message on node " << Ippl::Comm->myNode();
-p1.cpp: sprintf(name,"out.%d",Ippl::Comm->myNode());
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p1.cpp: Ippl::Comm->barrier();
-p2.cpp: sprintf(name,"out.%d",Ippl::Comm->myNode());
-p2.cpp: Ippl::Comm->barrier();
-p2.cpp: Ippl::Comm->barrier();
-p2.cpp: Ippl::Comm->barrier();
-p2.cpp: Ippl::Comm->barrier();
-p2.cpp: Ippl::Comm->barrier();
-p2.cpp: Ippl::Comm->barrier();
-p3.cpp: sprintf(name,"out.%d",Ippl::Comm->myNode());
-p3.cpp: Ippl::Comm->barrier();
-randomparallel.cpp: Ippl::Comm->barrier();
diff --git a/ippl/test/simple/stest1 b/ippl/test/simple/stest1
deleted file mode 100755
index 1f2d48885..000000000
--- a/ippl/test/simple/stest1
+++ /dev/null
@@ -1,12 +0,0 @@
-vnode = 0
- [-1][-1]=0 [0][-1]=0 [1][-1]=0 [2][-1]=0 [3][-1]=0 [4][-1]=0 [5][-1]=0 [6][-1]=0 [7][-1]=0 [8][-1]=0
- [-1][0]=0 [0][0]=0 [1][0]=1 [2][0]=2 [3][0]=3 [4][0]=4 [5][0]=5 [6][0]=6 [7][0]=7 [8][0]=0
- [-1][1]=0 [0][1]=0 [1][1]=1 [2][1]=2 [3][1]=3 [4][1]=4 [5][1]=5 [6][1]=6 [7][1]=7 [8][1]=0
- [-1][2]=0 [0][2]=0 [1][2]=1 [2][2]=2 [3][2]=3 [4][2]=4 [5][2]=5 [6][2]=6 [7][2]=7 [8][2]=0
- [-1][3]=0 [0][3]=0 [1][3]=1 [2][3]=2 [3][3]=3 [4][3]=4 [5][3]=5 [6][3]=6 [7][3]=7 [8][3]=0
- [-1][4]=0 [0][4]=0 [1][4]=1 [2][4]=2 [3][4]=3 [4][4]=4 [5][4]=5 [6][4]=6 [7][4]=7 [8][4]=0
- [-1][5]=0 [0][5]=0 [1][5]=1 [2][5]=2 [3][5]=3 [4][5]=4 [5][5]=5 [6][5]=6 [7][5]=7 [8][5]=0
- [-1][6]=0 [0][6]=0 [1][6]=1 [2][6]=2 [3][6]=3 [4][6]=4 [5][6]=5 [6][6]=6 [7][6]=7 [8][6]=0
- [-1][7]=0 [0][7]=0 [1][7]=1 [2][7]=2 [3][7]=3 [4][7]=4 [5][7]=5 [6][7]=6 [7][7]=7 [8][7]=0
- [-1][8]=0 [0][8]=0 [1][8]=0 [2][8]=0 [3][8]=0 [4][8]=0 [5][8]=0 [6][8]=0 [7][8]=0 [8][8]=0
-------------------
diff --git a/ippl/test/simple/stest2 b/ippl/test/simple/stest2
deleted file mode 100755
index 1f2d48885..000000000
--- a/ippl/test/simple/stest2
+++ /dev/null
@@ -1,12 +0,0 @@
-vnode = 0
- [-1][-1]=0 [0][-1]=0 [1][-1]=0 [2][-1]=0 [3][-1]=0 [4][-1]=0 [5][-1]=0 [6][-1]=0 [7][-1]=0 [8][-1]=0
- [-1][0]=0 [0][0]=0 [1][0]=1 [2][0]=2 [3][0]=3 [4][0]=4 [5][0]=5 [6][0]=6 [7][0]=7 [8][0]=0
- [-1][1]=0 [0][1]=0 [1][1]=1 [2][1]=2 [3][1]=3 [4][1]=4 [5][1]=5 [6][1]=6 [7][1]=7 [8][1]=0
- [-1][2]=0 [0][2]=0 [1][2]=1 [2][2]=2 [3][2]=3 [4][2]=4 [5][2]=5 [6][2]=6 [7][2]=7 [8][2]=0
- [-1][3]=0 [0][3]=0 [1][3]=1 [2][3]=2 [3][3]=3 [4][3]=4 [5][3]=5 [6][3]=6 [7][3]=7 [8][3]=0
- [-1][4]=0 [0][4]=0 [1][4]=1 [2][4]=2 [3][4]=3 [4][4]=4 [5][4]=5 [6][4]=6 [7][4]=7 [8][4]=0
- [-1][5]=0 [0][5]=0 [1][5]=1 [2][5]=2 [3][5]=3 [4][5]=4 [5][5]=5 [6][5]=6 [7][5]=7 [8][5]=0
- [-1][6]=0 [0][6]=0 [1][6]=1 [2][6]=2 [3][6]=3 [4][6]=4 [5][6]=5 [6][6]=6 [7][6]=7 [8][6]=0
- [-1][7]=0 [0][7]=0 [1][7]=1 [2][7]=2 [3][7]=3 [4][7]=4 [5][7]=5 [6][7]=6 [7][7]=7 [8][7]=0
- [-1][8]=0 [0][8]=0 [1][8]=0 [2][8]=0 [3][8]=0 [4][8]=0 [5][8]=0 [6][8]=0 [7][8]=0 [8][8]=0
-------------------
diff --git a/ippl/test/simple/t1.cpp b/ippl/test/simple/t1.cpp
deleted file mode 100644
index 21944e6de..000000000
--- a/ippl/test/simple/t1.cpp
+++ /dev/null
@@ -1,24 +0,0 @@
-/***************************************************************************
- *
- * The EXPDE Framework
- *
- *
- ***************************************************************************/
-
-#include "Ippl.h"
-
-int main(int argc, char *argv[])
-{
- Ippl ippl(argc,argv);
- Inform testmsg(argv[0]);
- Inform msg2all(argv[0],INFORM_ALL_NODES);
-
- testmsg << "Test" << endl;
- msg2all << "From all nodes " << endl;
- return 0;
-}
-/***************************************************************************
- * $RCSfile: t1.cpp,v $ $Author: adelmann $
- * $Revision: 1.1.1.1 $ $Date: 2003/03/28 12:47:26 $
- * IPPL_VERSION_ID: $Id: t1.cpp,v 1.1.1.1 2003/03/28 12:47:26 adelmann Exp $
- ***************************************************************************/
diff --git a/ippl/test/simple/text.correct.TestBC b/ippl/test/simple/text.correct.TestBC
deleted file mode 100755
index 2191f15f6..000000000
--- a/ippl/test/simple/text.correct.TestBC
+++ /dev/null
@@ -1,708 +0,0 @@
-++++++++BConds object cbc begin++++++++
-BConds:(
-NegReflectFace, Face=0 ,
-ZeroFace, Face=1 ,
-ParallelPeriodicFace, Face=2 ,
-ParallelPeriodicFace, Face=3
-)
-
-++++++++BConds object cbc end++++++++++
-++++++++++cA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00
-
---------------------------------------------------J = 1
-5.000e+00 6.000e+00 7.000e+00 8.000e+00 9.000e+00
-
---------------------------------------------------J = 2
-1.000e+01 1.100e+01 1.200e+01 1.300e+01 1.400e+01
-
---------------------------------------------------J = 3
-1.500e+01 1.600e+01 1.700e+01 1.800e+01 1.900e+01
-
---------------------------------------------------J = 4
-2.000e+01 2.100e+01 2.200e+01 2.300e+01 2.400e+01
-
-++++++++++cB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
--1.600e+01 -1.500e+01 1.500e+01 1.600e+01 1.700e+01
-
---------------------------------------------------J = 1
--2.100e+01 -2.000e+01 2.000e+01 2.100e+01 2.200e+01
-
---------------------------------------------------J = 2
--1.000e+00 0.000e+00 0.000e+00 1.000e+00 2.000e+00
-
---------------------------------------------------J = 3
--6.000e+00 -5.000e+00 5.000e+00 6.000e+00 7.000e+00
-
---------------------------------------------------J = 4
--1.100e+01 -1.000e+01 1.000e+01 1.100e+01 1.200e+01
-
-++++++++++vA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00
-
---------------------------------------------------J = 1
-5.000e+00 6.000e+00 7.000e+00 8.000e+00 9.000e+00
-
---------------------------------------------------J = 2
-1.000e+01 1.100e+01 1.200e+01 1.300e+01 1.400e+01
-
---------------------------------------------------J = 3
-1.500e+01 1.600e+01 1.700e+01 1.800e+01 1.900e+01
-
---------------------------------------------------J = 4
-0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00
-
-++++++++++vB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
--1.200e+01 -1.100e+01 1.000e+01 1.100e+01 1.200e+01
-
---------------------------------------------------J = 1
--1.700e+01 -1.600e+01 1.500e+01 1.600e+01 1.700e+01
-
---------------------------------------------------J = 2
--2.000e+00 -1.000e+00 0.000e+00 1.000e+00 2.000e+00
-
---------------------------------------------------J = 3
--7.000e+00 -6.000e+00 5.000e+00 6.000e+00 7.000e+00
-
---------------------------------------------------J = 4
--1.200e+01 -1.100e+01 1.000e+01 1.100e+01 1.200e+01
-
-++++++++++vcA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 1.00e+00 )
-( 2.00e+00 , 2.00e+00 ) ( 3.00e+00 , 3.00e+00 )
-( 4.00e+00 , 4.00e+00 )
-
---------------------------------------------------J = 1
-( 5.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 ) ( 8.00e+00 , 8.00e+00 )
-( 9.00e+00 , 9.00e+00 )
-
---------------------------------------------------J = 2
-( 1.00e+01 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 ) ( 1.30e+01 , 1.30e+01 )
-( 1.40e+01 , 1.40e+01 )
-
---------------------------------------------------J = 3
-( 1.50e+01 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 ) ( 1.80e+01 , 1.80e+01 )
-( 1.90e+01 , 1.90e+01 )
-
---------------------------------------------------J = 4
-( 2.00e+01 , 2.00e+01 ) ( 2.10e+01 , 2.10e+01 )
-( 2.20e+01 , 2.20e+01 ) ( 2.30e+01 , 2.30e+01 )
-( 2.40e+01 , 2.40e+01 )
-
-++++++++++vcB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( -1.00e+00 , 1.00e+00 ) ( 0.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 1.00e+00 )
-( 2.00e+00 , 2.00e+00 )
-
---------------------------------------------------J = 1
-( -6.00e+00 , 6.00e+00 ) ( -5.00e+00 , 5.00e+00 )
-( 5.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 )
-
---------------------------------------------------J = 2
-( -1.10e+01 , 1.10e+01 ) ( -1.00e+01 , 1.00e+01 )
-( 1.00e+01 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 )
-
---------------------------------------------------J = 3
-( -1.60e+01 , 1.60e+01 ) ( -1.50e+01 , 1.50e+01 )
-( 1.50e+01 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 )
-
---------------------------------------------------J = 4
-( -2.10e+01 , 2.10e+01 ) ( -2.00e+01 , 2.00e+01 )
-( 2.00e+01 , 2.00e+01 ) ( 2.10e+01 , 2.10e+01 )
-( 2.20e+01 , 2.20e+01 )
-
-++++++++++sA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1, 0:4:1) ~~~~~~~~
-==================================================K = 0
---------------------------------------------------J = 0
-(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 1
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 2
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 3
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 4
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
-==================================================K = 1
---------------------------------------------------J = 0
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 1
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 2
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 3
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 4
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
-==================================================K = 2
---------------------------------------------------J = 0
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 1
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 2
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 3
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 4
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
-==================================================K = 3
---------------------------------------------------J = 0
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 1
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 2
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 3
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
---------------------------------------------------J = 4
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)
-
-
-==================================================K = 4
---------------------------------------------------J = 0
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 1
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 2
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
---------------------------------------------------J = 3
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)
-
-
---------------------------------------------------J = 4
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)
-(1.20e+01 , 1.20e+01 , 1.20e+01)(1.20e+01 , 1.20e+01 , 1.20e+01)(1.20e+01 , 1.20e+01 , 1.20e+01)
-
-
-++++++++++sB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1, 0:4:1) ~~~~~~~~
-==================================================K = 0
---------------------------------------------------J = 0
-(1.00e+00 , -1.00e+00 , -1.00e+00)(-1.00e+00 , 1.00e+00 , -1.00e+00)(-1.00e+00 , -1.00e+00 , 1.00e+00)
-(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)
-(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-
-
---------------------------------------------------J = 1
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(1.00e+00 , -1.00e+00 , -1.00e+00)(-1.00e+00 , 1.00e+00 , -1.00e+00)(-1.00e+00 , -1.00e+00 , 1.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-
-
---------------------------------------------------J = 2
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 3
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 4
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
-==================================================K = 1
---------------------------------------------------J = 0
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(1.00e+00 , -1.00e+00 , -1.00e+00)(-1.00e+00 , 1.00e+00 , -1.00e+00)(-1.00e+00 , -1.00e+00 , 1.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-
-
---------------------------------------------------J = 1
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 2
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 3
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 4
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
-==================================================K = 2
---------------------------------------------------J = 0
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 1
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 2
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 3
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 4
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
-==================================================K = 3
---------------------------------------------------J = 0
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 1
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 2
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 3
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 4
-(8.00e+00 , -8.00e+00 , -8.00e+00)(-8.00e+00 , 8.00e+00 , -8.00e+00)(-8.00e+00 , -8.00e+00 , 8.00e+00)
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
-==================================================K = 4
---------------------------------------------------J = 0
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 1
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 2
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 3
-(8.00e+00 , -8.00e+00 , -8.00e+00)(-8.00e+00 , 8.00e+00 , -8.00e+00)(-8.00e+00 , -8.00e+00 , 8.00e+00)
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 4
-(9.00e+00 , -9.00e+00 , -9.00e+00)(-9.00e+00 , 9.00e+00 , -9.00e+00)(-9.00e+00 , -9.00e+00 , 9.00e+00)
-(8.00e+00 , -8.00e+00 , -8.00e+00)(-8.00e+00 , 8.00e+00 , -8.00e+00)(-8.00e+00 , -8.00e+00 , 8.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
-++++++++++vfzA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 2.00e+00 , 0.00e+00 ) ( 3.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 ) ( 8.00e+00 , 8.00e+00 )
-( 0.00e+00 , 8.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 ) ( 1.30e+01 , 1.30e+01 )
-( 0.00e+00 , 1.30e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 ) ( 1.80e+01 , 1.80e+01 )
-( 0.00e+00 , 1.80e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( 1.60e+01 , 0.00e+00 )
-( 1.70e+01 , 0.00e+00 ) ( 1.80e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vfzB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 6.00e+00 ) ( -6.00e+00 , 5.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 6.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.10e+01 ) ( -1.10e+01 , 1.00e+01 )
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 0.00e+00 , 1.10e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.60e+01 ) ( -1.60e+01 , 1.50e+01 )
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 0.00e+00 , 1.60e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -1.60e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.60e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vfzC+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , -6.00e+00 ) ( -1.00e+00 , -5.00e+00 )
-( 0.00e+00 , -5.00e+00 ) ( 1.00e+00 , -6.00e+00 )
-( 0.00e+00 , -6.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 6.00e+00 ) ( -6.00e+00 , 5.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 6.00e+00 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vzA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 2.00e+00 , 0.00e+00 ) ( 3.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 ) ( 8.00e+00 , 8.00e+00 )
-( 0.00e+00 , 9.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 ) ( 1.30e+01 , 1.30e+01 )
-( 0.00e+00 , 1.40e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 ) ( 1.80e+01 , 1.80e+01 )
-( 0.00e+00 , 1.90e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( 2.10e+01 , 0.00e+00 )
-( 2.20e+01 , 0.00e+00 ) ( 2.30e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vzB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 7.00e+00 ) ( -6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 7.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.20e+01 ) ( -1.10e+01 , 1.10e+01 )
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 0.00e+00 , 1.20e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.70e+01 ) ( -1.60e+01 , 1.60e+01 )
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 0.00e+00 , 1.70e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -2.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 2.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vzC+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , -7.00e+00 ) ( -6.00e+00 , -6.00e+00 )
-( 0.00e+00 , -5.00e+00 ) ( 6.00e+00 , -6.00e+00 )
-( 0.00e+00 , -7.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 7.00e+00 ) ( -6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 7.00e+00 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-
diff --git a/ippl/test/simple/text.test.TestBC b/ippl/test/simple/text.test.TestBC
deleted file mode 100755
index 2191f15f6..000000000
--- a/ippl/test/simple/text.test.TestBC
+++ /dev/null
@@ -1,708 +0,0 @@
-++++++++BConds object cbc begin++++++++
-BConds:(
-NegReflectFace, Face=0 ,
-ZeroFace, Face=1 ,
-ParallelPeriodicFace, Face=2 ,
-ParallelPeriodicFace, Face=3
-)
-
-++++++++BConds object cbc end++++++++++
-++++++++++cA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00
-
---------------------------------------------------J = 1
-5.000e+00 6.000e+00 7.000e+00 8.000e+00 9.000e+00
-
---------------------------------------------------J = 2
-1.000e+01 1.100e+01 1.200e+01 1.300e+01 1.400e+01
-
---------------------------------------------------J = 3
-1.500e+01 1.600e+01 1.700e+01 1.800e+01 1.900e+01
-
---------------------------------------------------J = 4
-2.000e+01 2.100e+01 2.200e+01 2.300e+01 2.400e+01
-
-++++++++++cB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
--1.600e+01 -1.500e+01 1.500e+01 1.600e+01 1.700e+01
-
---------------------------------------------------J = 1
--2.100e+01 -2.000e+01 2.000e+01 2.100e+01 2.200e+01
-
---------------------------------------------------J = 2
--1.000e+00 0.000e+00 0.000e+00 1.000e+00 2.000e+00
-
---------------------------------------------------J = 3
--6.000e+00 -5.000e+00 5.000e+00 6.000e+00 7.000e+00
-
---------------------------------------------------J = 4
--1.100e+01 -1.000e+01 1.000e+01 1.100e+01 1.200e+01
-
-++++++++++vA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00
-
---------------------------------------------------J = 1
-5.000e+00 6.000e+00 7.000e+00 8.000e+00 9.000e+00
-
---------------------------------------------------J = 2
-1.000e+01 1.100e+01 1.200e+01 1.300e+01 1.400e+01
-
---------------------------------------------------J = 3
-1.500e+01 1.600e+01 1.700e+01 1.800e+01 1.900e+01
-
---------------------------------------------------J = 4
-0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00
-
-++++++++++vB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
--1.200e+01 -1.100e+01 1.000e+01 1.100e+01 1.200e+01
-
---------------------------------------------------J = 1
--1.700e+01 -1.600e+01 1.500e+01 1.600e+01 1.700e+01
-
---------------------------------------------------J = 2
--2.000e+00 -1.000e+00 0.000e+00 1.000e+00 2.000e+00
-
---------------------------------------------------J = 3
--7.000e+00 -6.000e+00 5.000e+00 6.000e+00 7.000e+00
-
---------------------------------------------------J = 4
--1.200e+01 -1.100e+01 1.000e+01 1.100e+01 1.200e+01
-
-++++++++++vcA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 1.00e+00 )
-( 2.00e+00 , 2.00e+00 ) ( 3.00e+00 , 3.00e+00 )
-( 4.00e+00 , 4.00e+00 )
-
---------------------------------------------------J = 1
-( 5.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 ) ( 8.00e+00 , 8.00e+00 )
-( 9.00e+00 , 9.00e+00 )
-
---------------------------------------------------J = 2
-( 1.00e+01 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 ) ( 1.30e+01 , 1.30e+01 )
-( 1.40e+01 , 1.40e+01 )
-
---------------------------------------------------J = 3
-( 1.50e+01 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 ) ( 1.80e+01 , 1.80e+01 )
-( 1.90e+01 , 1.90e+01 )
-
---------------------------------------------------J = 4
-( 2.00e+01 , 2.00e+01 ) ( 2.10e+01 , 2.10e+01 )
-( 2.20e+01 , 2.20e+01 ) ( 2.30e+01 , 2.30e+01 )
-( 2.40e+01 , 2.40e+01 )
-
-++++++++++vcB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( -1.00e+00 , 1.00e+00 ) ( 0.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 1.00e+00 )
-( 2.00e+00 , 2.00e+00 )
-
---------------------------------------------------J = 1
-( -6.00e+00 , 6.00e+00 ) ( -5.00e+00 , 5.00e+00 )
-( 5.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 )
-
---------------------------------------------------J = 2
-( -1.10e+01 , 1.10e+01 ) ( -1.00e+01 , 1.00e+01 )
-( 1.00e+01 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 )
-
---------------------------------------------------J = 3
-( -1.60e+01 , 1.60e+01 ) ( -1.50e+01 , 1.50e+01 )
-( 1.50e+01 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 )
-
---------------------------------------------------J = 4
-( -2.10e+01 , 2.10e+01 ) ( -2.00e+01 , 2.00e+01 )
-( 2.00e+01 , 2.00e+01 ) ( 2.10e+01 , 2.10e+01 )
-( 2.20e+01 , 2.20e+01 )
-
-++++++++++sA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1, 0:4:1) ~~~~~~~~
-==================================================K = 0
---------------------------------------------------J = 0
-(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 1
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 2
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 3
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 4
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
-==================================================K = 1
---------------------------------------------------J = 0
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 1
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 2
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 3
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 4
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
-==================================================K = 2
---------------------------------------------------J = 0
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 1
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 2
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 3
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 4
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
-==================================================K = 3
---------------------------------------------------J = 0
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 1
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 2
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 3
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
---------------------------------------------------J = 4
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)
-
-
-==================================================K = 4
---------------------------------------------------J = 0
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 1
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 2
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
---------------------------------------------------J = 3
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)
-
-
---------------------------------------------------J = 4
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)(1.10e+01 , 1.10e+01 , 1.10e+01)
-(1.20e+01 , 1.20e+01 , 1.20e+01)(1.20e+01 , 1.20e+01 , 1.20e+01)(1.20e+01 , 1.20e+01 , 1.20e+01)
-
-
-++++++++++sB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1, 0:4:1) ~~~~~~~~
-==================================================K = 0
---------------------------------------------------J = 0
-(1.00e+00 , -1.00e+00 , -1.00e+00)(-1.00e+00 , 1.00e+00 , -1.00e+00)(-1.00e+00 , -1.00e+00 , 1.00e+00)
-(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)
-(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)(0.00e+00 , 0.00e+00 , 0.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-
-
---------------------------------------------------J = 1
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(1.00e+00 , -1.00e+00 , -1.00e+00)(-1.00e+00 , 1.00e+00 , -1.00e+00)(-1.00e+00 , -1.00e+00 , 1.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-
-
---------------------------------------------------J = 2
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 3
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 4
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
-==================================================K = 1
---------------------------------------------------J = 0
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(1.00e+00 , -1.00e+00 , -1.00e+00)(-1.00e+00 , 1.00e+00 , -1.00e+00)(-1.00e+00 , -1.00e+00 , 1.00e+00)
-(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)(1.00e+00 , 1.00e+00 , 1.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-
-
---------------------------------------------------J = 1
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 2
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 3
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 4
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
-==================================================K = 2
---------------------------------------------------J = 0
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(2.00e+00 , -2.00e+00 , -2.00e+00)(-2.00e+00 , 2.00e+00 , -2.00e+00)(-2.00e+00 , -2.00e+00 , 2.00e+00)
-(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)(2.00e+00 , 2.00e+00 , 2.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-
-
---------------------------------------------------J = 1
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 2
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 3
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 4
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
-==================================================K = 3
---------------------------------------------------J = 0
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(3.00e+00 , -3.00e+00 , -3.00e+00)(-3.00e+00 , 3.00e+00 , -3.00e+00)(-3.00e+00 , -3.00e+00 , 3.00e+00)
-(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)(3.00e+00 , 3.00e+00 , 3.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-
-
---------------------------------------------------J = 1
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 2
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 3
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 4
-(8.00e+00 , -8.00e+00 , -8.00e+00)(-8.00e+00 , 8.00e+00 , -8.00e+00)(-8.00e+00 , -8.00e+00 , 8.00e+00)
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
-==================================================K = 4
---------------------------------------------------J = 0
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(4.00e+00 , -4.00e+00 , -4.00e+00)(-4.00e+00 , 4.00e+00 , -4.00e+00)(-4.00e+00 , -4.00e+00 , 4.00e+00)
-(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)(4.00e+00 , 4.00e+00 , 4.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-
-
---------------------------------------------------J = 1
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(5.00e+00 , -5.00e+00 , -5.00e+00)(-5.00e+00 , 5.00e+00 , -5.00e+00)(-5.00e+00 , -5.00e+00 , 5.00e+00)
-(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)(5.00e+00 , 5.00e+00 , 5.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-
-
---------------------------------------------------J = 2
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(6.00e+00 , -6.00e+00 , -6.00e+00)(-6.00e+00 , 6.00e+00 , -6.00e+00)(-6.00e+00 , -6.00e+00 , 6.00e+00)
-(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)(6.00e+00 , 6.00e+00 , 6.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-
-
---------------------------------------------------J = 3
-(8.00e+00 , -8.00e+00 , -8.00e+00)(-8.00e+00 , 8.00e+00 , -8.00e+00)(-8.00e+00 , -8.00e+00 , 8.00e+00)
-(7.00e+00 , -7.00e+00 , -7.00e+00)(-7.00e+00 , 7.00e+00 , -7.00e+00)(-7.00e+00 , -7.00e+00 , 7.00e+00)
-(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)(7.00e+00 , 7.00e+00 , 7.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-
-
---------------------------------------------------J = 4
-(9.00e+00 , -9.00e+00 , -9.00e+00)(-9.00e+00 , 9.00e+00 , -9.00e+00)(-9.00e+00 , -9.00e+00 , 9.00e+00)
-(8.00e+00 , -8.00e+00 , -8.00e+00)(-8.00e+00 , 8.00e+00 , -8.00e+00)(-8.00e+00 , -8.00e+00 , 8.00e+00)
-(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)(8.00e+00 , 8.00e+00 , 8.00e+00)
-(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)(9.00e+00 , 9.00e+00 , 9.00e+00)
-(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)(1.00e+01 , 1.00e+01 , 1.00e+01)
-
-
-++++++++++vfzA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 2.00e+00 , 0.00e+00 ) ( 3.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 ) ( 8.00e+00 , 8.00e+00 )
-( 0.00e+00 , 8.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 ) ( 1.30e+01 , 1.30e+01 )
-( 0.00e+00 , 1.30e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 ) ( 1.80e+01 , 1.80e+01 )
-( 0.00e+00 , 1.80e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( 1.60e+01 , 0.00e+00 )
-( 1.70e+01 , 0.00e+00 ) ( 1.80e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vfzB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 6.00e+00 ) ( -6.00e+00 , 5.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 6.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.10e+01 ) ( -1.10e+01 , 1.00e+01 )
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 0.00e+00 , 1.10e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.60e+01 ) ( -1.60e+01 , 1.50e+01 )
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 0.00e+00 , 1.60e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -1.60e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.60e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vfzC+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , -6.00e+00 ) ( -1.00e+00 , -5.00e+00 )
-( 0.00e+00 , -5.00e+00 ) ( 1.00e+00 , -6.00e+00 )
-( 0.00e+00 , -6.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 6.00e+00 ) ( -6.00e+00 , 5.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 6.00e+00 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 6.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vzA+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 2.00e+00 , 0.00e+00 ) ( 3.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 7.00e+00 , 7.00e+00 ) ( 8.00e+00 , 8.00e+00 )
-( 0.00e+00 , 9.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 1.20e+01 , 1.20e+01 ) ( 1.30e+01 , 1.30e+01 )
-( 0.00e+00 , 1.40e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 1.70e+01 , 1.70e+01 ) ( 1.80e+01 , 1.80e+01 )
-( 0.00e+00 , 1.90e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( 2.10e+01 , 0.00e+00 )
-( 2.20e+01 , 0.00e+00 ) ( 2.30e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vzB+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , 7.00e+00 ) ( -6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 7.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 1.20e+01 ) ( -1.10e+01 , 1.10e+01 )
-( 0.00e+00 , 1.00e+01 ) ( 1.10e+01 , 1.10e+01 )
-( 0.00e+00 , 1.20e+01 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 1.70e+01 ) ( -1.60e+01 , 1.60e+01 )
-( 0.00e+00 , 1.50e+01 ) ( 1.60e+01 , 1.60e+01 )
-( 0.00e+00 , 1.70e+01 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -2.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 2.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-++++++++++vzC+++++++++++
-~~~~~~~~ field slice (0:4:1, 0:4:1) ~~~~~~~~
---------------------------------------------------J = 0
-( 0.00e+00 , 0.00e+00 ) ( -1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 1
-( 0.00e+00 , -7.00e+00 ) ( -6.00e+00 , -6.00e+00 )
-( 0.00e+00 , -5.00e+00 ) ( 6.00e+00 , -6.00e+00 )
-( 0.00e+00 , -7.00e+00 )
-
---------------------------------------------------J = 2
-( 0.00e+00 , 0.00e+00 ) ( -1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.00e+00 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
---------------------------------------------------J = 3
-( 0.00e+00 , 7.00e+00 ) ( -6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 5.00e+00 ) ( 6.00e+00 , 6.00e+00 )
-( 0.00e+00 , 7.00e+00 )
-
---------------------------------------------------J = 4
-( 0.00e+00 , 0.00e+00 ) ( -1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 ) ( 1.10e+01 , 0.00e+00 )
-( 0.00e+00 , 0.00e+00 )
-
-
diff --git a/ippl/test/simple/vtest1 b/ippl/test/simple/vtest1
deleted file mode 100755
index 0d418e0f1..000000000
--- a/ippl/test/simple/vtest1
+++ /dev/null
@@ -1,12 +0,0 @@
-vnode = 0
- [-1][-1]=( 0 , 0 ) [0][-1]=( 0 , 0 ) [1][-1]=( 0 , 0 ) [2][-1]=( 0 , 0 ) [3][-1]=( 0 , 0 ) [4][-1]=( 0 , 0 ) [5][-1]=( 0 , 0 ) [6][-1]=( 0 , 0 ) [7][-1]=( 0 , 0 ) [8][-1]=( 0 , 0 )
- [-1][0]=( 0 , 0 ) [0][0]=( 0 , 0 ) [1][0]=( 1 , 0 ) [2][0]=( 2 , 0 ) [3][0]=( 3 , 0 ) [4][0]=( 4 , 0 ) [5][0]=( 5 , 0 ) [6][0]=( 6 , 0 ) [7][0]=( 7 , 0 ) [8][0]=( 0 , 0 )
- [-1][1]=( 0 , 0 ) [0][1]=( 0 , 1 ) [1][1]=( 1 , 1 ) [2][1]=( 2 , 1 ) [3][1]=( 3 , 1 ) [4][1]=( 4 , 1 ) [5][1]=( 5 , 1 ) [6][1]=( 6 , 1 ) [7][1]=( 7 , 1 ) [8][1]=( 0 , 0 )
- [-1][2]=( 0 , 0 ) [0][2]=( 0 , 2 ) [1][2]=( 1 , 2 ) [2][2]=( 2 , 2 ) [3][2]=( 3 , 2 ) [4][2]=( 4 , 2 ) [5][2]=( 5 , 2 ) [6][2]=( 6 , 2 ) [7][2]=( 7 , 2 ) [8][2]=( 0 , 0 )
- [-1][3]=( 0 , 0 ) [0][3]=( 0 , 3 ) [1][3]=( 1 , 3 ) [2][3]=( 2 , 3 ) [3][3]=( 3 , 3 ) [4][3]=( 4 , 3 ) [5][3]=( 5 , 3 ) [6][3]=( 6 , 3 ) [7][3]=( 7 , 3 ) [8][3]=( 0 , 0 )
- [-1][4]=( 0 , 0 ) [0][4]=( 0 , 4 ) [1][4]=( 1 , 4 ) [2][4]=( 2 , 4 ) [3][4]=( 3 , 4 ) [4][4]=( 4 , 4 ) [5][4]=( 5 , 4 ) [6][4]=( 6 , 4 ) [7][4]=( 7 , 4 ) [8][4]=( 0 , 0 )
- [-1][5]=( 0 , 0 ) [0][5]=( 0 , 5 ) [1][5]=( 1 , 5 ) [2][5]=( 2 , 5 ) [3][5]=( 3 , 5 ) [4][5]=( 4 , 5 ) [5][5]=( 5 , 5 ) [6][5]=( 6 , 5 ) [7][5]=( 7 , 5 ) [8][5]=( 0 , 0 )
- [-1][6]=( 0 , 0 ) [0][6]=( 0 , 6 ) [1][6]=( 1 , 6 ) [2][6]=( 2 , 6 ) [3][6]=( 3 , 6 ) [4][6]=( 4 , 6 ) [5][6]=( 5 , 6 ) [6][6]=( 6 , 6 ) [7][6]=( 7 , 6 ) [8][6]=( 0 , 0 )
- [-1][7]=( 0 , 0 ) [0][7]=( 0 , 7 ) [1][7]=( 1 , 7 ) [2][7]=( 2 , 7 ) [3][7]=( 3 , 7 ) [4][7]=( 4 , 7 ) [5][7]=( 5 , 7 ) [6][7]=( 6 , 7 ) [7][7]=( 7 , 7 ) [8][7]=( 0 , 0 )
- [-1][8]=( 0 , 0 ) [0][8]=( 0 , 0 ) [1][8]=( 0 , 0 ) [2][8]=( 0 , 0 ) [3][8]=( 0 , 0 ) [4][8]=( 0 , 0 ) [5][8]=( 0 , 0 ) [6][8]=( 0 , 0 ) [7][8]=( 0 , 0 ) [8][8]=( 0 , 0 )
-------------------
diff --git a/ippl/test/simple/vtest2 b/ippl/test/simple/vtest2
deleted file mode 100755
index 0d418e0f1..000000000
--- a/ippl/test/simple/vtest2
+++ /dev/null
@@ -1,12 +0,0 @@
-vnode = 0
- [-1][-1]=( 0 , 0 ) [0][-1]=( 0 , 0 ) [1][-1]=( 0 , 0 ) [2][-1]=( 0 , 0 ) [3][-1]=( 0 , 0 ) [4][-1]=( 0 , 0 ) [5][-1]=( 0 , 0 ) [6][-1]=( 0 , 0 ) [7][-1]=( 0 , 0 ) [8][-1]=( 0 , 0 )
- [-1][0]=( 0 , 0 ) [0][0]=( 0 , 0 ) [1][0]=( 1 , 0 ) [2][0]=( 2 , 0 ) [3][0]=( 3 , 0 ) [4][0]=( 4 , 0 ) [5][0]=( 5 , 0 ) [6][0]=( 6 , 0 ) [7][0]=( 7 , 0 ) [8][0]=( 0 , 0 )
- [-1][1]=( 0 , 0 ) [0][1]=( 0 , 1 ) [1][1]=( 1 , 1 ) [2][1]=( 2 , 1 ) [3][1]=( 3 , 1 ) [4][1]=( 4 , 1 ) [5][1]=( 5 , 1 ) [6][1]=( 6 , 1 ) [7][1]=( 7 , 1 ) [8][1]=( 0 , 0 )
- [-1][2]=( 0 , 0 ) [0][2]=( 0 , 2 ) [1][2]=( 1 , 2 ) [2][2]=( 2 , 2 ) [3][2]=( 3 , 2 ) [4][2]=( 4 , 2 ) [5][2]=( 5 , 2 ) [6][2]=( 6 , 2 ) [7][2]=( 7 , 2 ) [8][2]=( 0 , 0 )
- [-1][3]=( 0 , 0 ) [0][3]=( 0 , 3 ) [1][3]=( 1 , 3 ) [2][3]=( 2 , 3 ) [3][3]=( 3 , 3 ) [4][3]=( 4 , 3 ) [5][3]=( 5 , 3 ) [6][3]=( 6 , 3 ) [7][3]=( 7 , 3 ) [8][3]=( 0 , 0 )
- [-1][4]=( 0 , 0 ) [0][4]=( 0 , 4 ) [1][4]=( 1 , 4 ) [2][4]=( 2 , 4 ) [3][4]=( 3 , 4 ) [4][4]=( 4 , 4 ) [5][4]=( 5 , 4 ) [6][4]=( 6 , 4 ) [7][4]=( 7 , 4 ) [8][4]=( 0 , 0 )
- [-1][5]=( 0 , 0 ) [0][5]=( 0 , 5 ) [1][5]=( 1 , 5 ) [2][5]=( 2 , 5 ) [3][5]=( 3 , 5 ) [4][5]=( 4 , 5 ) [5][5]=( 5 , 5 ) [6][5]=( 6 , 5 ) [7][5]=( 7 , 5 ) [8][5]=( 0 , 0 )
- [-1][6]=( 0 , 0 ) [0][6]=( 0 , 6 ) [1][6]=( 1 , 6 ) [2][6]=( 2 , 6 ) [3][6]=( 3 , 6 ) [4][6]=( 4 , 6 ) [5][6]=( 5 , 6 ) [6][6]=( 6 , 6 ) [7][6]=( 7 , 6 ) [8][6]=( 0 , 0 )
- [-1][7]=( 0 , 0 ) [0][7]=( 0 , 7 ) [1][7]=( 1 , 7 ) [2][7]=( 2 , 7 ) [3][7]=( 3 , 7 ) [4][7]=( 4 , 7 ) [5][7]=( 5 , 7 ) [6][7]=( 6 , 7 ) [7][7]=( 7 , 7 ) [8][7]=( 0 , 0 )
- [-1][8]=( 0 , 0 ) [0][8]=( 0 , 0 ) [1][8]=( 0 , 0 ) [2][8]=( 0 , 0 ) [3][8]=( 0 , 0 ) [4][8]=( 0 , 0 ) [5][8]=( 0 , 0 ) [6][8]=( 0 , 0 ) [7][8]=( 0 , 0 ) [8][8]=( 0 , 0 )
-------------------
diff --git a/ippl/test/sindex/compsindex.cpp b/ippl/test/sindex/compsindex.cpp
deleted file mode 100644
index 7a1876e9f..000000000
--- a/ippl/test/sindex/compsindex.cpp
+++ /dev/null
@@ -1,169 +0,0 @@
-// -*- C++ -*-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by PSI.
- * All rights in the program are reserved by PSI.
- * Neither PSI nor the author(s)
- * makes any warranty, express or implied, or assumes any liability or
- * responsibility for the use of this software
- *
- * Visit http://www.acl.lanl.gov/POOMS for more details
- *
- ***************************************************************************/
-
-// -*- C++ -*-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by PSI.
- * All rights in the program are reserved by PSI.
- * Neither PSI nor the author(s)
- * makes any warranty, express or implied, or assumes any liability or
- * responsibility for the use of this software
- *
- * Visit http://www.acl.lanl.gov/POOMS for more details
- *
- ***************************************************************************/
-
-#include "Ippl.h"
-
-// Benchmarks where() against SIndex for a particular example, where sparsness
-// and compression should be high. For testing Bill's new code making SIndex
-// compression-smarter.
-
-//-----------------------------------------------------------------------------
-// User-inserted prototypes to get debugger access (examples for user ref.):
-#ifdef __MWERKS__
-// Avoid bogus warning about no prototypes in CW4:
-void dfp3(BareField<double,3U> &f);
-void defp3(BareField<double,3U> &f, int i, int j, int k);
-void dsfp3(BareField<double,3U> &f, int,int,int, int,int,int, int,int,int);
-#endif // __MWERKS__
-// Scalar (double):
-void dfp3(BareField<double,3U> &f) {fp3(f);}
-void defp3(BareField<double,3U> &f, int i, int j, int k) {efp3(f,i,j,k);}
-void dsfp3(BareField<double,3U> &f,
- int base1, int bound1, int stride1,
- int base2, int bound2, int stride2,
- int base3, int bound3, int stride3) {
- sfp3(f,base1,bound1,stride1,base2,bound2,stride2,base3,bound3,stride3);}
-//-----------------------------------------------------------------------------
-
-
-int main(int argc, char *argv[]) {
-
- Ippl ippl(argc,argv);
- Inform pout(NULL,0);
-
- Timer t; t.start(); t.stop(); t.clear();
- float tcpu, tclk;
- float tcreatecpu, tcreateclk;
- float tsetupcpu, tsetupclk;
- float tassigncpu, tassignclk;
- float time1, time2;
-
- // Various counters, constants, etc:
- const unsigned D=3U;
- double hightail;
- int d;
- int tag = Ippl::Comm->next_tag(tag);
- unsigned ngrid[D]; // grid sizes sizes
-
- // Layout information:
- unsigned vnodes; // number of vnodes; input at cmd line
- e_dim_tag serialParallel[D]; // Specifies SERIAL, PARALLEL dims
- for (d=0; d<D; d++) serialParallel[d] = PARALLEL;
-
- // read in vnodes etc. from command line
- if (argc != (D+3)) {
- pout << "Usage: " << argv[0] << " <nx> <ny> <nz> <vnodes> <percent>";
- pout << "\n, where <percent> = percentage of field not used in where";
- pout << endl;
- } else {
- for (d=0; d < D; ++d)
- ngrid[d] = atoi(argv[d+1]);
- vnodes = atoi(argv[d+1]);
- hightail = 0.01 * atof(argv[d+2]);
- Index I(ngrid[0]), J(ngrid[1]), K(ngrid[2]);
- NDIndex<3U> ndi(I,J,K);
- FieldLayout<D> layout(ndi, serialParallel, vnodes);
- Field<double,D> A(layout), B(layout);
-
- // time how long the SIndex creation takes
- t.start();
- SIndex<D> SI(layout);
- t.stop(); tcreatecpu=t.cpu_time(); tcreateclk=t.clock_time(); t.clear();
-
- // Assign A to be random in the lower eighth of the 3D box, zero elsewhere,
- // which gives good compression with enough vnodes:
- A = 0.0;
- NDIndex<D> lowerEighth;
- for (d=0; d<D; d++) lowerEighth[d] = Index(ngrid[d]/8);
- A[lowerEighth] = IpplRandom;
- pout << endl << "A.CompressedFraction() = "
- << A.CompressedFraction() << endl;
-
- // Assign B to zero, so it starts fully compressed; should next uncompress
- // only lower eighth of it:
- B = 0.0;
-
- // pout << "sum of A at start = " << sum(A) << endl;
- // pout << "sum of B at start = " << sum(B) << endl;
-
- // time how long the SIndex assignment takes
- t.start();
- SI = gt(A, hightail);
- t.stop(); tsetupcpu=t.cpu_time(); tsetupclk=t.clock_time(); t.clear();
-
- // time how long a regular where-statement assignment takes
- t.start();
- B[I][J][K] = where((gt(A[I][J][K],hightail)), A[I][J][K], B[I][J][K]);
- t.stop(); tcpu= t.cpu_time(); tclk=t.clock_time(); t.clear();
- time1 = tclk;
-
- pout << endl << "time for where(gt(A[I][J][K]," << hightail
- << ")) assign = " << tclk << " clock, "
- << tcpu << " cpu " << endl;
- pout << "sum of B after where assign = " << sum(B) << endl;
-
- // Reassign B to zero, so it starts fully compressed; should next
- // uncompress only lower eighth of it:
- B = 0.0;
-
- // time how long a sparse index field assignment takes
- t.start();
- B[SI] = A[SI];
- t.stop(); tassigncpu=t.cpu_time(); tassignclk=t.clock_time(); t.clear();
- time2 = tcreateclk + tsetupclk + tassignclk;
-
- pout << endl << "time for SIndex creation = " << tcreateclk;
- pout << " clock, " << tcreatecpu << " cpu " << endl;
- pout << "time for SIndex setup (SI = gt(A," << hightail << ")) = "
- << tsetupclk;
- pout << " clock, " << tsetupcpu << " cpu " << endl;
- pout << "time for [SIndex] assign " << tassignclk << " clock, "
- << tassigncpu << " cpu " << endl;
- pout << "sum of B after SIndex assign = " << sum(B) << endl;
-
- pout << endl << "========================================" << endl << endl;
-
- pout << "where() timing: " <<time1 << " clock" << endl;
-
- pout << "SIndex timing: create+setup+assign: " <<time2 << " clock" << endl;
-
- pout << "where()/SIndex Ratio: "
- << (time2 > 0.0 ? (time1/time2) : 0.0) << endl << endl;
- }
-
- return 0;
-}
-
-/***************************************************************************
- * $RCSfile: addheaderfooter,v $ $Author: adelmann $
- * $Revision: 1.1.1.1 $ $Date: 2003/01/23 07:40:17 $
- * IPPL_VERSION_ID: $Id: addheaderfooter,v 1.1.1.1 2003/01/23 07:40:17 adelmann Exp $
- ***************************************************************************/
-
diff --git a/ippl/test/sindex/sindex.cpp b/ippl/test/sindex/sindex.cpp
deleted file mode 100644
index 8a4d63a41..000000000
--- a/ippl/test/sindex/sindex.cpp
+++ /dev/null
@@ -1,120 +0,0 @@
-// -*- C++ -*-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by PSI.
- * All rights in the program are reserved by PSI.
- * Neither PSI nor the author(s)
- * makes any warranty, express or implied, or assumes any liability or
- * responsibility for the use of this software
- *
- * Visit http://www.acl.lanl.gov/POOMS for more details
- *
- ***************************************************************************/
-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by the Regents of the University of
- * Visit http://people.web.psi.ch/adelmann/ for more details
- *
- ***************************************************************************/
-
-#include "Ippl.h"
-#include "Index/SIndex.h"
-
-/***************************************************************************
- A simple program to test the capabilities of the SIndex and SOffset classes.
- ***************************************************************************/
-
-int main(int argc, char *argv[]) {
-
- Ippl ippl(argc,argv);
- Inform testmsg(argv[0], INFORM_ALL_NODES);
-
- const unsigned Dim=2;
- Index I(25);
- Index J(50);
- Index I2(2);
- Index J2(5);
- NDIndex<Dim> NDX(I,J);
- int IP[Dim];
- IP[0] = -1;
- IP[1] = -1;
-
- SOffset<Dim> so1(1,1);
- SOffset<Dim> so2(IP);
- testmsg << "Created SOffset so1 = " << so1 << endl;
- testmsg << "Created SOffset so2 = " << so2 << endl;
- testmsg << "Adding IP to so1 = " << so1 + IP << endl;
- testmsg << "Adding so1 to IP = " << IP + so1 << endl;
- testmsg << "Adding so2 to so1 = " << so1 + so2 << endl;
- testmsg << "Adding so1 to NDX " << NDX << " = " << NDX + so1 << endl;
- testmsg << "Adding NDX " << NDX << " to so1 = " << so1 + NDX << endl;
- testmsg << "Multipplying NDX " << NDX << " by so2 = " << NDX * so2 << endl;
- testmsg << "Multipplying so2 by NDX " << NDX << " = " << so2 * NDX << endl;
- so1 -= IP;
- so2 += IP;
- testmsg << "Accumulated IP from so1 = " << so1 << endl;
- testmsg << "Accumulated IP into so2 = " << so2 << endl;
-
- FieldLayout<Dim> layout(I, J, PARALLEL, PARALLEL, 2*Ippl::getNodes());
- Field<double,Dim> A(layout);
- Field<bool,Dim> B(layout);
-
- SIndex<Dim> s1(layout);
- SIndex<Dim> s2 = s1(1,-1);
- SIndex<Dim> x3 = s1(IP);
-
- testmsg << "Created s1 = " << s1 << endl;
- testmsg << "Created s2 = " << s2 << endl;
- testmsg << "Created x3 = " << x3 << endl;
-
- s1.addIndex(NDIndex<Dim>(Index(2), Index(3)));
- s2.addIndex(SOffset<Dim>(0,0));
- s2.addIndex(NDIndex<Dim>(Index(20,23), Index(45,46)));
-
- testmsg << "Added new points, s1 = "; s1.printDebug(testmsg);
- testmsg << "Added new points, s2 = "; s2.printDebug(testmsg);
-
- SIndex<Dim> s3(s1);
- s3 = NDIndex<Dim>(Index(1,5), Index(2,4));
-
- testmsg << "Created s3 = " << s3 << endl;
-
- s3 &= s1;
-
- testmsg << "Intersection of s3 and s1 = " << s3 << endl;
-
- // now, test assigment of a Field expression to an SIndex
- A[I][J] = I + J;
- B = lt(A,10);
- s3 = lt(A,10);
-
- // do a union with a slightly different condition
- s3 |= (gt(A,15) && lt(A,20));
- testmsg << "union of s3 and where 15 < A < 20 = " << s3 << endl;
-
- // do an intersection with an overlapping condition
- s3 &= (gt(A,8) && lt(A,12));
- testmsg << "intersection of s3 and where 8 < A < 12 = " << s3 << endl;
-
- // do an indexed assignment to s3
- s3[I2][J2] = (lt(A[I2 + 1][J2+1], 5) && gt(A[I2 + 2][J2 + 2], 0));
- testmsg << "s3[I2][J2] = expr ==> s3 = " << s3 << endl;
-
- // now use s3 in a Field expression
- testmsg << "Originally, A = I + J ... after A[s3] = A[s3(1,1)]: A = ";
- A[s3] = A[s3(1,1)] + 10;
- testmsg << A << endl;
-
- return 0;
-}
-
-/***************************************************************************
- * $RCSfile: sindex.cpp,v $ $Author: adelmann $
- * $Revision: 1.1.1.1 $ $Date: 2003/01/23 07:40:40 $
- * IPPL_VERSION_ID: $Id: sindex.cpp,v 1.1.1.1 2003/01/23 07:40:40 adelmann Exp $
- ***************************************************************************/
diff --git a/ippl/test/sindex/subfield.cpp b/ippl/test/sindex/subfield.cpp
deleted file mode 100644
index 694c4839e..000000000
--- a/ippl/test/sindex/subfield.cpp
+++ /dev/null
@@ -1,233 +0,0 @@
-// -*- C++ -*-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by PSI.
- * All rights in the program are reserved by PSI.
- * Neither PSI nor the author(s)
- * makes any warranty, express or implied, or assumes any liability or
- * responsibility for the use of this software
- *
- * Visit http://www.acl.lanl.gov/POOMS for more details
- *
- ***************************************************************************/
-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by the Regents of the University of
- * Visit http://people.web.psi.ch/adelmann/ for more details
- *
- ***************************************************************************/
-
-#include "Ippl.h"
-
-/***************************************************************************
- A simple program to test SubField assignments using sparse index operations
- (with SIndex, SOffset, and NDIndex objects as the subsetting mechanism).
- ***************************************************************************/
-
-int main(int argc, char *argv[]) {
-
- Ippl ippl(argc,argv);
- Inform testmsg(argv[0], INFORM_ALL_NODES);
-
- const unsigned Dim=2;
- Index I(4);
- Index I2(2);
- Index J(4);
- Index K(2);
- FieldLayout<Dim> layout(I, J, PARALLEL, PARALLEL, 4*Ippl::getNodes());
- FieldLayout<3> layout3(I, J, K, PARALLEL, PARALLEL, PARALLEL,
- 2*Ippl::getNodes());
- BConds<double,Dim> dbc;
- for (int f=0; f < 2*Dim; f++) dbc[f] = new ZeroFace<double,Dim>(f);
- Field<double,Dim> A(layout, dbc, GuardCellSizes<Dim>(1));
- A = 0.0;
- BConds<bool,Dim> bbc;
- for (int f=0; f < 2*Dim; f++) bbc[f] = new ZeroFace<bool,Dim>(f);
- Field<bool,Dim> B(layout, bbc, GuardCellSizes<Dim>(1));
- B = true;
- Field<double,Dim> C(layout, dbc, GuardCellSizes<Dim>(1));
- C = 0.0;
-
- BConds<float,3> fbc;
- for (int f=0; f < 2*3; f++) fbc[f] = new ZeroFace<float,3>(f);
- Field<float,3> A3(layout3, fbc, GuardCellSizes<3>(2));
- A3 = 0.0;
- Field<float,3> B3(layout3, fbc, GuardCellSizes<3>(2));
- B3 = 0.0;
- SIndex<Dim> s1(layout);
- s1.addIndex(NDIndex<Dim>(Index(1,1), Index(1,2)));
-
- SIndex<Dim> s2 = s1(1,-1);
- s2.addIndex(SOffset<Dim>(0,0));
-
- testmsg << "\n************ testing 2D SubParticleAttrib<> ************" << endl;
- A[I][J] = (I+1) + (J*10);
- testmsg << "s1 = " << s1 << endl;
- testmsg << " A = " << A << endl;
- ParticleAttrib<double> PA;
- PA[s1] = A[s1];
- testmsg << "Result of PA[s1] = A[s1] : PA = " << endl;
- testmsg << PA[s1] << endl;
-
- testmsg << "s2 = " << s2 << endl;
- PA[s2] = A[s2] + A[s2];
- testmsg << "Result of PA[s2] = A[s2] + A[s2] : PA = " << endl;
- testmsg << PA[s2] << endl;
-
- PA[s2] *= (A[s2] + PA[s2]);
- testmsg << "Result of PA[s2] *= A[s2] + PA[s2] : PA = " << endl;
- testmsg << PA[s2] << endl;
-
- A = 1.0;
- SIndex<Dim> sX = s1;
- sX.addIndex(SOffset<Dim>(0,0));
- testmsg << "Created sX = " << sX << endl;
-
- A[sX] += PA[sX];
- testmsg << "Result of A[sX] = 1 + PA[sX] : A = " << A << endl;
-
- A = 0.0;
- C = 1.0;
- A[sX] = C[sX];
- testmsg << "Result of A[sX] = C[sX], after C = 1.0: A = " << A << endl;
-
- ParticleAttrib<Vektor<double,2> > PA2;
- PA2(0)[sX] = A[sX] + PA[sX];
- PA2(1)[sX] = -(A[sX] + PA[sX]);
- testmsg << "Result of PA2[sX] = +- (A[sX] + PA[sX]):" << endl;
- testmsg << " A = " << A << endl;
- testmsg << " PA = " << PA[sX] << endl;
- testmsg << "PA2 = " << PA2[sX] << endl;
-
- testmsg << "\n************ testing 2D SubBareField<NDIndex> ************" << endl;
-
- Index IX(0,2);
- Index JX(2,2);
- NDIndex<Dim> sub1(IX, JX);
- SubBareField<double,Dim,NDIndex<Dim> > SA1(A, sub1);
- testmsg << "Created SubBareField SA1 from A and NDIndex " << sub1 << ":" <<endl;
- testmsg << " SA1 = " << SA1 << endl;
-
- SA1 = 3.0;
- testmsg << "SA1 set to 3:" << endl;
- testmsg << " A = " << A << endl;
-
- NDIndex<Dim> sub2(IX + 1, JX);
- SubBareField<double,Dim,NDIndex<Dim> > SA2(A, sub2);
- testmsg << "Created SubBareField SA2 from A and NDIndex " << sub2 << ":" <<endl;
- testmsg << " SA2 = " << SA2 << endl;
-
- SA1 = SA2 + 4.0;
- testmsg << "SA1 on " << sub1 << " set to (SA2 + 4):" << endl;
- testmsg << " A = " << A << endl;
-
- testmsg << "\n************ testing 2D SubBareField<SIndex> ************" << endl;
-
- SubBareField<double,Dim,SIndex<Dim> > SB1(A, s1);
- testmsg << "Created SubBareField SB1 from A and SIndex " << s1 << ":" << endl;
- testmsg << " SB1 = " << SB1 << endl;
-
- SB1 = -2.0;
- testmsg << "SB1 set to -2:" << endl;
- testmsg << " A = " << A << endl;
-
- SubBareField<double,Dim,SIndex<Dim> > SB2(A, s1(1,-1));
- testmsg << "Created SubBareField SB2 from A and SIndex " << SB2.getDomain();
- testmsg << ":" << endl;
- testmsg << " SB2 = " << SB2 << endl;
-
- SB1 = SB2 + 4.0;
- testmsg << "SB1 set to (SB2 + 4):" << endl;
- testmsg << " A = " << A << endl;
-
- testmsg << "\n************ testing 2D SubBareField<SOffset> ************" << endl;
-
- SubBareField<double,Dim,SOffset<Dim> > SC1(A, SOffset<Dim>(2,1));
- testmsg << "Created SubBareField SC1 from A and SOffset " << SC1.getDomain();
- testmsg << ":" << endl;
- testmsg << " SC1 = " << SC1 << endl;
-
- SC1 = 10.0;
- testmsg << "SC1 set to 10:" << endl;
- testmsg << " A = " << A << endl;
-
- SubBareField<double,Dim,SOffset<Dim> > SC2(A, SOffset<Dim>(2,2));
- SubBareField<double,Dim,SOffset<Dim> > SC3(A, SOffset<Dim>(3,2));
- SC1 = SC1 + (SC2 + SC3) * 10;
- testmsg << "SC1 set to 10 + [(2,2) + (3,2)]*10:" << endl;
- testmsg << " A = " << A << endl;
-
- testmsg << "\n************ testing 2D Field[SIndex] ************" << endl;
-
- testmsg << "Current SIndex s1 = " << s1 << endl;
- testmsg << "Results of A[s1] = -1000:" << endl;
- A[s1] = -1000.0;
- testmsg << " A = " << A << endl;
- B[s1] = ne(A[s1], A[s1]);
- testmsg << "Results of B[s1] = (A[s1] != A[s1]):" << endl;
- testmsg << " B = " << B << endl;
- C[s1] = 1.0;
- testmsg << "Results of C[s1] = 1:" << endl;
- testmsg << " C = " << C << endl;
- C[s1] = A[s1(-1,-1)] - 2000;
- testmsg << "Results of C[s1] = (A[s1(-1,-1)] - 2000) :" << endl;
- testmsg << " C = " << C << endl;
-
- testmsg << "\n************ testing 2D Field[SIndex,compressed] ****" << endl;
-
- Field<double,Dim> A2(layout);
- A2 = 1.0;
- testmsg << " A2 = " << A2 << endl;
- testmsg << " A2.compressedFraction = " << A2.CompressedFraction() << endl;
- A2[I2][I2] = -10.0;
- testmsg << "Initial settings:" << endl;
- testmsg << " A2 = " << A2 << endl;
- testmsg << " A2.compressedFraction = " << A2.CompressedFraction() << endl;
- s1 = ne(A2, 0.0);
- testmsg << " s1 = " << s1 << endl;
- testmsg << " A2 = " << A2 << endl;
- testmsg << " A2.compressedFraction = " << A2.CompressedFraction() << endl;
- A2[s1] = A2[s1] + 100.0;
- testmsg << "Results of A2[s1] = A2[s1] + 100.0 :" << endl;
- testmsg << " A2 = " << A2 << endl;
- testmsg << " A2.compressedFraction = " << A2.CompressedFraction() << endl;
- testmsg << " s1 = " << s1 << endl;
-
- testmsg << "\n************ testing 3D Field[SIndex] ************" << endl;
-
- A3[I][J][K] = I + J + K;
- B3 = 0.0;
-
- SIndex<3> sindx3(layout3);
- sindx3 = eq(A3, 2.0);
-
- SOffset<3> offset;
- offset[0] = 1;
- offset[1] = 0;
- offset[2] = (-1);
-
- testmsg << "Current 3D SIndex sindx3 = " << sindx3 << endl;
- testmsg << "Current 3D A3 = " << A3 << endl;
- testmsg << "A3 after A3[sindx3] = -1:" << endl;
- A3[sindx3] = -1.0;
- testmsg << " A3 = " << A3 << endl;
- B3[sindx3] = 1.0;
- testmsg << "Results of B3[sindx3] = 1:" << endl;
- testmsg << " B3 = " << B3 << endl;
- B3[sindx3] = A3[sindx3(offset)] - 2000.0;
- testmsg << "Results of B3[sindx3] = A3[sindx3(1,0,-1)] - 2000.0:" << endl;
- testmsg << " B3 = " << B3 << endl;
-
- return 0;
-}
-
-/***************************************************************************
- * $RCSfile: subfield.cpp,v $ $Author: adelmann $
- * $Revision: 1.1.1.1 $ $Date: 2003/01/23 07:40:40 $
- * IPPL_VERSION_ID: $Id: subfield.cpp,v 1.1.1.1 2003/01/23 07:40:40 adelmann Exp $
- ***************************************************************************/
diff --git a/ippl/test/sindex/vertcell.cpp b/ippl/test/sindex/vertcell.cpp
deleted file mode 100644
index 971104e44..000000000
--- a/ippl/test/sindex/vertcell.cpp
+++ /dev/null
@@ -1,66 +0,0 @@
-// -*- C++ -*-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by PSI.
- * All rights in the program are reserved by PSI.
- * Neither PSI nor the author(s)
- * makes any warranty, express or implied, or assumes any liability or
- * responsibility for the use of this software
- *
- * Visit http://www.acl.lanl.gov/POOMS for more details
- *
- ***************************************************************************/
-
-/***************************************************************************
- *
- * The IPPL Framework
- *
- * This program was prepared by the Regents of the University of
- * Visit http://people.web.psi.ch/adelmann/ for more details
- *
- ***************************************************************************/
-
-#include "Ippl.h"
-
-/***************************************************************************
- A simple program to test SubField assignments using two fields, one on
- a cell-centered layout, the other on a vertex-centered (but otherwise
- similar) layout.
- ***************************************************************************/
-
-int main(int argc, char *argv[]) {
-
- Ippl ippl(argc,argv);
- Inform testmsg(argv[0], INFORM_ALL_NODES);
-
- const unsigned Dim=2;
- Index I(4), J(4), U(5), V(5);
- FieldLayout<Dim> LC(I, J, PARALLEL, PARALLEL, 4);
- FieldLayout<Dim> LV(U, V, PARALLEL, PARALLEL, 4);
-
- Field<int,Dim,UniformCartesian<Dim>,Cell> A(LC, GuardCellSizes<Dim>(1));
- Field<int,Dim,UniformCartesian<Dim>,Vert> B(LV, GuardCellSizes<Dim>(1));
-
- A[I][J] = J + 1;
- B[U][V] = -V - 1;
- testmsg << "A at start: " << A << endl;
- testmsg << "B at start: " << B << endl;
-
- SIndex<Dim> SI(LC);
- SI = eq(A,2) || eq(A,4);
- testmsg << "SI at start: " << SI << endl;
-
- A[SI] = A[SI(0,-1)] + B[SI(1,1)];
- testmsg << "A after assignment: " << A << endl;
- testmsg << "B after assignment: " << B << endl;
-
- return 0;
-}
-
-/***************************************************************************
- * $RCSfile: vertcell.cpp,v $ $Author: adelmann $
- * $Revision: 1.1.1.1 $ $Date: 2003/01/23 07:40:40 $
- * IPPL_VERSION_ID: $Id: vertcell.cpp,v 1.1.1.1 2003/01/23 07:40:40 adelmann Exp $
- ***************************************************************************/
diff --git a/tests/ippl_src/Field/BCond.cpp b/tests/ippl_src/Field/BCond.cpp
index 0c3f725e4..ca9c5d5ee 100644
--- a/tests/ippl_src/Field/BCond.cpp
+++ b/tests/ippl_src/Field/BCond.cpp
@@ -12,10 +12,11 @@
#include <fstream>
#include <iostream>
+#include <string>
namespace {
- void hardCodedOutput(char* filename); // Prototype of function defined below.
- bool thediff(char* filename1, char* filename2);
+ void hardCodedOutput(std::string filename); // Prototype of function defined below.
+ bool thediff(std::string filename1, std::string filename2);
}
constexpr unsigned Dim = 2;
@@ -317,7 +318,7 @@ namespace {
//-----------------------------------------------------------------------------
// Mock up the Unix "diff" utility to compare two files:
//-----------------------------------------------------------------------------
-bool thediff(char* filename1, char* filename2)
+bool thediff(std::string filename1, std::string filename2)
{
bool same = true;
char ch1, ch2;
@@ -335,7 +336,7 @@ bool thediff(char* filename1, char* filename2)
}
//-----------------------------------------------------------------------------
-void hardCodedOutput(char* filename)
+void hardCodedOutput(std::string filename)
{
std::ofstream of(filename);
of << "++++++++BConds object cbc begin++++++++" << std::endl;
diff --git a/tests/ippl_src/Field/FieldDebug.cpp b/tests/ippl_src/Field/FieldDebug.cpp
index aaf18392a..782341507 100644
--- a/tests/ippl_src/Field/FieldDebug.cpp
+++ b/tests/ippl_src/Field/FieldDebug.cpp
@@ -11,8 +11,8 @@
#include <fstream>
namespace {
- void hardCodedOutput(char* filename); // Prototype of function defined below.
- bool thediff(char* filename1, char* filename2);
+ void hardCodedOutput(std::string filename); // Prototype of function defined below.
+ bool thediff(std::string filename1, std::string filename2);
}
constexpr double roundOffError = 1e-10;
@@ -99,7 +99,7 @@ namespace {
//-----------------------------------------------------------------------------
// Mock up the Unix "diff" utility to compare two files:
//-----------------------------------------------------------------------------
- bool thediff(char* filename1, char* filename2)
+ bool thediff(std::string filename1, std::string filename2)
{
bool same = true;
char ch1, ch2;
@@ -116,7 +116,7 @@ namespace {
}
//-----------------------------------------------------------------------------
- void hardCodedOutput(char* filename)
+ void hardCodedOutput(std::string filename)
{
std::ofstream of(filename);
of << std::endl
diff --git a/tests/ippl_src/Index/Index.cpp b/tests/ippl_src/Index/Index.cpp
index 1ebaf07dd..6f66973b9 100644
--- a/tests/ippl_src/Index/Index.cpp
+++ b/tests/ippl_src/Index/Index.cpp
@@ -1,7 +1,11 @@
#include "gtest/gtest.h"
+#include "opal_test_utilities/SilenceTest.h"
+
#include "Index/SOffset.h"
+#include <iostream>
+
constexpr unsigned Dim = 2;
constexpr double roundOffError = 1e-10;
@@ -55,4 +59,346 @@ TEST(Index, SOffset)
NDIndex<Dim> N1(Index(0,1), Index(0,1));
EXPECT_TRUE(! A.inside(N1));
EXPECT_TRUE( B.inside(N1));
+}
+
+#include "Index/SIndex.h"
+
+TEST(Index, SIndex)
+{
+ OpalTestUtilities::SilenceTest silencer;
+
+ Index I(25);
+ Index J(50);
+ Index I2(2);
+ Index J2(5);
+ NDIndex<Dim> NDX(I,J);
+ int IP[Dim];
+ IP[0] = -1;
+ IP[1] = -1;
+
+ SOffset<Dim> so1(1,1);
+ SOffset<Dim> so2(IP);
+ std::cout << "Created SOffset so1 = " << so1 << std::endl;
+ std::cout << "Created SOffset so2 = " << so2 << std::endl;
+ std::cout << "Adding IP to so1 = " << so1 + IP << std::endl;
+ EXPECT_EQ((so1+IP)[0],0);
+ std::cout << "Adding so1 to IP = " << IP + so1 << std::endl;
+ EXPECT_EQ((IP+so1)[1],0);
+ std::cout << "Adding so2 to so1 = " << so1 + so2 << std::endl;
+ EXPECT_EQ((so1+so2)[0],0);
+ std::cout << "Adding so1 to NDX " << NDX << " = " << NDX + so1 << std::endl;
+ std::cout << "Adding NDX " << NDX << " to so1 = " << so1 + NDX << std::endl;
+ std::cout << "Multiplying NDX " << NDX << " by so2 = " << NDX * so2 << std::endl;
+ std::cout << "Multiplying so2 by NDX " << NDX << " = " << so2 * NDX << std::endl;
+ so1 -= IP;
+ so2 += IP;
+ std::cout << "Accumulated IP from so1 = " << so1 << std::endl;
+ EXPECT_EQ((so1)[0],2);
+ std::cout << "Accumulated IP into so2 = " << so2 << std::endl;
+ EXPECT_EQ((so2)[0],-2);
+
+ FieldLayout<Dim> layout(I, J, PARALLEL, PARALLEL, 2*Ippl::getNodes());
+ Field<double,Dim> A(layout);
+ Field<bool,Dim> B(layout);
+
+ SIndex<Dim> s1(layout);
+ SIndex<Dim> s2 = s1(1,-1);
+ SIndex<Dim> x3 = s1(IP);
+
+ std::cout << "Created s1 = " << s1 << std::endl;
+ std::cout << "Created s2 = " << s2 << std::endl;
+ std::cout << "Created x3 = " << x3 << std::endl;
+
+ s1.addIndex(NDIndex<Dim>(Index(2), Index(3)));
+ s2.addIndex(SOffset<Dim>(0,0));
+ s2.addIndex(NDIndex<Dim>(Index(20,23), Index(45,46)));
+
+ // testmsg << "Added new points, s1 = "; s1.printDebug(testmsg);
+ // testmsg << "Added new points, s2 = "; s2.printDebug(testmsg);
+
+ SIndex<Dim> s3(s1);
+ s3 = NDIndex<Dim>(Index(1,5), Index(2,4));
+
+ std::cout << "Created s3 = " << s3 << std::endl;
+
+ s3 &= s1;
+
+ std::cout << "Intersection of s3 and s1 = " << s3 << std::endl;
+
+ // now, test assigment of a Field expression to an SIndex
+ A[I][J] = I + J;
+ B = lt(A,10);
+ s3 = lt(A,10);
+
+ // do a union with a slightly different condition
+ s3 |= (gt(A,15) && lt(A,20));
+ std::cout << "union of s3 and where 15 < A < 20 = " << s3 << std::endl;
+
+ // do an intersection with an overlapping condition
+ s3 &= (gt(A,8) && lt(A,12));
+ std::cout << "intersection of s3 and where 8 < A < 12 = " << s3 << std::endl;
+
+ // do an indexed assignment to s3
+ s3[I2][J2] = (lt(A[I2 + 1][J2+1], 5) && gt(A[I2 + 2][J2 + 2], 0));
+ std::cout << "s3[I2][J2] = expr ==> s3 = " << s3 << std::endl;
+
+ // now use s3 in a Field expression
+ std::cout << "Originally, A = I + J ... after A[s3] = A[s3(1,1)]: A = ";
+ A[s3] = A[s3(1,1)] + 10;
+ std::cout << A << std::endl;
+
+ EXPECT_EQ(sum(A), 45685);
+}
+
+TEST(Index, SubField)
+{
+ OpalTestUtilities::SilenceTest silencer;
+
+ Index I(4);
+ Index I2(2);
+ Index J(4);
+ Index K(2);
+ FieldLayout<Dim> layout(I, J, PARALLEL, PARALLEL, 4*Ippl::getNodes());
+ FieldLayout<3> layout3(I, J, K, PARALLEL, PARALLEL, PARALLEL,
+ 2*Ippl::getNodes());
+ BConds<double,Dim> dbc;
+ for (unsigned int f=0; f < 2*Dim; f++) dbc[f] = new ZeroFace<double,Dim>(f);
+ Field<double,Dim> A(layout, dbc, GuardCellSizes<Dim>(1));
+ A = 0.0;
+ BConds<bool,Dim> bbc;
+ for (unsigned int f=0; f < 2*Dim; f++) bbc[f] = new ZeroFace<bool,Dim>(f);
+ Field<bool,Dim> B(layout, bbc, GuardCellSizes<Dim>(1));
+ B = true;
+ Field<double,Dim> C(layout, dbc, GuardCellSizes<Dim>(1));
+ C = 0.0;
+
+ BConds<float,3> fbc;
+ for (unsigned int f=0; f < 2*3; f++) fbc[f] = new ZeroFace<float,3>(f);
+ Field<float,3> A3(layout3, fbc, GuardCellSizes<3>(2));
+ A3 = 0.0;
+ Field<float,3> B3(layout3, fbc, GuardCellSizes<3>(2));
+ B3 = 0.0;
+ SIndex<Dim> s1(layout);
+ s1.addIndex(NDIndex<Dim>(Index(1,1), Index(1,2)));
+
+ SIndex<Dim> s2 = s1(1,-1);
+ s2.addIndex(SOffset<Dim>(0,0));
+
+ std::cout << "\n************ testing 2D SubParticleAttrib<> ************" << std::endl;
+ A[I][J] = (I+1) + (J*10);
+ std::cout << "s1 = " << s1 << std::endl;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),280);
+ // not working anymore
+ // ParticleAttrib<double> PA;
+ // PA[s1] = A[s1];
+ // std::cout << "Result of PA[s1] = A[s1] : PA = " << std::endl;
+ // std::cout << PA[s1] << std::endl;
+
+ // std::cout << "s2 = " << s2 << std::endl;
+ // PA[s2] = A[s2] + A[s2];
+ // std::cout << "Result of PA[s2] = A[s2] + A[s2] : PA = " << std::endl;
+ // std::cout << PA[s2] << std::endl;
+
+ // PA[s2] *= (A[s2] + PA[s2]);
+ // std::cout << "Result of PA[s2] *= A[s2] + PA[s2] : PA = " << std::endl;
+ // std::cout << PA[s2] << std::endl;
+
+ // A = 1.0;
+ SIndex<Dim> sX = s1;
+ sX.addIndex(SOffset<Dim>(0,0));
+ std::cout << "Created sX = " << sX << std::endl;
+
+ // A[sX] += PA[sX];
+ // std::cout << "Result of A[sX] = 1 + PA[sX] : A = " << A << std::endl;
+
+ A = 0.0;
+ C = 1.0;
+ A[sX] = C[sX];
+ std::cout << "Result of A[sX] = C[sX], after C = 1.0: A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),3);
+
+ // ParticleAttrib<Vektor<double,2> > PA2;
+ // PA2(0)[sX] = A[sX] + PA[sX];
+ // PA2(1)[sX] = -(A[sX] + PA[sX]);
+ // std::cout << "Result of PA2[sX] = +- (A[sX] + PA[sX]):" << std::endl;
+ // std::cout << " A = " << A << std::endl;
+ // std::cout << " PA = " << PA[sX] << std::endl;
+ // std::cout << "PA2 = " << PA2[sX] << std::endl;
+
+ std::cout << "\n************ testing 2D SubBareField<NDIndex> ************" << std::endl;
+
+ Index IX(0,2);
+ Index JX(2,2);
+ NDIndex<Dim> sub1(IX, JX);
+ SubBareField<double,Dim,NDIndex<Dim> > SA1(A, sub1);
+ std::cout << "Created SubBareField SA1 from A and NDIndex " << sub1 << ":" <<std::endl;
+ std::cout << " SA1 = " << SA1 << std::endl;
+
+ SA1 = 3.0;
+ std::cout << "SA1 set to 3:" << std::endl;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),11);
+
+ NDIndex<Dim> sub2(IX + 1, JX);
+ SubBareField<double,Dim,NDIndex<Dim> > SA2(A, sub2);
+ std::cout << "Created SubBareField SA2 from A and NDIndex " << sub2 << ":" <<std::endl;
+ std::cout << " SA2 = " << SA2 << std::endl;
+
+ SA1 = SA2 + 4.0;
+ std::cout << "SA1 on " << sub1 << " set to (SA2 + 4):" << std::endl;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),20);
+
+ std::cout << "\n************ testing 2D SubBareField<SIndex> ************" << std::endl;
+
+ SubBareField<double,Dim,SIndex<Dim> > SB1(A, s1);
+ std::cout << "Created SubBareField SB1 from A and SIndex " << s1 << ":" << std::endl;
+ std::cout << " SB1 = " << SB1 << std::endl;
+
+ SB1 = -2.0;
+ std::cout << "SB1 set to -2:" << std::endl;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),5);
+
+ SubBareField<double,Dim,SIndex<Dim> > SB2(A, s1(1,-1));
+ std::cout << "Created SubBareField SB2 from A and SIndex " << SB2.getDomain();
+ std::cout << ":" << std::endl;
+ std::cout << " SB2 = " << SB2 << std::endl;
+
+ SB1 = SB2 + 4.0;
+ std::cout << "SB1 set to (SB2 + 4):" << std::endl;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),23);
+
+ std::cout << "\n************ testing 2D SubBareField<SOffset> ************" << std::endl;
+
+ SubBareField<double,Dim,SOffset<Dim> > SC1(A, SOffset<Dim>(2,1));
+ std::cout << "Created SubBareField SC1 from A and SOffset " << SC1.getDomain();
+ std::cout << ":" << std::endl;
+ std::cout << " SC1 = " << SC1 << std::endl;
+
+ SC1 = 10.0;
+ std::cout << "SC1 set to 10:" << std::endl;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),33);
+
+ SubBareField<double,Dim,SOffset<Dim> > SC2(A, SOffset<Dim>(2,2));
+ SubBareField<double,Dim,SOffset<Dim> > SC3(A, SOffset<Dim>(3,2));
+ SC1 = SC1 + (SC2 + SC3) * 10;
+ std::cout << "SC1 set to 10 + [(2,2) + (3,2)]*10:" << std::endl;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),73);
+
+ std::cout << "\n************ testing 2D Field[SIndex] ************" << std::endl;
+
+ std::cout << "Current SIndex s1 = " << s1 << std::endl;
+ std::cout << "Results of A[s1] = -1000:" << std::endl;
+ A[s1] = -1000.0;
+ std::cout << " A = " << A << " " << sum(A) << std::endl;
+ EXPECT_EQ(sum(A),-2939);
+ B[s1] = ne(A[s1], A[s1]);
+ std::cout << "Results of B[s1] = (A[s1] != A[s1]):" << std::endl;
+ std::cout << " B = " << B << " " << sum(B) << std::endl;
+ EXPECT_EQ(sum(B),1);
+ C[s1] = 1.0;
+ std::cout << "Results of C[s1] = 1:" << std::endl;
+ std::cout << " C = " << C << " " << sum(C) << std::endl;
+ EXPECT_EQ(sum(C),16);
+ C[s1] = A[s1(-1,-1)] - 2000;
+ std::cout << "Results of C[s1] = (A[s1(-1,-1)] - 2000) :" << std::endl;
+ std::cout << " C = " << C << " " << sum(C) << std::endl;
+ EXPECT_EQ(sum(C),-6987);
+
+ std::cout << "\n************ testing 2D Field[SIndex,compressed] ****" << std::endl;
+
+ Field<double,Dim> A2(layout);
+ A2 = 1.0;
+ std::cout << " A2 = " << A2 << " " << sum(A2) << std::endl;
+ EXPECT_EQ(sum(A2),16);
+ std::cout << " A2.compressedFraction = " << A2.CompressedFraction() << std::endl;
+ A2[I2][I2] = -10.0;
+ std::cout << "Initial settings:" << std::endl;
+ std::cout << " A2 = " << A2 << " " << sum(A2) << std::endl;
+ EXPECT_EQ(sum(A2),-28);
+ std::cout << " A2.compressedFraction = " << A2.CompressedFraction() << std::endl;
+ s1 = ne(A2, 0.0);
+ std::cout << " s1 = " << s1 << std::endl;
+ std::cout << " A2 = " << A2 << " " << sum(A2) << std::endl;
+ EXPECT_EQ(sum(A2),-28);
+ std::cout << " A2.compressedFraction = " << A2.CompressedFraction() << std::endl;
+ A2[s1] = A2[s1] + 100.0;
+ std::cout << "Results of A2[s1] = A2[s1] + 100.0 :" << std::endl;
+ std::cout << " A2 = " << A2 << " " << sum(A2) << std::endl;
+ EXPECT_EQ(sum(A2),1572);
+ std::cout << " A2.compressedFraction = " << A2.CompressedFraction() << std::endl;
+ std::cout << " s1 = " << s1 << std::endl;
+
+ std::cout << "\n************ testing 3D Field[SIndex] ************" << std::endl;
+
+ A3[I][J][K] = I + J + K;
+ B3 = 0.0;
+
+ SIndex<3> sindx3(layout3);
+ sindx3 = eq(A3, 2.0);
+
+ SOffset<3> offset;
+ offset[0] = 1;
+ offset[1] = 0;
+ offset[2] = (-1);
+
+ std::cout << "Current 3D SIndex sindx3 = " << sindx3 << std::endl;
+ std::cout << "Current 3D A3 = " << A3 << " " << sum(A3) << std::endl;
+ EXPECT_EQ(sum(A3),112);
+ std::cout << "A3 after A3[sindx3] = -1:" << std::endl;
+ A3[sindx3] = -1.0;
+ std::cout << " A3 = " << A3 << " " << sum(A3) << std::endl;
+ EXPECT_EQ(sum(A3),97);
+ B3[sindx3] = 1.0;
+ std::cout << "Results of B3[sindx3] = 1:" << std::endl;
+ std::cout << " B3 = " << B3 << " " << sum(B3) << std::endl;
+ EXPECT_EQ(sum(B3),5);
+ B3[sindx3] = A3[sindx3(offset)] - 2000.0;
+ std::cout << "Results of B3[sindx3] = A3[sindx3(1,0,-1)] - 2000.0:" << std::endl;
+ std::cout << " B3 = " << B3 << " " << sum(B3) << std::endl;
+ EXPECT_EQ(sum(B3),-10002);
+}
+
+#include "Field/Field.h"
+#include "FieldLayout/FieldLayout.h"
+
+/***************************************************************************
+ A simple program to test SubField assignments using two fields, one on
+ a cell-centered layout, the other on a vertex-centered (but otherwise
+ similar) layout.
+***************************************************************************/
+
+TEST(Index, VertCell)
+{
+ OpalTestUtilities::SilenceTest silencer;
+
+
+ Index I(4), J(4), U(5), V(5);
+ FieldLayout<Dim> LC(I, J, PARALLEL, PARALLEL, 4);
+ FieldLayout<Dim> LV(U, V, PARALLEL, PARALLEL, 4);
+
+ Field<int,Dim,UniformCartesian<Dim>,Cell> A(LC, GuardCellSizes<Dim>(1));
+ Field<int,Dim,UniformCartesian<Dim>,Vert> B(LV, GuardCellSizes<Dim>(1));
+
+ A[I][J] = J + 1;
+ B[U][V] = -V - 1;
+ std::cout << "A at start: " << A << std::endl;
+ std::cout << "B at start: " << B << std::endl;
+
+ SIndex<Dim> SI(LC);
+ SI = eq(A,2) || eq(A,4);
+ std::cout << "SI at start: " << SI << std::endl;
+
+ A[SI] = A[SI(0,-1)] + B[SI(1,1)];
+
+ std::cout << "A at end: " << A << std::endl;
+
+ int s = sum(A);
+
+ EXPECT_EQ(s, 0);
}
\ No newline at end of file
diff --git a/tests/ippl_src/Meshes/CartesianCentering.cpp b/tests/ippl_src/Meshes/CartesianCentering.cpp
index 85adb3570..33cedbe5f 100644
--- a/tests/ippl_src/Meshes/CartesianCentering.cpp
+++ b/tests/ippl_src/Meshes/CartesianCentering.cpp
@@ -12,8 +12,8 @@
// forward declarations
namespace {
- void hardCodedOutput(char* filename); // Prototype of function defined below.
- bool thediff(char* filename1, char* filename2);
+ void hardCodedOutput(std::string filename); // Prototype of function defined below.
+ bool thediff(std::string filename1, std::string filename2);
extern const CenteringEnum zz[2] = {CELL, VERTEX};
}
@@ -84,7 +84,7 @@ namespace {
//-----------------------------------------------------------------------------
// Mock up the Unix "diff" utility to compare two files:
//-----------------------------------------------------------------------------
-bool thediff(char* filename1, char* filename2)
+bool thediff(std::string filename1, std::string filename2)
{
bool same = true;
char ch1, ch2;
@@ -102,7 +102,7 @@ bool thediff(char* filename1, char* filename2)
}
//-----------------------------------------------------------------------------
-void hardCodedOutput(char* filename)
+void hardCodedOutput(std::string filename)
{
std::ofstream of(filename);
of << "CartesianCentering: no specialized name (yet) for this case" << std::endl
--
GitLab