Equation with "if statements"

[engrmessi]

Hi all, I am writing an equation (expression) that comprises "if statements" in my model as shown in the attached file. I just want to know if the equation is written correctly. if not I need correction on how best to write it in MOOSE. Thanks!

[GiudGiud]

hello

no it's not. I dont know if && works so I d do this:

expression = 'if (t > 0, if(t<to, L0 * (1 + accslip / 0.0111), if (t<to + tf, L0 * (1 + accslip / 0.0111) * exp(-k * (t - to)), 0)), 0)'

[GiudGiud]

regardless of if && works, you need to keep this if (a, b, c) structure, even when nesting ifs: if (a1, if(a2, b, d), c)

[engrmessi]

Thank you for your response
But the first part of second inequality condition "to < t" was not accounted for. Is it Ok to ignore it? Also some of the inequality conditions has something like "less than or equal to"

[GiudGiud]

can you paste the rest of the input? this does not look like an error on this line, but on two lines in the input

[engrmessi]

Please is there a way I can send the input script to you privately since the project is confidential to my research lab.?

[GiudGiud]

sure. guillaume.giudicelli@inl.gov

[GiudGiud]

I looked at the inputs. You have in the _repassivat.i input several [/] instead of [] for the closing block. [/] is causing an error.

I did not find any issue in diffcorr_c.i though I cannot run it because I do not have the two txt files necessary (one slip systems and one PropertyReadFile)

[engrmessi]

My code is working now
Thank you very much.

[ArashPartow]

@engrmessi / @GiudGiud : If the following PR #13194 ever gets merged, you can in the future have the following expression:

" L :=                                                                       "
"   switch                                                                   "
"   {                                                                        "
"      case t < 0         : 0;                                               "
"      case t < to        : L0 * (1 + accslip / 0.0111);                     "
"      case t < (to + tf) : L0 * (1 + accslip / 0.0111) * exp(-k * (t - to); "
"      default            : 0;                                               "
"   };                                                                       "
"                                                                            ";