00001 00008 // Additional Sources or Losses ( 00010 // Well... It needs to be 4D array, so it's 1d array of 3d arrays. 00011 00012 #ifndef SourcesAndLosses_H 00013 #define SourcesAndLosses_H 00014 00015 #include "../Matrix/Matrix.h" 00016 #include "../Parameters/Parameters.h" 00017 #include "../Grid/Grid.h" 00018 00019 class SourcesAndLosses : public vector < Matrix3D<double> > { 00020 public: 00021 // Parameters 00022 ParamStructure::SL SL_parameters; 00023 00024 // Matrix1D<Matrix3D<double>> SL; 00025 00026 SourcesAndLosses(){} 00027 SourcesAndLosses(ParamStructure::SL parameters, Grid &grid, int numberOfIterations, double timeStep); 00028 void Initialize(ParamStructure::SL parameters, Grid &grid, int numberOfIterations, double timeStep); 00029 00030 }; 00031 00032 00033 #endif