High-Performance Design Patterns for Modern Fortran

Listing 1

class(tensor):: u_t, u

real:: nu = 1.0

u_t = nu (.laplacian.u) − (u.dot.(.grad.u))