A question in cfdAssembleStressTerm

[liuqi45]

Look at the code for calculating geoDiff_f
In the 78 row,''n = cfdUnit(Sf);''
In the 79 row,''e = cfdUnit(CF);''
In the 81 row, ''CF_limited = max(dot(n',CF'), 0.05*cfdMag(CF)')';''
In the 84 row,''magCF = cfdMag(CF_limited);''
In the 105 row,''geoDiff_f = cfdMag(Ef)./magCF;''

I think the CF_limited is wrong, it should be ''CF_limited = max(dot(e',CF'), 0.05*cfdMag(CF)')';''

[ruohai0925]

Hi @liuqi45 , it seems that the authors are not actively maintaining this repo now. I am also investigating this code these days. If you can send me an email to zdsjtu@gmail.com, maybe we can try to learn the codes together and exchange some experiences.

[shanjierenyidp]

HI liuqi,
I think you are correct. The original code gives CF_limited as a scalar, which is wrong. You correction makes it back to an array.