src/minimizers/CGFamilyMinimizer.cpp

#include <cmath>

#include "minimizers/OOPSEMinimizer.hpp"
#include "utils/Utility.hpp"

CGFamilyMinmizer::CGFamilyMinmizer(SimInfo *theInfo
                                                               MinimizerParameterSet * param)
                         :OOPSEMinimizer(theInfo, param){
  prevG.resize(ndim);
  prevX.resize(ndim);
  direction.resize(ndim);

  stepSize = paramSet->getStepSize();

}
int CGFamilyMinmizer::checkConvg(){
  double fTol;
  double relativeFTol;  // relative tolerance
  double deltaF;
  double gTol;
  double relativeGTol;
  double gnorm;
  

  // test function tolerance test
  fTol =paramSet->getFTol();
  relativeFTol = fTol * std::max(1.0,fabs(curF));  // relative tolerance
  deltaF = prevF - curF;
  
  if (fabs(deltaF) <= relativeFTol) {

    if (bVerbose){
      cout << "function value tolerance test passed" << std::endl;
      cout << "ftol = " << fTol
             << "\tdeltaf = " << deltaF<< std::endl;
    }
    return CONVG_FTOL;
  }
  
//gradient tolerance test
  gTol = paramSet->getGTol();
  relativeGTol = gTol * std::max(1.0,fabs(curF));

#ifndef IS_MPI
  gnorm = sqrt(dotProduct(curG, curG));
#else
  double localDP;
  double globalDP;

  localDP = dotProduct(curG, curG);
  MPI_Allreduce(&localDP, &globalDP, 1, MPI_DOUBLE,MPI_SUM, MPI_COMM_WORLD);  
  gnorm  = sqrt(globalDP);
#endif

  if (gnorm <= relativeGTol) {
      cout << "gradient tolerance test" << std::endl;
      cout << "gnorm = " << gnorm
             << "\trelativeGTol = " << relativeGTol<< std::endl;
    return CONVG_GTOL;
  }
  
  //absolute gradient tolerance test

  if (gnorm <= gTol) {
      cout << "absolute gradient tolerance test" << std::endl;
      cout << "gnorm = " << gnorm
             << "\tgTol = " << gTol<< std::endl;
    return CONVG_ABSGTOL;
  }

  // did not converge yet
  return CONVG_UNCONVG;
}

Revision:	1884
Committed:	Tue Dec 14 19:08:44 2004 UTC (20 years, 10 months ago) by tim
File size:	1848 byte(s)
Log Message:	more fix in MPI version
#	User	Rev	Content
1	tim	1884	#include <cmath>
2	gezelter	1490
3	tim	1492	#include "minimizers/OOPSEMinimizer.hpp"
4			#include "utils/Utility.hpp"
5	gezelter	1490
6	tim	1884	CGFamilyMinmizer::CGFamilyMinmizer(SimInfo *theInfo
7	gezelter	1490	MinimizerParameterSet * param)
8	tim	1884	:OOPSEMinimizer(theInfo, param){
9	gezelter	1490	prevG.resize(ndim);
10			prevX.resize(ndim);
11			direction.resize(ndim);
12
13			stepSize = paramSet->getStepSize();
14
15			}
16			int CGFamilyMinmizer::checkConvg(){
17			double fTol;
18			double relativeFTol; // relative tolerance
19			double deltaF;
20			double gTol;
21			double relativeGTol;
22			double gnorm;
23
24
25			// test function tolerance test
26			fTol =paramSet->getFTol();
27			relativeFTol = fTol * std::max(1.0,fabs(curF)); // relative tolerance
28			deltaF = prevF - curF;
29
30			if (fabs(deltaF) <= relativeFTol) {
31
32			if (bVerbose){
33	tim	1830	cout << "function value tolerance test passed" << std::endl;
34	gezelter	1490	cout << "ftol = " << fTol
35	tim	1830	<< "\tdeltaf = " << deltaF<< std::endl;
36	gezelter	1490	}
37			return CONVG_FTOL;
38			}
39
40			//gradient tolerance test
41			gTol = paramSet->getGTol();
42			relativeGTol = gTol * std::max(1.0,fabs(curF));
43
44			#ifndef IS_MPI
45			gnorm = sqrt(dotProduct(curG, curG));
46			#else
47			double localDP;
48			double globalDP;
49
50			localDP = dotProduct(curG, curG);
51			MPI_Allreduce(&localDP, &globalDP, 1, MPI_DOUBLE,MPI_SUM, MPI_COMM_WORLD);
52			gnorm = sqrt(globalDP);
53			#endif
54
55			if (gnorm <= relativeGTol) {
56	tim	1830	cout << "gradient tolerance test" << std::endl;
57	gezelter	1490	cout << "gnorm = " << gnorm
58	tim	1830	<< "\trelativeGTol = " << relativeGTol<< std::endl;
59	gezelter	1490	return CONVG_GTOL;
60			}
61
62			//absolute gradient tolerance test
63
64			if (gnorm <= gTol) {
65	tim	1830	cout << "absolute gradient tolerance test" << std::endl;
66	gezelter	1490	cout << "gnorm = " << gnorm
67	tim	1830	<< "\tgTol = " << gTol<< std::endl;
68	gezelter	1490	return CONVG_ABSGTOL;
69			}
70
71			// did not converge yet
72			return CONVG_UNCONVG;
73			}
74