Dear Professor,
I am so sorry for my frequent problems.
According to your previous advice, I have switched my code to a single varf. However, the results from the code with a single varf are still wrong.
I find that the results between the two codes with solve and single are same when the time is 1, while the results of iteration over time are wrong. Therefore, I think that the codes of the variational form are correctly written. I guess that the code about iteration may be wrong. However, the code about iteration is totally same between the codes with solve and varf. I suspect there is an intrinsic error in my code with varf, which I do not understand. Therefore, please give me some more advice.
The two codes with “solve” and “varf” are attached to this post, respectively.
Sequential_2D_Solve.edp (2.2 KB)
Sequential_2D_Varf1.edp (2.2 KB)
Best,
Liu