#include <darcy_flow_mh_output.hh>

Collaboration diagram for DarcyFlowMHOutput:

[legend]

Classes
struct	DiffData

struct	FEData

class	OutputFields
	Standard quantities for output in DarcyFlowMH. More...

class	OutputSpecificFields
	Specific quantities for output in DarcyFlowMH - error estimates etc. More...

Public Member Functions
	DarcyFlowMHOutput (DarcyFlowInterface *flow, Input::Record in_rec)

virtual	~DarcyFlowMHOutput ()

void	output ()
	Calculate values for output. More...

Static Public Member Functions
static const Input::Type::Instance &	get_input_type (FieldSet &eq_data, const std::string &equation_name)

static const Input::Type::Instance &	get_input_type_specific ()

Protected Types
typedef const vector< unsigned int > &	ElementSetRef

Protected Member Functions
virtual void	prepare_output (Input::Record in_rec)

virtual void	prepare_specific_output (Input::Record in_rec)

void	output_internal_flow_data ()

void	compute_l2_difference ()

void	l2_diff_local (DHCellAccessor dh_cell, FEValues< 3 > &fe_values, FEValues< 3 > &fv_rt, FieldPython< 3, FieldValue< 3 >::Vector > &anal_sol, DiffData &result)
	Computes L2 error on an element. More...

Protected Attributes
DarcyFlowInterface *	darcy_flow

Mesh *	mesh_

bool	compute_errors_
	Specific experimental error computing. More...

VectorMPI	corner_pressure

VectorMPI	ele_pressure

VectorMPI	ele_piezo_head

std::vector< double >	l2_diff_pressure

std::vector< double >	l2_diff_velocity

std::vector< double >	l2_diff_divergence

std::shared_ptr< DOFHandlerMultiDim >	dh_

std::shared_ptr< DiscreteSpace >	ds

OutputFields	output_fields

OutputSpecificFields	output_specific_fields

std::shared_ptr< OutputTime >	output_stream

ofstream	raw_output_file
	Raw data output file. More...

bool	is_output_specific_fields
	Output specific field stuff. More...

struct DarcyFlowMHOutput::DiffData	diff_data

FEData	fe_data

Detailed Description

Actually this class only collect former code from postprocess.* This code depends on values stored in mesh and has to be changed to use fields or other data provided by interface to darcy_flow. We have to relay on the interface in order to allow different implementation of darcy_flow.

Principal functionalities of current postprocess are:

move computed values into mesh - this can be removed after we get new output classes, other dependent code can use directly field interface
compute interpolation to nodes - this can be external and operate directly on fields
compute water balances over materials and boundary parts - is feasible through field interface

Further functionality of this class should be:

prepare data for general output classes

Definition at line 74 of file darcy_flow_mh_output.hh.

Member Typedef Documentation

typedef const vector<unsigned int>& DarcyFlowMHOutput::ElementSetRef

protected

Definition at line 111 of file darcy_flow_mh_output.hh.

Constructor & Destructor Documentation

DarcyFlowMHOutput::DarcyFlowMHOutput	(	DarcyFlowInterface *	flow,
		Input::Record	in_rec
	)

Definition at line 127 of file darcy_flow_mh_output.cc.

DarcyFlowMHOutput::~DarcyFlowMHOutput ( )

virtual

Definition at line 233 of file darcy_flow_mh_output.cc.

Member Function Documentation

void DarcyFlowMHOutput::compute_l2_difference ( )

protected

Temporary hack. Calculate approximation of L2 norm for: 1) difference between regularized pressure and analytical solution (using FunctionPython) 2) difference between RT velocities and analytical solution 3) difference of divergence

TODO: 1) implement field objects 2) implement DG_P2 finite elements 3) implement pressure postprocessing (result is DG_P2 field) 4) implement calculation of L2 norm for two field (compute the norm and values on individual elements as P0 field)

Definition at line 490 of file darcy_flow_mh_output.cc.