Skip to content
Snippets Groups Projects
Commit b6f7dff4 authored by Tom Fischer's avatar Tom Fischer
Browse files

Fixed include dependencies.

parent 8afc6023
No related branches found
No related tags found
No related merge requests found
......@@ -13,9 +13,6 @@
*/
#include "FEMCondition.h"
#include "ProcessInfo.h"
#include "GEOObjects.h" //for SourceTerm
FEMCondition::FEMCondition(const std::string &geometry_name, CondType t)
: _type(t), _geoName("[unspecified]"), _associated_geometry(geometry_name)
......@@ -26,16 +23,20 @@ FEMCondition::FEMCondition(const std::string &geometry_name, CondType t)
this->setProcessDistributionType(FiniteElement::INVALID_DIS_TYPE);
}
FEMCondition::FEMCondition(const std::string &geometry_name, FiniteElement::ProcessType pt,
FiniteElement::PrimaryVariable pv, GeoLib::GEOTYPE gt, const std::string &gn,
FiniteElement::DistributionType dt, CondType ct)
FEMCondition::FEMCondition(const std::string &geometry_name,
FiniteElement::ProcessType pt,
FiniteElement::PrimaryVariable pv,
GeoLib::GEOTYPE gt,
const std::string &gn,
FiniteElement::DistributionType dt,
CondType ct)
: ProcessInfo(pt, pv/*, NULL*/), GeoInfo(gt, NULL), DistributionInfo(dt), _type(ct),
_geoName(gn), _associated_geometry(geometry_name)
{
}
FEMCondition::FEMCondition(const FEMCondition &cond, CondType t)
: ProcessInfo(cond.getProcessType(), cond.getProcessPrimaryVariable()/*, NULL*/),
: ProcessInfo(cond.getProcessType(), cond.getProcessPrimaryVariable() /*, NULL*/),
GeoInfo(cond.getGeomType(), cond.getGeoObj()),
DistributionInfo(cond.getProcessDistributionType()),
_type(t),
......@@ -67,6 +68,6 @@ void FEMCondition::setDisValues(const std::vector< std::pair<size_t, double> > &
nodes.push_back(dis_values[i].first);
values.push_back(dis_values[i].second);
}
this->_disNodes=nodes;
this->_disValues=values;
this->_disNodes = nodes;
this->_disValues = values;
}
......@@ -15,22 +15,16 @@
#ifndef FEMCONDITION_H
#define FEMCONDITION_H
#include <vector>
#include "DistributionInfo.h"
#include "GeoInfo.h"
#include "GeoObject.h"
#include "Point.h"
#include "ProcessInfo.h"
#include <vector>
class CBoundaryCondition;
class CInitialCondition;
class CSourceTerm;
//class GEOObjects;
class GridAdapter;
#include "GEOObjects.h"
/**
* \brief Adapter class for handling FEM Conditions in the user Interface
*/
......@@ -47,11 +41,12 @@ public:
FEMCondition(const std::string &geometry_name, CondType t = UNSPECIFIED);
FEMCondition(const std::string &geometry_name,
FiniteElement::ProcessType pt = FiniteElement::INVALID_PROCESS,
FiniteElement::PrimaryVariable pv = FiniteElement::INVALID_PV,
GeoLib::GEOTYPE gt = GeoLib::INVALID,
const std::string &gn = "[unspecified]",
FiniteElement::DistributionType dt = FiniteElement::INVALID_DIS_TYPE, CondType ct = UNSPECIFIED);
FiniteElement::ProcessType pt = FiniteElement::INVALID_PROCESS,
FiniteElement::PrimaryVariable pv = FiniteElement::INVALID_PV,
GeoLib::GEOTYPE gt = GeoLib::INVALID,
const std::string &gn = "[unspecified]",
FiniteElement::DistributionType dt = FiniteElement::INVALID_DIS_TYPE,
CondType ct = UNSPECIFIED);
FEMCondition(const FEMCondition &cond, CondType t);
~FEMCondition() {}
......@@ -75,7 +70,8 @@ public:
void setConstantDisValue(double disValue) {_disValues.clear(); _disValues.push_back(disValue); }
/// Sets a vector of values specifying the distribution.
void setDisValues(const std::vector<std::size_t> &disNodes, const std::vector<double> &disValues)
void setDisValues(const std::vector<std::size_t> &disNodes,
const std::vector<double> &disValues)
{
_disNodes = disNodes;
_disValues = disValues;
......@@ -86,7 +82,7 @@ public:
void setDisValues(const std::vector< std::pair<std::size_t, double> > &dis_values);
/// Removes all distribution values
void clearDisValues() { _disValues.resize(0); };
void clearDisValues() { _disValues.resize(0); }
/// Sets the name of the geo-object the condition is assigned to.
void setGeoName(std::string geoName) { _geoName = geoName; }
......
......@@ -16,6 +16,8 @@
// GeoLib
#include "AABB.h"
#include "Color.h"
#include "Polyline.h"
#include "Surface.h"
#include "FEMCondition.h"
#include "VtkConditionSource.h"
......@@ -29,9 +31,9 @@
#include <vtkPoints.h>
#include <vtkPolyData.h>
#include <vtkPolygon.h>
#include <vtkProperty.h>
#include <vtkSmartPointer.h>
#include <vtkStreamingDemandDrivenPipeline.h>
#include <vtkProperty.h>
#include <vtkLookupTable.h>
......@@ -131,7 +133,8 @@ int VtkConditionSource::RequestData( vtkInformation* request,
dis_type_value = static_cast<vtkIdType>(_dis_type_map.size());
_dis_type_map.insert(std::pair<FiniteElement::DistributionType, size_t>(type, dis_type_value));
}
else dis_type_value = it->second;
else
dis_type_value = it->second;
if ((*_cond_vec)[n]->getGeomType() == GeoLib::POINT)
{
......@@ -166,7 +169,8 @@ int VtkConditionSource::RequestData( vtkInformation* request,
}
else if ((*_cond_vec)[n]->getGeomType() == GeoLib::POLYLINE)
{
const GeoLib::Polyline* ply = static_cast<const GeoLib::Polyline*>((*_cond_vec)[n]->getGeoObj());
const GeoLib::Polyline* ply =
static_cast<const GeoLib::Polyline*>((*_cond_vec)[n]->getGeoObj());
const int nPoints = ply->getNumberOfPoints();
newLines->InsertNextCell(nPoints);
double value (0);
......@@ -182,7 +186,7 @@ int VtkConditionSource::RequestData( vtkInformation* request,
scalars->InsertNextValue(dis_values[0]);
else if (type == FiniteElement::LINEAR || type == FiniteElement::LINEAR_NEUMANN)
{
for (size_t j = 0; j < dis_values.size(); j ++)
for (size_t j = 0; j < dis_values.size(); j++)
{
if (static_cast<int>(dis_nodes[j]) == i)
value = dis_values[j];
......@@ -216,7 +220,8 @@ int VtkConditionSource::RequestData( vtkInformation* request,
if (point_idx_map[point_index] == -1)
{
point_idx_map[point_index] = pnt_id;
newPoints->InsertNextPoint((*_points)[point_index]->getCoords());
newPoints->InsertNextPoint(
(*_points)[point_index]->getCoords());
aPolygon->GetPointIds()->SetId(j, pnt_id);
distypes->InsertNextValue(dis_type_value);
......@@ -246,7 +251,7 @@ int VtkConditionSource::RequestData( vtkInformation* request,
else if ((*_cond_vec)[n]->getGeomType() == GeoLib::INVALID)
{
size_t nValues = dis_values.size();
for (size_t i=0; i<nValues; i++)
for (size_t i = 0; i < nValues; i++)
{
//vtkIdType pid = newPoints->InsertNextPoint((*_points)[dis_nodes[i]]->getData());
vtkIdType pid = newPoints->InsertNextPoint((*_points)[i]->getCoords());
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment