PrintHamiltonianAverage again

TestSuite/inputs/input1801.inp

@@ -18,6 +18,7 @@ LadderLeg=2
 Model=Heisenberg
 HeisenbergTwiceS=1
 
+PrintHamiltonianAverage=s==c
 SolverOptions=twositedmrg,restart,TargetingAncilla
 Version=version
 OutputFile=data1801.txt

TestSuite/inputs/input1811.inp

@@ -19,6 +19,7 @@ Connectors 1 0.0
 Model=HeisenbergAncillaC
 HeisenbergTwiceS=1
 
+PrintHamiltonianAverage=s==c
 SolverOptions=twositedmrg,restart,TargetingAncilla
 Version=version
 OutputFile=data1811.txt

TestSuite/inputs/input1812.inp

@@ -19,6 +19,7 @@ Connectors 1 0.0
 Model=HeisenbergAncillaC
 HeisenbergTwiceS=1
 
+PrintHamiltonianAverage=s==c
 SolverOptions=twositedmrg,restart,TargetingAncilla,neverNormalizeVectors
 Version=version
 OutputFile=data1812

TestSuite/inputs/input1851.inp

@@ -30,6 +30,7 @@ Model=TjAncillaC
 potentialV	 12 0 0 0 0 0 0
                     0 0 0 0 0 0
 
+PrintHamiltonianAverage=s==c
 SolverOptions=restart,TargetingAncilla,twositedmrg
 Version=version
 OutputFile=data1851.txt

TestSuite/inputs/input1861.inp

@@ -30,6 +30,7 @@ Model=TjAncillaC2
 potentialV	 12 0 0 0 0 0 0 
                     0 0 0 0 0 0 
 
+PrintHamiltonianAverage=s==c
 SolverOptions=restart,TargetingAncilla,twositedmrg
 Version=version
 OutputFile=data1861.txt

TestSuite/inputs/input1951.inp

@@ -28,6 +28,7 @@ FiniteLoops 4
 RepeatFiniteLoopsFrom=0
 RepeatFiniteLoopsTimes=10
 
+PrintHamiltonianAverage=s==c
 TargetElectronsTotal=8
 TargetSzPlusConst=4
 TargetExtra 2 2 2

TestSuite/inputs/input2101.inp

@@ -36,6 +36,7 @@ FiniteLoops 2
 
 RepeatFiniteLoopsTimes=10
 RecoverySave=l%2,keep,M=100
+PrintHamiltonianAverage=s==c
 
 TruncationTolerance=1e-7
 TridiagEps=1e-7

TestSuite/inputs/input2102.inp

@@ -30,7 +30,7 @@ InfiniteLoopKeptStates=data2101
 FiniteLoops 2
 10 400 2  -10 400 2
 
-
+PrintHamiltonianAverage=s==c
 SolverOptions=CorrectionVectorTargeting,restart,twositedmrg
 CorrectionA=0
 Version=version

src/Engine/TargetingMetts.h

@@ -278,8 +278,8 @@ public:
 		this->common().cocoon(block1, direction, doBorderIfBorder);
 
 		PsimagLite::String predicate = model_.params().printHamiltonianAverage;
-		const SizeType linSize = model_.geometry().numberOfSites();
-		PsimagLite::PredicateAwesome<>::replaceAll(predicate, "c", ttos(linSize));
+		const SizeType center = model_.geometry().numberOfSites()/2;
+		PsimagLite::PredicateAwesome<>::replaceAll(predicate, "c", ttos(center));
 		PsimagLite::PredicateAwesome<> pAwesome(predicate);
 		assert(block1.size() > 0);
 		if (pAwesome.isTrue("s", block1[0]))

src/Engine/TargetingTimeStep.h

@@ -275,8 +275,8 @@ private:
 		this->common().cocoon(block1, direction, doBorderIfBorder);
 
 		PsimagLite::String predicate = this->model().params().printHamiltonianAverage;
-		const SizeType linSize = this->model().geometry().numberOfSites();
-		PsimagLite::PredicateAwesome<>::replaceAll(predicate, "c", ttos(linSize));
+		const SizeType center = this->model().geometry().numberOfSites()/2;
+		PsimagLite::PredicateAwesome<>::replaceAll(predicate, "c", ttos(center));
 		PsimagLite::PredicateAwesome<> pAwesome(predicate);
 		assert(block1.size() > 0);
 		if (pAwesome.isTrue("s", block1[0]))