Provide state/error handling for linear algebra

doc/specs/index.md

@@ -22,6 +22,7 @@ This is an index/directory of the specifications (specs) for each new module/fea
  - [io](./stdlib_io.html) - Input/output helper & convenience
  - [kinds](./stdlib_kinds.html) - Kind parameters
  - [linalg](./stdlib_linalg.html) - Linear Algebra
+ - [linalg_state_type](./stdlib_linalg_state_type.html) - Linear Algebra state and error handling 
  - [logger](./stdlib_logger.html) - Runtime logging system
  - [math](./stdlib_math.html) - General purpose mathematical functions
  - [optval](./stdlib_optval.html) - Fallback value for optional arguments

doc/specs/stdlib_linalg_state_type.md

@@ -0,0 +1,64 @@
+---
+title: linalg_state_type
+---
+
+# Linear Algebra -- State and Error Handling Module
+
+[TOC]
+
+## Introduction
+
+The `stdlib_linalg_state` module provides a derived type holding information on the  
+state of linear algebra operations, and procedures for expert control of linear algebra workflows. 
+All linear algebra procedures are engineered to support returning an optional `linalg_state_type` 
+variable to holds such information, as a form of expert API. If the user does not require state 
+information, but fatal errors are encountered during the execution of linear algebra routines, the 
+program will undergo a hard stop.
+Instead, if the state argument is present, the program will never stop, but will return detailed error 
+information into the state handler. 
+
+## Derived types provided
+
+<!-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -->
+### The `linalg_state_type` derived type
+
+The `linalg_state_type` is defined as a derived type containing an integer error flag, and 
+fixed-size character strings to store an error message and the location of the error state change. 
+Fixed-size string storage was chosen to facilitate the compiler's memory allocation and ultimately 
+ensure maximum computational performance.  
+
+A similarly named generic interface, `linalg_state_type`, is provided to allow the developer to 
+create diagnostic messages and raise error flags easily. The call starts with an error flag or 
+the location of the event, and is followed by an arbitrary list of `integer`, `real`, `complex` or 
+`character` variables. Numeric variables may be provided as either scalars or rank-1 (array) inputs. 
+
+#### Type-bound procedures
+
+The following convenience type-bound procedures are provided: 
+- `print()` returns an allocatable character string containing state location, message, and error flag; 
+- `print_message()` returns an allocatable character string containing the state message; 
+- `ok()` returns a `logical` flag that is `.true.` in case of successful state (`flag==LINALG_SUCCESS`);
+- `error()` returns a `logical` flag that is `.true.` in case of error state (`flag/=LINALG_SUCCESS`).
+
+#### Status
+
+Experimental
+
+#### Example
+
+```fortran
+{!example/linalg/example_state1.f90!}
+```
+
+## Error flags provided
+
+The module provides the following state flags: 
+- `LINALG_SUCCESS`: Successful execution
+- `LINALG_VALUE_ERROR`: Numerical errors (such as infinity, not-a-number, range bounds) are encountered.
+- `LINALG_ERROR`: Linear Algebra errors are encountered, such as: non-converging iterations, impossible operations, etc.
+- `LINALG_INTERNAL_ERROR`: Provided as a developer safeguard for internal errors that should never occur.
+
+## Comparison operators provided
+
+The module provides overloaded comparison operators for all comparisons of a `linalg_state_type` variable 
+with an integer error flag: `<`, `<=`, `==`, `>=`, `>`, `/=`.

example/linalg/example_state1.f90

@@ -0,0 +1,18 @@
+program example_state1
+  use stdlib_linalg_state
+  implicit none
+  type(linalg_state_type) :: err
+
+  ! Create a state flag
+  err = linalg_state_type(LINALG_VALUE_ERROR,'just an example with scalar ',&
+                          'integer=',1,'real=',2.0,'complex=',(3.0,1.0),'and array ',[1,2,3],'inputs')
+
+  ! Print flag
+  print *, err%print()
+
+  ! Check success
+  print *, 'Check error: ',err%error()
+  print *, 'Check flag : ',err /= LINALG_SUCCESS
+
+
+end program example_state1

src/CMakeLists.txt

@@ -24,6 +24,7 @@ set(fppFiles
     stdlib_linalg_outer_product.fypp
     stdlib_linalg_kronecker.fypp
     stdlib_linalg_cross_product.fypp
+    stdlib_linalg_state.fypp 
     stdlib_optval.fypp
     stdlib_selection.fypp
     stdlib_sorting.fypp

src/stdlib_linalg_constants.fypp

@@ -7,7 +7,6 @@ module stdlib_linalg_constants
      public
 
 
-
      ! Integer size support for ILP64 builds should be done here
      integer, parameter :: ilp = int32
      private            :: int32, int64