Skip to content

Conversation

angsch
Copy link
Collaborator

@angsch angsch commented Sep 17, 2024

Description

The documentation and the code of GESVJ deviate in two spots with respect to the usage of .GT. and .GE.:

*> The number of rows of the input matrix A. 1/DLAMCH('E') > M >= 0.

IF( DBLE( M )*EPSLN.GE.ONE ) THEN

The PR takes the version in the code and updates the documentation.

*> It is required that CTOL >= ONE, i.e. it is not

ELSE IF( UCTOL .AND. ( WORK( 1 ).LE.ONE ) ) THEN

This deviation is fixed in the code, which now literally permits EPS as threshold.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants