Albany_SolutionMaxValueResponseFunction.hpp

Go to the documentation of this file.

//*****************************************************************//
//    Albany 2.0:  Copyright 2012 Sandia Corporation               //
//    This Software is released under the BSD license detailed     //
//    in the file "license.txt" in the top-level Albany directory  //
//*****************************************************************//

#ifndef ALBANY_SOLUTIONMAXVALUERESPONSEFUNCTION_HPP
#define ALBANY_SOLUTIONMAXVALUERESPONSEFUNCTION_HPP

#include "Albany_SamplingBasedScalarResponseFunction.hpp"

namespace Albany {

  class SolutionMaxValueResponseFunction : 
    public SamplingBasedScalarResponseFunction {
  public:
  
    SolutionMaxValueResponseFunction(
      const Teuchos::RCP<const Epetra_Comm>& comm, 
      int neq = 1, int eq = 0, bool interleavedOrdering=true);

    virtual ~SolutionMaxValueResponseFunction();

    virtual unsigned int numResponses() const;

    virtual void 
    evaluateResponse(const double current_time,
         const Epetra_Vector* xdot,
         const Epetra_Vector* xdotdot,
         const Epetra_Vector& x,
         const Teuchos::Array<ParamVec>& p,
         Epetra_Vector& g);

    virtual void 
    evaluateTangent(const double alpha, 
        const double beta,
        const double omega,
        const double current_time,
        bool sum_derivs,
        const Epetra_Vector* xdot,
        const Epetra_Vector* xdotdot,
        const Epetra_Vector& x,
        const Teuchos::Array<ParamVec>& p,
        ParamVec* deriv_p,
        const Epetra_MultiVector* Vxdot,
        const Epetra_MultiVector* Vxdotdot,
        const Epetra_MultiVector* Vx,
        const Epetra_MultiVector* Vp,
        Epetra_Vector* g,
        Epetra_MultiVector* gx,
        Epetra_MultiVector* gp);

    virtual void 
    evaluateGradient(const double current_time,
         const Epetra_Vector* xdot,
         const Epetra_Vector* xdotdot,
         const Epetra_Vector& x,
         const Teuchos::Array<ParamVec>& p,
         ParamVec* deriv_p,
         Epetra_Vector* g,
         Epetra_MultiVector* dg_dx,
         Epetra_MultiVector* dg_dxdot,
         Epetra_MultiVector* dg_dxdotdot,
         Epetra_MultiVector* dg_dp);

  private:

    SolutionMaxValueResponseFunction(const SolutionMaxValueResponseFunction&);
    
    SolutionMaxValueResponseFunction& operator=(const SolutionMaxValueResponseFunction&);

  protected:

    int neq;

    int eq;

    bool interleavedOrdering;

    void computeMaxValue(const Epetra_Vector& x, double& val, int& index);

  };

}

#endif // ALBANY_SOLUTIONMAXVALUERESPONSEFUNCTION_HPP