# Link Uqlab and Castem

**URL:** <https://uqworld.org/t/link-uqlab-and-castem/959>\
**Category:** Community Q&A and How To\
**Created:** [December 1, 2020, 8:46am UTC](https://uqworld.org/t/link-uqlab-and-castem/959 "2020-12-01T08:46:32Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![IHM](https://uqworld.org/letter_avatar_proxy/v4/letter/i/fbc32d/32.png) [@IHM](https://uqworld.org/u/IHM)\
**Post date:** [December 1, 2020, 8:46am UTC](https://uqworld.org/t/link-uqlab-and-castem/959/1 "2020-12-01T08:46:32Z")

</div>

Hello, I am a beginner with the Uqlab software and I want to perform the sensitivity analysis of a finite element digital model implemented in the Castem software. My problem is that I can’t make the connection between the two software I read the Uqlink user manual I tried to follow the examples but I have no value. However, I was able to link my Castem to MatLab with the following code and the calculation goes correctly.  
n\_sim = 10;  
 F\_ct Ft.m = 3.5e6; Ft.cov = 0.13; Ft.type = 'Lognormal'; Mod\_E  
mod\_E.m = 35e9;  
mod\_E.cov = 0.12;  
mod\_E.type = ‘Lognormal’;  
% Initiate the vector to store results

Ft.data = generRV(n\_sim,Ft);  
mod\_E.data = generRV(n\_sim,mod\_E);

sigma = zeros(1,n\_sim);  
fleche = zeros(1,n\_sim);  
timeCPU = 0;

for i\_sim = 1 : n\_sim

```
%% Create vector x for ith simulation
x = [Ft.data(i_sim) mod_E.data(i_sim)];

%% Write vector x to file
fid = fopen('inputstat.txt','w');
%
n = max(size(x));
for i = 1:n    
    fprintf(fid,'%20.16E\n',x(i));
end
fclose(fid);

%% Call Castem to calculate
tic    

dos('castem18 essais_stat_OK.dgibi');
  
timeCPU = timeCPU + toc;

%% Get Output
File = 'resultstat.out';            
A = textread(File,'%f');      

 fleche(i_sim) = A(1,1);
 sigma(i_sim) = A(2,1);      

```

end

Here is my uqlab code.

clearvars  
rng(1,‘twister’)  
uqlab

ModelOpts.Type = ‘UQLink’;  
ModelOpts.Name = ‘CASTEM Wrapper’;

EXECUTABLE = ‘C:\Cast3M\PCW\_18\bin’;  
INPUTFILE = ‘essais\_stat\_OK.dgibi’;

COMMANDLINE = sprintf(’%s %s’,EXECUTABLE,INPUTFILE);

ModelOpts.Command = COMMANDLINE;

ModelOpts.Template = ‘essais\_stat\_OK.dgibi.tpl’;

ModelOpts.Output.FileName = ‘resultstat.out’;  
ModelOpts.Output.Parser = ‘readresultOutput’;

ModelOpts.ExecutionPath = fullfile(uq\_rootPath,…  
‘Examples’,‘UQLink’,‘sens\_stat’) ;

ModelOpts.Format = {’%1.8f’};

ModelOpts.Display = ‘quiet’;

monModelCASTEM = uq\_createModel(ModelOpts);

InputOpts.Marginals(1).Name = ‘E’;  
InputOpts.Marginals(1).Type = ‘Lognormal’;  
InputOpts.Marginals(1).Moments = [35e9 4.2e9]; % (Pa)

InputOpts.Marginals(2).Name = ‘Ft’;  
InputOpts.Marginals(2).Type = ‘Lognormal’;  
InputOpts.Marginals(2).Moments = [3.5e6 0.455e6];

myInput = uq\_createInput(InputOpts);

X = uq\_getSample(10,‘LHS’);

Y = uq\_evalModel(monModelCASTEM,X);

Thank you in advance for your return

---

<div class="post-metadata">

**Author:** ![damarginal](https://uqworld.org/user_avatar/uqworld.org/damarginal/32/73_2.png) [@damarginal](https://uqworld.org/u/damarginal)\
**Post date:** [January 15, 2021, 3:10pm UTC](https://uqworld.org/t/link-uqlab-and-castem/959/2 "2021-01-15T15:10:34Z")

</div>

This is a duplicate of the following (resolved) topic:

[https://uqworld.org/t/linking-uqlab-and-castem](https://uqworld.org/t/linking-uqlab-and-castem)

This topic is now closed.

---

<div class="post-metadata">

**Author:** ![damarginal](https://uqworld.org/user_avatar/uqworld.org/damarginal/32/73_2.png) [@damarginal](https://uqworld.org/u/damarginal)\
**Post date:** [January 15, 2021, 3:10pm UTC](https://uqworld.org/t/link-uqlab-and-castem/959/3 "2021-01-15T15:10:53Z")

</div>