anchor.diff:         solver_parameters = solver_parameters_in
       ELSE
          solver_parameters = FixedSolverParameters_t()
  +    END IF
  +
  +    !! Print out parameters
  +    IF (solver_parameters%be_verbose) THEN
  +       CALL WriteHeader("Linear Solver")
  +       CALL EnterSubLog
  +       CALL WriteElement(key="Method", text_value_in="Cholesky Decomposition")
  +       CALL PrintFixedSolverParameters(solver_parameters)
       END IF
   
       !! First get the local matrix in a dense recommendation for quick lookup
  
positive.code:     CALL UnpackCholesky(values_per_column_l, index_l, values_l, LMat)
  
      !! Cleanup
      DEALLOCATE(pivot_vector)
      DEALLOCATE(diag)
      DEALLOCATE(dense_a)
  
negative.code:     CALL UnpackCholesky(values_per_column_l, index_l, values_l, LMat)
  
      !! Cleanup
      IF (solver_parameters%be_verbose) THEN
         CALL PrintMatrixInformation(LMat)
         CALL ExitSubLog
      END IF
      DEALLOCATE(dense_a)
      DEALLOCATE(values_per_column_l)
      DEALLOCATE(index_l)
  