Skip to content
Snippets Groups Projects
Commit cfcf391d authored by Thomas Toulorge's avatar Thomas Toulorge
Browse files

Fixed problem with initial step length calculation in ALGLIB CG optimizer

parent 5ff81146
No related branches found
No related tags found
No related merge requests found
...@@ -3449,9 +3449,9 @@ INPUT PARAMETERS: ...@@ -3449,9 +3449,9 @@ INPUT PARAMETERS:
A - matrix, array[N,N] A - matrix, array[N,N]
IsUpper - (optional) storage type: IsUpper - (optional) storage type:
* if True, symmetric matrix A is given by its upper * if True, symmetric matrix A is given by its upper
triangle, and the lower triangle isnt used triangle, and the lower triangle isn�t used
* if False, symmetric matrix A is given by its lower * if False, symmetric matrix A is given by its lower
triangle, and the upper triangle isnt used triangle, and the upper triangle isn�t used
* if not given, both lower and upper triangles must be * if not given, both lower and upper triangles must be
filled. filled.
...@@ -3496,9 +3496,9 @@ INPUT PARAMETERS: ...@@ -3496,9 +3496,9 @@ INPUT PARAMETERS:
A - matrix, array[N,N] A - matrix, array[N,N]
IsUpper - (optional) storage type: IsUpper - (optional) storage type:
* if True, symmetric matrix A is given by its upper * if True, symmetric matrix A is given by its upper
triangle, and the lower triangle isnt used triangle, and the lower triangle isn�t used
* if False, symmetric matrix A is given by its lower * if False, symmetric matrix A is given by its lower
triangle, and the upper triangle isnt used triangle, and the upper triangle isn�t used
* if not given, both lower and upper triangles must be * if not given, both lower and upper triangles must be
filled. filled.
...@@ -6702,7 +6702,7 @@ lbl_24: ...@@ -6702,7 +6702,7 @@ lbl_24:
* to the norm of DK * to the norm of DK
*/ */
v = ae_v_dotproduct(&state->dk.ptr.p_double[0], 1, &state->dk.ptr.p_double[0], 1, ae_v_len(0,n-1)); v = ae_v_dotproduct(&state->dk.ptr.p_double[0], 1, &state->dk.ptr.p_double[0], 1, ae_v_len(0,n-1));
state->laststep = ae_sqrt(v, _state); state->laststep = 1.0/ae_sqrt(v, _state);
} }
else else
{ {
...@@ -11553,9 +11553,9 @@ INPUT PARAMETERS: ...@@ -11553,9 +11553,9 @@ INPUT PARAMETERS:
A - matrix, array[N,N] A - matrix, array[N,N]
IsUpper - (optional) storage type: IsUpper - (optional) storage type:
* if True, symmetric matrix A is given by its upper * if True, symmetric matrix A is given by its upper
triangle, and the lower triangle isnt used triangle, and the lower triangle isn�t used
* if False, symmetric matrix A is given by its lower * if False, symmetric matrix A is given by its lower
triangle, and the upper triangle isnt used triangle, and the upper triangle isn�t used
* if not given, both lower and upper triangles must be * if not given, both lower and upper triangles must be
filled. filled.
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment