Albany_SolutionTwoNormResponseFunction.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_SOLUTIONTWONORMRESPONSEFUNCTION_HPP
#define ALBANY_SOLUTIONTWONORMRESPONSEFUNCTION_HPP

#include "Albany_SamplingBasedScalarResponseFunction.hpp"

namespace Albany {

  class SolutionTwoNormResponseFunction : 
    public SamplingBasedScalarResponseFunction {
  public:
  
    SolutionTwoNormResponseFunction(
      const Teuchos::RCP<const Epetra_Comm>& comm);

    virtual ~SolutionTwoNormResponseFunction();

    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:

    SolutionTwoNormResponseFunction(const SolutionTwoNormResponseFunction&);
    
    SolutionTwoNormResponseFunction& operator=(const SolutionTwoNormResponseFunction&);

  };

}

#endif // ALBANY_SOLUTIONTWONORMRESPONSEFUNCTION_HPP