APPSPACK_Directions.hpp

Go to the documentation of this file.

// $Source: /space/CVS-Acro/acro/packages/appspack/appspack/src/APPSPACK_Directions.hpp,v $ 

//@HEADER
// ************************************************************************
// 
//          APPSPACK: Asynchronous Parallel Pattern Search
//                 Copyright (2003) Sandia Corporation
// 
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
// license for use of this work by or on behalf of the U.S. Government.
// 
// This library is free software; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2.1 of the
// License, or (at your option) any later version.
//  
// This library is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//                                                                                 
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA.                                                                           .
// 
// Questions? Contact Tammy Kolda (tgkolda@sandia.gov) 
// 
// ************************************************************************
//@HEADER

#ifndef APPSPACK_DIRECTIONS_HPP
#define APPSPACK_DIRECTIONS_HPP

#include "APPSPACK_Point.hpp"
#include "APPSPACK_Parameter_List.hpp"
#include "APPSPACK_Constraints_Linear.hpp"

namespace APPSPACK
{


class Directions
{

public:


  Directions(Parameter::List& params, const Constraints::Linear& constraints_in);

  ~Directions();
  


  void  getDirectionIndices(vector<int>& idvector) const;

  const Vector& getDirection(int i) const;

  double getStep(int i) const;

  void setStepConverged(int i);


  bool isStepConverged() const;

  bool empty() const;



  void computeNewDirections(const Point& newPoint);


  void appendNewDirections();

  void setTrueStepAndTag(int i, double trueStep_in, int tag_in);


  void reduceStep(int i);



  void print(const string label) const;

  double getSmallestStep() const;

private:


  void buildNormalCone(Matrix& VpT, Matrix& VlT) const;


  void buildTangentCone(const Matrix& VpT, const Matrix& VlT,
                          Matrix& T);
         

  void buildWithNothing(Matrix& D);
                       

  bool buildWithCddlib(const Matrix& VpT, const Matrix& VlT,
                          Matrix& T);


  bool buildWithLapack(const Matrix& VpT, const Matrix& VlT,
                          Matrix& T);

  

  void generateUnconstrained(Matrix& D);
 
  void generateForLinear(Matrix& D);
                         

  void addNormalDirections(const Matrix& VpT, const Matrix& VlT,
                           Matrix& D);


  void addCompassDirections(const Matrix& VlT, Matrix& D);
                           

  bool updateConstraintState(double newStep);
  

  void updateDirectionInfo(double newStep, bool isAppend=false);

private:
  
  const Constraints::Linear& constraints;

  const int nDimensions;

  const Vector zero;

  const double stepTolerance;

  const double minStep;

  const double theta;

  const double epsilonMax;

  int nDirections;

  Matrix direction;

  Vector step;


  Vector trueStep;


  vector<int> tag;

  Vector tmpVector;

  mutable vector<int> idxVector;


  map< vector<Constraints::ActiveType>, Matrix>  directionCache;

  map< vector<Constraints::ActiveType>, Matrix>::iterator  cacheIter;
  

  vector<Constraints::ActiveType> constraintState;

  double epsilonMin;
  
  bool withNormals;  

  bool withCompass;  

  int nCached;
  
  int nLapack;
  
  int nCddlib;

  int nMaxDirections;

  int nAppend;

  Vector xDistance;

};

}

#endif

 

© Sandia Corporation | Site Contact | Privacy and Security

Generated on Fri Feb 16 10:33:35 2007 for APPSPACK 5.0.1 by 1.3.9.1 written by Dimitri van Heesch, © 1997-2002