Albany_AggregateScalarResponseFunction.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_AGGREGATE_SCALAR_RESPONSE_FUNCTION_HPP
#define ALBANY_AGGREGATE_SCALAR_RESPONSE_FUNCTION_HPP

#include "Albany_SamplingBasedScalarResponseFunction.hpp"
#include "Teuchos_Array.hpp"

namespace Albany {

  class AggregateScalarResponseFunction : 
    public SamplingBasedScalarResponseFunction {
  public:
  
    AggregateScalarResponseFunction(
      const Teuchos::RCP<const Epetra_Comm>& comm,
      const Teuchos::Array< Teuchos::RCP<ScalarResponseFunction> >& responses);

    virtual void setup();

    virtual ~AggregateScalarResponseFunction();

    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:

    AggregateScalarResponseFunction(const AggregateScalarResponseFunction&);
    
    AggregateScalarResponseFunction& operator=(const AggregateScalarResponseFunction&);

  protected:

    Teuchos::Array< Teuchos::RCP<ScalarResponseFunction> > responses;

  };

}

#endif // ALBANY_AGGREGATE_SCALAR_RESPONSE_FUNCTION_HPP