Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in
Toggle navigation
Menu
Open sidebar
Mpp
MLUQ
Commits
bef6855c
Commit
bef6855c
authored
Mar 26, 2020
by
niklas.baumgarten
Browse files
added experimental changes
parent
882c69cd
Changes
9
Hide whitespace changes
Inline
Side-by-side
mlmc/src/MainProgram.cpp
View file @
bef6855c
...
...
@@ -10,10 +10,10 @@ void MainProgram::Initialize() {
problem
=
createStochasticProblem
();
assemble
=
createAssemble
();
plots
=
make_unique
<
MultilevelPlotter
>
(
*
meshes
);
mlmc
=
make_unique
<
MultilevelMonteCarlo
>
(
initLevels
,
initSampleAmount
,
meshes
,
stochField
,
assemble
,
plots
);
//
plots = make_unique<MultilevelPlotter>(*meshes);
//
//
mlmc = make_unique<MultilevelMonteCarlo>(initLevels, initSampleAmount, meshes,
//
stochField, assemble, plots);
}
void
MainProgram
::
Run
()
{
...
...
mlmc/src/MultilevelMonteCarlo.cpp
View file @
bef6855c
#include
"MultilevelMonteCarlo.h"
#include
"MultilevelMonteCarlo.h
pp
"
using
namespace
std
;
...
...
mlmc/src/MultilevelMonteCarlo.hpp
View file @
bef6855c
#ifndef MULTILEVELMONTECARLO_H
#define MULTILEVELMONTECARLO_H
#ifndef MULTILEVELMONTECARLO_H
PP
#define MULTILEVELMONTECARLO_H
PP
#include
"utils/IndentedLogger.hpp"
#include
"mc/MonteCarloElliptic.h"
#include
"mc/MonteCarloTransport.h"
#include
"mc/MonteCarloElliptic.h
pp
"
#include
"mc/MonteCarloTransport.h
pp
"
class
MultilevelMonteCarlo
{
...
...
@@ -114,4 +114,4 @@ public:
void
ShowResultsOverEpsilon
(
const
vector
<
double
>
&
epsLst
);
};
#endif //MULTILEVELMONTECARLO_H
#endif //MULTILEVELMONTECARLO_H
PP
mlmc/src/mc/MonteCarlo.cpp
View file @
bef6855c
#include
"MonteCarlo.h"
#include
"MonteCarlo.h
pp
"
void
MonteCarlo
::
Method
()
{
...
...
mlmc/src/mc/MonteCarlo.hpp
View file @
bef6855c
...
...
@@ -56,8 +56,8 @@ public:
config
.
get
(
"MCPlotting"
,
plotting
);
config
.
get
(
"functional"
,
functional
);
coarsePlot
=
plots
.
GetPlot
(
l
-
1
);
finePlot
=
plots
.
GetPlot
(
l
);
//
coarsePlot = plots.GetPlot(l - 1);
//
finePlot = plots.GetPlot(l);
config
.
get
(
"MCVerbose"
,
verbose
);
logger
=
IndentedLogger
::
GetInstance
();
...
...
mlmc/src/mc/MonteCarloElliptic.cpp
View file @
bef6855c
#include
"MonteCarloElliptic.h"
#include
"MonteCarloElliptic.h
pp
"
using
namespace
std
;
...
...
mlmc/src/mc/MonteCarloElliptic.hpp
View file @
bef6855c
#ifndef MLMC_MONTECARLOELLIPTIC_H
#define MLMC_MONTECARLOELLIPTIC_H
#include
"MonteCarlo.h"
#include
"MonteCarlo.h
pp
"
#include
"assemble/EllipticAssemble.h"
...
...
mlmc/src/mc/MonteCarloTransport.cpp
View file @
bef6855c
#include
"MonteCarloTransport.h"
#include
"MonteCarloTransport.h
pp
"
using
namespace
std
;
...
...
mlmc/src/mc/MonteCarloTransport.hpp
View file @
bef6855c
#ifndef MLMC_MONTECARLOTRANSPORT_HPP
#define MLMC_MONTECARLOTRANSPORT_HPP
#include
"MonteCarlo.h"
#include
"MonteCarlo.h
pp
"
#include
"assemble/DGTransportAssemble.h"
#include
"timestepping/TimeIntegrator.h"
//#include "timestepping/TimeIntegrator2.h"
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment