Implementation of the Householder matrix diagonalization algorithm for PDL