forAll 起到什么作用呢？

[mohui]

forAll(dgdt, celli)
        {
            if (dgdt[celli] > 0.0 && alpha1[celli] > 0.0)
            {
                Sp[celli] -= dgdt[celli]*alpha1[celli];
                Su[celli] += dgdt[celli]*alpha1[celli];
                
            }
            else if (dgdt[celli] < 0.0 && alpha1[celli] < 1.0)
            {
                Sp[celli] += dgdt[celli]*max(1.0 - alpha1[celli],scalar(0));
               
            }
        }

既然sp和dgdt都是体标量场，为什么不直接相加呢？比如用Sp+=dgdt()：这样子形式表达。两者表达形式上直接有什么区别吗

[李东岳]

if (dgdt[celli] > 0.0 && alpha1[celli] > 0.0)
主要是if限制了应用。

[mohui]

@李东岳 哦哦，知道了，那我直接相加一个量也是可以的吧。

[李东岳]

可以啊，Why not？