# Vectorized mfile-based model with an if statement

**URL:** <https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199>\
**Category:** Community Q&A and How To\
**Tags:** model, metamodeling, reliability\
**Created:** [May 25, 2021, 3:26pm UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199 "2021-05-25T15:26:28Z")\
**Posts on this page:** 9\
**Page:** 1

<div class="post-metadata">

**Author:** ![Neryvaldo\_Galvao](https://uqworld.org/user_avatar/uqworld.org/neryvaldo_galvao/32/597_2.png) [@Neryvaldo\_Galvao](https://uqworld.org/u/Neryvaldo_Galvao)\
**Post date:** [May 25, 2021, 3:26pm UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/1 "2021-05-25T15:26:28Z")

</div>

Dear UQLab Community

I have a mfile-based model with this if statement that needs to be true so I can evaluate a performance function according to different metamodel previously established.

It is my understanding that the mfile is vectorized, improving the performance of the code by several orders of magnitude, since it allows the evaluation of N different realizations in a single call.

Nevertheless, I notice that the if statement that I am using the mfile is never true, but it should be for some runs. When the mfile was not vectorized I could confirm this.

Is this something that can be done? If so, please let me know what I am missing.  
Please find enclosed a short example of the code in the file I am using.  
[mfilebased\_model.m](https://uqworld.org/uploads/short-url/qTJJM0GWeNNXtgP17Mk7o7SA90g.m) (308 Bytes)

Furthermore, I should say that I am using subset simulation as the reliability analysis method.

I thank you in advance for your time and support  
Best Regards,  
Neryvaldo Galvão

---

<div class="post-metadata">

**Author:** ![styfen.schaer](https://uqworld.org/user_avatar/uqworld.org/styfen.schaer/32/849_2.png) [@styfen.schaer](https://uqworld.org/u/styfen.schaer)\
**Post date:** [May 25, 2021, 5:26pm UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/2 "2021-05-25T17:26:22Z")

</div>

Dear Neryvaldo Galvão,

Since your code is short, I take the liberty of posting it below to make it more convenient for other users.

The fact that your IF statement never evaluates to true after you vectorized your code makes me wonder if your condition is correct. Your current condition is equivalent to:  
`all(X(:,18) <= X(:,19))`  
However, you probably wanted the following condition:  
`any(X(:,18) <= X(:,19))`  
Can you confirm that your condition is correct?

Best regards,  
Styfen

```matlab
function Y = myHRAEventTree(X)

load(Metamodel_1)
load(Metamodel_2)
load(Metamodel_3)

if X(:,18)<=X(:,19)
    Y = uq_evalModel(Metamodel_1,X(:,1:16));
elseif X(:,18)<=X(:,19)+X(:,20)
    Y = uq_evalModel(Metamodel_2,X(:,1:16));
else    
    Y = uq_evalModel(Metamodel_3,X(:,1:16));
end

end

```

---

<div class="post-metadata">

**Author:** ![Neryvaldo\_Galvao](https://uqworld.org/user_avatar/uqworld.org/neryvaldo_galvao/32/597_2.png) [@Neryvaldo\_Galvao](https://uqworld.org/u/Neryvaldo_Galvao)\
**Post date:** [May 26, 2021, 9:06am UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/3 "2021-05-26T09:06:21Z")

</div>

Dear Styfen Schaer

Thank you for your support

all(X(:,18) \<= X(:,19)) is not what I meant to do ☹. Nevertheless, I not sure if any(X(:,18) \<= X(:,19)) will work as well 🤔.

What I need is: for any realization N of a single call, if X(18)\<=X(19) the Metamodel\_1 should be used to get the output of interest i.e., Y for that single realization where the condition “if” is verified.

The non-vectorized code that I had before was something like the following and it was working fine, but the performance was terrible. 🥵

```
function Y = myHRAEventTree(X)

load(Metamodel_1)
load(Metamodel_2)
load(Metamodel_3)

for i=1:size(X,1)
if X(i,18)<=X(i,19)
    Y(i) = uq_evalModel(Metamodel_1,X(i,1:16));
elseif X(i,18)<=X(i,19)+X(i,20)
    Y(i) = uq_evalModel(Metamodel_2,X(i,1:16));
else   
    Y(i) = uq_evalModel(Metamodel_3,X(i,1:16));
end
end

     
end

```

How can I convert this code to a vectorized mfile-based model? 🙏

Looking forward to hearing from you.

Thank you in advance for your support.  
Best regards,  
Neryvaldo Galvão

---

<div class="post-metadata">

**Author:** ![styfen.schaer](https://uqworld.org/user_avatar/uqworld.org/styfen.schaer/32/849_2.png) [@styfen.schaer](https://uqworld.org/u/styfen.schaer)\
**Post date:** [May 26, 2021, 9:45am UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/4 "2021-05-26T09:45:07Z")

</div>

Dear Neryvaldo Galvão;

any() will not work in this case. However, I think you’re looking for something like below:

Best regards,  
Styfen

```matlab
function Y = myHRAEventTree(X)

load(Metamodel_1)
load(Metamodel_2)
load(Metamodel_3)

cond1 = X(:,18) <= X(:,19);
X1 = X(cond1,16);

cond2 = X(:,18) <= X(:,19) + X(:,20);
X2 = X(cond2,16);

cond3 = ~cond1 & ~cond2;
X3 = X(cond3,16);

Y1 = uq_evalModel(Metamodel_1, X1);
Y2 = uq_evalModel(Metamodel_2, X2);
Y3 = uq_evalModel(Metamodel_3, X3);

Y = zeros(size(X,1),1);
Y(cond1) = Y1;
Y(cond2) = Y2;
Y(cond3) = Y3;
     
end

```

---

<div class="post-metadata">

**Author:** ![Neryvaldo\_Galvao](https://uqworld.org/user_avatar/uqworld.org/neryvaldo_galvao/32/597_2.png) [@Neryvaldo\_Galvao](https://uqworld.org/u/Neryvaldo_Galvao)\
**Post date:** [May 26, 2021, 11:15am UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/5 "2021-05-26T11:15:59Z")

</div>

Dear Styfen Schaer

Thank you very much 😀. It looks like the code you suggested is working fine 👍.

Nevertheless, I need to understand a little more about it. Where can I read more about the approach you suggesting?

Furthermore, lets say I need to do the following multiplication X(i,6) = X(i,17) \* X(i,15) when cond1 = X(:,18) \<= X(:,19) is true. How can I do this in this vectorized approach you proposing? I have tried a couple of things, but sadly it is not working.

Many thanks for your support  
Best Regards,  
Neryvaldo Galvão

---

<div class="post-metadata">

**Author:** ![styfen.schaer](https://uqworld.org/user_avatar/uqworld.org/styfen.schaer/32/849_2.png) [@styfen.schaer](https://uqworld.org/u/styfen.schaer)\
**Post date:** [May 26, 2021, 11:42am UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/6 "2021-05-26T11:42:07Z")

</div>

Dear Neryvaldo Galvão,

I’m glad it worked.  
The proposed approach is not directly related to UQLab but to Matlab.  
However, the following link shows some examples of how to use logical operators to select submatrizes:

[https://ch.mathworks.com/help/matlab/matlab\_prog/find-array-elements-that-meet-a-condition.html](https://ch.mathworks.com/help/matlab/matlab_prog/find-array-elements-that-meet-a-condition.html)

Regarding your problem you may have missed to use the elementwise operator.  
The following should work:

```matlab
cond1 = X(:,18) <= X(:,19);
X(cond1,6) = X(cond1,17) .* X(cond1,15);

```

Best regards,  
Styfen

---

<div class="post-metadata">

**Author:** ![Neryvaldo\_Galvao](https://uqworld.org/user_avatar/uqworld.org/neryvaldo_galvao/32/597_2.png) [@Neryvaldo\_Galvao](https://uqworld.org/u/Neryvaldo_Galvao)\
**Post date:** [May 26, 2021, 4:47pm UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/7 "2021-05-26T16:47:02Z")

</div>

Dear Styfen Schaer

Once again thank you very much. Almost everything is working fine 😀.

But it seems to me that these lines:

```
Y(cond1) = Y1;
Y(cond2) = Y2;
Y(cond3) = Y3;

```

are overwriting each other. At the end Y(cond3)=Y3 overwrites the values provided by Y(cond2) = Y2 and this last one overwrite the Y(cond1) = Y1. Making my efforts to separate the estimation of Y with different metamodel inconsequential 😰.

```
Y1 = uq_evalModel(Metamodel_1, X1);
Y2 = uq_evalModel(Metamodel_2, X2);
Y3 = uq_evalModel(Metamodel_3, X3);

```

Please let me know what are your thoughts on this.

Best Regards,  
Neryvaldo Galvão

---

<div class="post-metadata">

**Author:** ![styfen.schaer](https://uqworld.org/user_avatar/uqworld.org/styfen.schaer/32/849_2.png) [@styfen.schaer](https://uqworld.org/u/styfen.schaer)\
**Post date:** [May 26, 2021, 6:02pm UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/8 "2021-05-26T18:02:30Z")

</div>

> [@Neryvaldo\_Galvao](#):
>
> Please let me know what are your thoughts on this.

Well, I think that’s a tiny mistake on my part. 😉  
The second condition should read:

```matlab
cond2 = X(:,18) <= X(:,19) + X(:,20) & ~cond1;

```

You should double check that all conditions are correct now.

Best regards,  
Styfen

---

<div class="post-metadata">

**Author:** ![Neryvaldo\_Galvao](https://uqworld.org/user_avatar/uqworld.org/neryvaldo_galvao/32/597_2.png) [@Neryvaldo\_Galvao](https://uqworld.org/u/Neryvaldo_Galvao)\
**Post date:** [May 27, 2021, 9:48am UTC](https://uqworld.org/t/vectorized-mfile-based-model-with-an-if-statement/1199/9 "2021-05-27T09:48:47Z")

</div>

Great!! Looks like that everything is working fine. I still have to run a few more tests though.

Thank you very much for the super support 💪 💪. It is great to know I have you on the other side.

Regards,  
Neryvaldo Galvão