Homec4science
Diffusion Guggisberg-Liardon - Numerical integration (master) Browse

R7917/DomainCartesian.cpp
0f8ddfe1ee20master

Last Change
Raw File

DomainCartesian.cpp
View Options

#include "DomainCartesian.hpp"
// Constructor and destructor
DomainCartesian::DomainCartesian() {}
// Standard constructor
DomainCartesian::DomainCartesian(vector<vector<double>> boundaries_, vector<int> subdivision_) {
boundaries = boundaries_;
dimension = (int)boundaries.size();
subdivision = subdivision_;
// Exception : the amount of subdivisions must be over 0.
for(int i = 0; i < dimension; i++){
if (subdivision[i]<1){
OutOfRangeException error("The amount of subdivisions must be over 0.");
error.PrintDebug();
exit(0);
}
}
// Exception : only 1D and 2D integration methods have been implemented
if (dimension > 2){
OutOfRangeException error("The integration methods for dimensions over 2 have not been implemented.");
error.PrintDebug();
exit(0);
}
}
// Simplified constructor for 1D domain
DomainCartesian::DomainCartesian(double a_, double b_, int subdivision_) {
// Transform the input in the standard form of the cartesian domain
vector<double> dim1;
dim1.push_back(a_);
dim1.push_back(b_);
boundaries.push_back(dim1);
dimension = 1;
subdivision.push_back(subdivision_);
// Exception : the amount of subdivisions must be over 0.
if (subdivision[0]<1){
OutOfRangeException error("The amount of subdivisions must be over 0.");
error.PrintDebug();
exit(0);
}
}
// Get number of points in the discretization of the domain
int DomainCartesian::GetNumberPoints() {
// Iteration for on each dimension
int total = 1;
for (int i = 0; i < (int)subdivision.size(); i++) {
total *= subdivision[i]; // Multiplication of the points from previous dimensions with those in the i dimension
}
return total;
}
c4science · Help