20
20
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
22
22
* @author: $LastChangedBy: dutka $
23
* @date: $LastChangedDate: 2008-06-26 13:50:17 +0200 (jeu, 26 jun 2008) $
24
* Id: $Id: t_Study_save.cxx 862 2008-06-26 11:50:17Z dutka $
23
* @date: $LastChangedDate: 2008-10-31 11:52:04 +0100 (ven 31 oct 2008) $
24
* Id: $Id: t_Study_save.cxx 995 2008-10-31 10:52:04Z dutka $
26
26
#include <iostream>
93
93
// Add another PersistentObject to the Study (here a NumericalSample)
94
NumericalSample s1(3, 2, "mySample");
94
NumericalSample s1(3, 2);
95
s1.setName("mySample");
95
96
NumericalPoint p2(2, 0.); p2[0] = 100.; p2[1] = 200.; s1[0] = p2; p2.setName("One");
96
97
NumericalPoint p3(2, 0.); p3[0] = 101.; p3[1] = 201.; s1[1] = p3; p3.setName("Two");
97
98
NumericalPoint p4(2, 0.); p4[0] = 102.; p4[1] = 202.; s1[2] = p4; p4.setName("Three");
129
130
myStudy.add(point);
131
132
// Create a NearestPointAlgorithm::Result
132
NearestPointAlgorithm::Result nearestPointAlgorithmResult(NumericalPoint(4, 1.0), 15, 1e-10, 2e-10, 3e-10, 4e-10, "nearestPointAlgorithmResult");
133
NearestPointAlgorithm::Result nearestPointAlgorithmResult(NumericalPoint(4, 1.0), 15, 1e-10, 2e-10, 3e-10, 4e-10);
134
nearestPointAlgorithmResult.setName("nearestPointAlgorithmResult");
133
135
myStudy.add(nearestPointAlgorithmResult);
135
137
// Create a Simulation::Result
136
Simulation::Result simulationResult(0.5, 0.01, 150, 4, "simulationResult");
138
Simulation::Result simulationResult(0.5, 0.01, 150, 4);
139
simulationResult.setName("simulationResult");
137
140
myStudy.add(simulationResult);
139
142
// Create a Beta distribution
213
216
formulas[0] = "a+b+c";
214
217
formulas[1] = "a-b*c";
215
218
formulas[2] = "(a+2*b^2+3*c^3)/6";
216
NumericalMathFunction analytical(input, output, formulas, "analytical");
219
NumericalMathFunction analytical(input, output, formulas);
220
analytical.setName("analytical");
217
221
myStudy.add(analytical);
219
223
// Create a QuadraticCumul algorithm
220
QuadraticCumul quadraticCumul(RandomVector(analytical, RandomVector(IndependentCopula(analytical.getInputNumericalPointDimension()), "antecedent"), "composite"), "quadraticCumul");
224
RandomVector antecedent(IndependentCopula(analytical.getInputNumericalPointDimension()));
225
antecedent.setName("antecedent");
226
RandomVector composite(analytical, antecedent);
227
composite.setName("composite");
228
QuadraticCumul quadraticCumul(composite);
229
quadraticCumul.setName("quadraticCumul");
221
230
quadraticCumul.getMeanFirstOrder();
222
231
quadraticCumul.getMeanSecondOrder();
223
232
quadraticCumul.getCovariance();
225
234
myStudy.add(quadraticCumul);
227
// Create a FORM::Result
236
// Create a FORMResult
228
237
Description input2(2);
232
241
output2[0] = "d";
233
242
Description formula2(1);
234
243
formula2[0] = "y^2-x";
235
NumericalMathFunction model(input2, output2, formula2, "sum");
236
Event event(RandomVector(model, RandomVector(Normal(2), "input"), "output"), Greater(), 1.0, "failureEvent");
244
NumericalMathFunction model(input2, output2, formula2);
245
model.setName("sum");
246
RandomVector input3(Normal(2));
247
input3.setName("input");
248
RandomVector output3(model, input3);
249
output3.setName("output");
250
Event event(output3, Greater(), 1.0);
251
event.setName("failureEvent");
237
252
NumericalPoint designPoint(2, 0.0);
238
253
designPoint[0] = 1.0;
239
FORM::Result formResult(NumericalPoint(2, 1.0), event, false, "formResult");
254
FORMResult formResult(NumericalPoint(2, 1.0), event, false);
255
formResult.setName("formResult");
240
256
formResult.getImportanceFactors();
241
257
formResult.getEventProbabilitySensitivity();
242
258
myStudy.add(formResult);
244
// Create a SORM::Result
245
SORM::Result sormResult(NumericalPoint(2, 1.0), event, false, "sormResult");
260
// Create a SORMResult
261
SORMResult sormResult(NumericalPoint(2, 1.0), event, false);
262
sormResult.setName("sormResult");
246
263
sormResult.getEventProbabilityBreitung();
247
264
sormResult.getEventProbabilityHohenBichler();
248
265
sormResult.getEventProbabilityTvedt();