trace.f90 Source File

submodule(specialmatrices_tridiagonal) tridiagonal_trace
   implicit none(type, external)
contains
   module procedure trace_rdp
   tr = sum(A%dv)
   end procedure
end submodule