Line 142: if (Ld(j) > tol) bool(2)=1. I compared the code with 'https://github.com/tcew/nodal-dg.git', I think it should be bool(3)=1. Finally, in the process of learning DG method, your code has given me a lot of help. Thank you much.