problems with initialisation of piv in symmetric solver