[Solved] VSPRINC in VUMAT

Hi All,
I implemented a damage model in VUMAT subrooutine, and I would like to obtain tha maximum principal strain using the utility routine call vsprinc( nblock, s, eigVal, ndir, nshr ).
In my code, I assigned the following utility routine :
call VSPRINC(nblock,STRAIN,eigVal,NDIR,NSHR)
Prin_strain=ZERO
do p=1,3
Prin_strain=max(Prin_strain,eigVal(i,p)) ! Prin_strain is the principal strain, and STRAIN is the strain tensor
enddo

Can anyone please respond to me if this definition if correct?
if not, how can one assign the utility routine to obtain the maximum strain ?

Thank you in advance.

I dont use the VSPRINC function. I have written my own version (so I can use it in UMATS for other FE solvers) 😎
Best of luck,
Jorgen