Commit 124a7e07 authored by Claus Kleinwort's avatar Claus Kleinwort
Browse files

Fix for abort on singular QL decomposition of constraints matrix

git-svn-id: http://svnsrv.desy.de/public/MillepedeII/trunk@153 3547b9b0-65b8-46d3-b95d-921b3f43af62
parent 047f0ab4
......@@ -460,10 +460,10 @@ END SUBROUTINE qlpssq
!> Get eigenvalues.
!!
!! Get smallest and largest eingenvalue of L.
!! Get smallest and largest |eigenvalue| of L.
!!
!! \param [out] emin smallest eigenvalue
!! \param [out] emax largest eigenvalue
!! \param [out] emin eigenvalue with smallest absolute value
!! \param [out] emax eigenvalue with largest absolute value
!!
SUBROUTINE qlgete(emin,emax)
USE mpqldec
......
......@@ -52,7 +52,7 @@
!! 1. Download the software package from the DESY \c svn server to
!! \a target directory, e.g.:
!!
!! svn checkout http://svnsrv.desy.de/public/MillepedeII/tags/V04-03-05 target
!! svn checkout http://svnsrv.desy.de/public/MillepedeII/tags/V04-03-06 target
!!
!! 2. Create **Pede** executable (in \a target directory):
!!
......@@ -1318,11 +1318,11 @@ SUBROUTINE feasma
CALL qldec(matConstraintsT)
! check eignevalues of L
CALL qlgete(evmin,evmax)
PRINT *, ' largest eigenvalue of L: ', evmax
PRINT *, ' smallest eigenvalue of L: ', evmin
IF (evmax == 0.0_mpd) THEN
PRINT *, ' largest |eigenvalue| of L: ', evmax
PRINT *, ' smallest |eigenvalue| of L: ', evmin
IF (evmin == 0.0_mpd) THEN
CALL peend(27,'Aborted, singular QL decomposition of constraints matrix')
STOP 'FEASMA: stopping due singular QL decomposition of constraints matrix'
STOP 'FEASMA: stopping due to singular QL decomposition of constraints matrix'
END IF
END IF
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment