Hi,

According to the derivation for implicit solve and the formula in the comment, it seems to me that

// k3 = rhs3 - dt curl k2
k3 = RHS3;
Curl.AddMult(k2, k3, -dt);

instead? Otherwise the - dt curl k2 term will not be updated to k3.

Please let me know if my understanding is correct. Thanks.