GRINS::SteadyVisualization Class Reference

#include <steady_visualization.h>

Inheritance diagram for GRINS::SteadyVisualization:

Collaboration diagram for GRINS::SteadyVisualization:

Public Member Functions
	SteadyVisualization (const GetPot &input, const libMesh::Parallel::Communicator &comm LIBMESH_CAN_DEFAULT_TO_COMMWORLD)
	~SteadyVisualization ()
virtual void	output_residual (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system, const unsigned int time_step, const libMesh::Real time)
virtual void	output_residual_sensitivities (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system, const libMesh::ParameterVector &params, const unsigned int time_step, const libMesh::Real time)
virtual void	output_adjoint (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system, const unsigned int time_step, const libMesh::Real time)
virtual void	output_solution_sensitivities (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system, const libMesh::ParameterVector &params, const unsigned int time_step, const libMesh::Real time)
Public Member Functions inherited from GRINS::Visualization
	Visualization (const GetPot &input, const libMesh::Parallel::Communicator &comm LIBMESH_CAN_DEFAULT_TO_COMMWORLD)
virtual	~Visualization ()
void	output (SharedPtr< libMesh::EquationSystems > equation_system)
void	output (SharedPtr< libMesh::EquationSystems > equation_system, const unsigned int time_step, const libMesh::Real time)
void	output_residual (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system)
void	output_residual_sensitivities (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system, const libMesh::ParameterVector &params)
void	output_adjoint (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system)
void	output_solution_sensitivities (SharedPtr< libMesh::EquationSystems > equation_system, GRINS::MultiphysicsSystem *system, const libMesh::ParameterVector &params)
void	dump_visualization (SharedPtr< libMesh::EquationSystems > equation_system, const std::string &filename_prefix, const libMesh::Real time)

Additional Inherited Members
Protected Attributes inherited from GRINS::Visualization
std::string	_vis_output_file_prefix
std::vector< std::string >	_output_format

Detailed Description

Definition at line 33 of file steady_visualization.h.

Constructor & Destructor Documentation

GRINS::SteadyVisualization::SteadyVisualization	(	const GetPot &	input,
		const libMesh::Parallel::Communicator &comm	LIBMESH_CAN_DEFAULT_TO_COMMWORLD
	)

Definition at line 41 of file steady_visualization.C.

    : Visualization(input, comm)
  {
    return;
  }

GRINS::SteadyVisualization::~SteadyVisualization

(

)

Definition at line 48 of file steady_visualization.C.

  {
    return;
  }

Member Function Documentation

void GRINS::SteadyVisualization::output_adjoint ( SharedPtr< libMesh::EquationSystems >  equation_system, GRINS::MultiphysicsSystem *  system, const unsigned int  time_step, const libMesh::Real  time  )

virtual

Implements GRINS::Visualization.

Definition at line 104 of file steady_visualization.C.

References GRINS::Visualization::_vis_output_file_prefix, GRINS::Visualization::dump_visualization(), GRINS::CompositeQoI::get_qoi(), GRINS::CompositeQoI::n_qois(), and GRINS::QoIBase::name().

  {
    const libMesh::DifferentiableQoI* raw_qoi = system->get_qoi();
    const CompositeQoI* qoi = dynamic_cast<const CompositeQoI*>( raw_qoi );

    unsigned int n_qois = qoi->n_qois();

    for( unsigned int q = 0; q < n_qois; q++ )
      {
        libMesh::NumericVector<libMesh::Number>& dual_solution = system->get_adjoint_solution(q);

        const std::string& qoi_name = qoi->get_qoi(q).name();
        std::string filename = this->_vis_output_file_prefix+"_adjoint_"+qoi_name;

        system->solution->swap( dual_solution );
        equation_system->update();

        this->dump_visualization( equation_system, filename, 0.0 );

        // Now swap back and reupdate
        system->solution->swap( dual_solution );
        equation_system->update();
      }
  }

void GRINS::SteadyVisualization::output_residual ( SharedPtr< libMesh::EquationSystems >  equation_system, GRINS::MultiphysicsSystem *  system, const unsigned int  time_step, const libMesh::Real  time  )

virtual

Implements GRINS::Visualization.

Definition at line 53 of file steady_visualization.C.

References GRINS::Visualization::_vis_output_file_prefix, and GRINS::Visualization::dump_visualization().

  {
    std::string filename = this->_vis_output_file_prefix+"_residual";

    // Idea is that this->rhs stashes the residual. Thus, when we swap
    // with the solution, we should be dumping the residual. Then, we swap
    // back once we're done outputting.

    // Swap solution with computed residual
    system->solution->swap( *(system->rhs) );
    equation_system->update();

    this->dump_visualization( equation_system, filename, 0.0 );

    // Now swap back and reupdate
    system->solution->swap( *(system->rhs) );
    equation_system->update();

    return;
  }

void GRINS::SteadyVisualization::output_residual_sensitivities ( SharedPtr< libMesh::EquationSystems >  equation_system, GRINS::MultiphysicsSystem *  system, const libMesh::ParameterVector &  params, const unsigned int  time_step, const libMesh::Real  time  )

virtual

Implements GRINS::Visualization.

Definition at line 78 of file steady_visualization.C.

  {
    for (unsigned int p=0; p != params.size(); ++p)
      {
        std::stringstream pstr;
        pstr << p;

        std::string filename =
          this->_vis_output_file_prefix+"_dRdp"+pstr.str();

        // Swap solution with precomputed sensitivity rhs
        system->solution->swap(system->get_sensitivity_rhs(p));
        equation_system->update();

        this->dump_visualization( equation_system, filename, 0.0 );

        // Now swap back and reupdate
        system->solution->swap(system->get_sensitivity_rhs(p));
        equation_system->update();
      }
  }

void GRINS::SteadyVisualization::output_solution_sensitivities ( SharedPtr< libMesh::EquationSystems >  equation_system, GRINS::MultiphysicsSystem *  system, const libMesh::ParameterVector &  params, const unsigned int  time_step, const libMesh::Real  time  )

virtual

Implements GRINS::Visualization.

Definition at line 133 of file steady_visualization.C.

  {
    for (unsigned int p=0; p != params.size(); ++p)
      {
        std::stringstream pstr;
        pstr << p;

        std::string filename =
          this->_vis_output_file_prefix+"_dudp"+pstr.str();

        // Swap solution with precomputed sensitivity solution
        system->solution->swap(system->get_sensitivity_solution(p));
        equation_system->update();

        this->dump_visualization( equation_system, filename, 0.0 );

        // Now swap back and reupdate
        system->solution->swap(system->get_sensitivity_solution(p));
        equation_system->update();
      }
  }

---

The documentation for this class was generated from the following files:

• src/visualization/include/grins/steady_visualization.h
• src/visualization/src/steady_visualization.C