use PDL::LinearAlgebra::Real ... tricpy(zeros($size,$size), 0, my $sltm = ones($size, $size));