opustc_couplex — Returns the output for the OPUS - COUPLEX test case.
Y = opustc_couplex ( X ) Y = opustc_couplex ( X , D )
a np-by-15 matrix, where np is the number of experiments to perform
a parameter list with the determinist parameters
the subdirectory of the temporary directory where the simulation must be performed (default = "couplex")
set to true to make the simulation verbose (default %f)
the name of the CASTEM executable (default "castem09")
a np-by-168 matrix
Launch Cast3m and perform the COUPLEX simulation.
Prerequirements: The Cast3m software must be installed on the machine.
http://www-cast3m.cea.fr/cast3m/index.jsp
In order to check that cast3m is correctly installed on your system, open a terminal and type, on Windows:
castem09.bat
on Linux:
castem09
The data associated with this benchmark is provided in the "data/couplex" subdirectory of this module. This directory contains the following files:
couplex.dgibi.ref (input file of CASTEM): contains the template couplex input file for CASTEM. The content of this file is modified, so that the actual data correspond to the input X.
filePlanExperience.txt (input file of CASTEM): the design of experiments data file. Each row of this file correspond to an experiment to perform. Each column of this file correspond to an input parameter.
iode.dat (input file of CASTEM): the temporal evolution of the iode, the source term of the partial differential equation.
fileReponse.txt (output file of CASTEM): the output of the simulation. Each row of this file correspond to an experiment which has been performed. Each column of this file correspond to an output parameter.
resu_conc.dat (output file of CASTEM): the concentration of different species depending on the time and the point of measure (P1 to P8).
resu_flux.dat (output file of CASTEM): the flux of iode through different surfaces of the mesh.
See the data/couplex/opuscouplex.pdf document for a more complete manual of this benchmark.
CASTEM prints messages in the console. In some cases, CASTEM cannot read the data files. In this case, CASTEM reads from the console interactively, which obviously fails. To detect these situations, open the console where Scilab runs in with the consolebox function or the "F12" Scilab shortcut.
// Check that the "CASTEM" environment variable is well defined getenv("CASTEM") // CASTEM prins messages in the console. consolebox("on") // or F12. // Perform one experiment X = [25.0 3.0e-6 6.0 3.0e-5 1e-6 5.0e-4 50. 50. 50. 50. 21. 51. 3. 110. 160. ]; Y = opustc_couplex ( X ) // Enable verbose logging: D = init_param(); D = add_param(D,"verbose",%t); Y = opustc_couplex ( X , D ); // Customize directory D = init_param(); D = add_param(D,"directory","couplex_1"); D = add_param(D,"verbose",%t); Y = opustc_couplex ( X , D ); // See the data associated with the couplex benchmark. // This data is in the "data/couplex" subdirectory of this module. path = opustc_getpath ( ); ls(fullfile(path,"data","couplex"))
"OPUS, Cas test CEA, Propagation des incertitudes dans le probleme COUPLEX-1", JM. Martinez, Ph. Montarnal, M. Lamoureux, O. Le Maitre, 6 avril 2010
"Analyse d'incertitudes : application de la theorie des polynomes du chaos a l'analyse de surete du stockage de dechets nucleaires", Ecole Centrale Paris, Paul Balanca, Anne-Celine Boulanger, Marc Dupin, Cyril Labbe, Etude encadree par Florian De Vuyst