Changeset e00f47

Timestamp:

Apr 21, 2008, 2:19:24 PM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

961b75

Parents:

64fa9e

git-author:

Frederik Heber <heber@…> (04/18/08 15:14:39)

git-committer:

Frederik Heber <heber@…> (04/21/08 14:19:24)

Message:

GramSchNormalize() added

Location:

pcp/src

Files:

• gramsch.c (modified) (2 diffs)
• gramsch.h (modified) (1 diff)

pcp/src/gramsch.c ¶

@@ -168 +168 @@
  * \return Squared norm of wave function
  */
-static double GramSchNormalize(const struct Problem *P, struct LatticeLevel *Lev, fftw_complex *LPsiDat, double PsiSP) {
+double GramSchNormalize(const struct Problem *P, struct LatticeLevel *Lev, fftw_complex *LPsiDat, double PsiSP) {
   double LocalSP=0.0;
   int i,s = 0;
@@ -182 +182 @@
     MPI_Allreduce ( &LocalSP, &PsiSP, 1, MPI_DOUBLE, MPI_SUM, P->Par.comm_ST_Psi);
   } 
-  if (PsiSP < MYEPSILON) fprintf(stderr,"GramSchNormalize: PsiSP = %lg\n",PsiSP);
+  if ((PsiSP < MYEPSILON) && (P->Call.out[PsiOut])) fprintf(stderr,"GramSchNormalize: PsiSP = %lg\n",PsiSP);
   PsiSP = sqrt(PsiSP); // take square root
   for (i=0; i < Lev->MaxG; i++) { // and divide each coefficient by the norm

pcp/src/gramsch.h ¶

@@ -19 +19 @@
 void FreeMPI_OnePsiElement();
 void FirstInitGramSchData(struct Problem *P, struct Psis *Psi);
+void GramSch(struct Problem *P, struct LatticeLevel *Lev, struct Psis *Psi, enum PsiGramSchToDoType ToDo);
+double GramSchNormalize(const struct Problem *P, struct LatticeLevel *Lev, fftw_complex *LPsiDat, double PsiSP);
 double GramSchGetNorm2(const struct Problem *P, struct LatticeLevel *Lev, fftw_complex *LPsiDat);
-void GramSch(struct Problem *P, struct LatticeLevel *Lev, struct Psis *Psi, enum PsiGramSchToDoType ToDo);
 void SetGramSchExtraPsi(const struct Problem *P, struct Psis *Psi, enum PsiGramSchStatusType PsiGramSchStatus);
 void SetGramSchActualPsi(const struct Problem *P, struct Psis *Psi, enum PsiGramSchStatusType PsiGramSchStatus);