A slightly different approach for a strictly lower triangular matrix:

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