trace.f90 Source File

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