FEniCS Course - FEniCS Project
FEniCS Course - FEniCS Project
FEniCS Course - FEniCS Project
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Detailed time-stepping algorithm for the heatequationDefine the boundary conditionCompute u 0 as the projection of the given initial valueDefine the forms a and LAssemble the matrix A from the bilinear form at ← ∆twhile t T doAssemble the vector b from the linear form LApply the boundary conditionSolve the linear system AU = b for U and store in u 1t ← t + ∆tu 0 ← u 1 (get ready for next step)end while6 / 13