fixed the initialize x, added template for inserting D matrix

the compressive sensing (sparse regression) test example works with brgn of the smoothed L1 form

debug cs1 with intermediate output

add compressive sensing test example cs1.c

tao-brgn add feature smoothed L1 norm algorithm, after meeting with Todd

brgn modified with approximation of L1 regularizer

Merge branch 'denera/tao-brgn-matlab-segfault-fix' [PR #1184]

* denera/tao-brgn-matlab-segfault-fix:
MATLAB Leastsquares examples segfault fix for BRGN.
Trust region initialization was producing a s

Merge branch 'denera/tao-brgn-matlab-segfault-fix' [PR #1184]

* denera/tao-brgn-matlab-segfault-fix:
MATLAB Leastsquares examples segfault fix for BRGN.
Trust region initialization was producing a segfault on MatMult with the Jacobian.

Evaluating the Jacobian in MATLAB using a SEQAIJ PETSc type, but then converting it to SEQDENSE before using it in Tao operations resolved the problem.

show more ...

adding the matlab test to the automated test suite

TaoCreate_XYZ() docstring updates

Change CHKERRQ(ierr) after PetscInitialize() calls in examples to if (ierr) return ierr;

Since the PETSc initialization may not have been completed in the failed PetscInitialize() it is not
propper

Change CHKERRQ(ierr) after PetscInitialize() calls in examples to if (ierr) return ierr;

Since the PETSc initialization may not have been completed in the failed PetscInitialize() it is not
propper to call the PETSc error handling routines

Commit-type: bug-fix, testing-fix, example
Reported-by: jed Brown <jed@jedbrown.org>

show more ...

Fix MatSeqAIJFromMatlab() it assumed that the inode structure could change
only if the number of nonzeros changed in the matrix but, in fact, the struture
of the matrix could change while the number

Fix MatSeqAIJFromMatlab() it assumed that the inode structure could change
only if the number of nonzeros changed in the matrix but, in fact, the struture
of the matrix could change while the number of nonzeros remained the same.

Fixed memory leak in example and missing CHKERRQ(ierr);

Commit-type: bug-fix, testing-fix
Reported-by: Alp Dener <adener@anl.gov>

show more ...

matlab LS problem residual jacobian API updated to latest

Using dense Jacobians instead of SeqAIJ in matlab LS problems

Remove outdated requires: from top of several examples

Commit-type: testing-fix

Fixes applied per Lisandro's comments on commit 737f463:
- Order of INCREF and destroy operations fixed, and other Tao setter functions reviewed for correctness.
- Function name changed from `TaoSetR

Fixes applied per Lisandro's comments on commit 737f463:
- Order of INCREF and destroy operations fixed, and other Tao setter functions reviewed for correctness.
- Function name changed from `TaoSetResidualJacobianRoutine()` to `TaoSetJacobianResidualRoutine()` to match the pattern of other Jacobian setters.

show more ...

Regularization changed to x-x_old, user hook added into optimization loops, eliminated explicit MatTransposeMatMult option in Hessian computation

Added support for testing bounded regularized gauss-newton methods
by including Jacobian information.

Added the gauss-newton algorithm wrapper

Pounders updated to use the new LS residual interface

Changed VecScatterCreate() to take a communicator and return scatter as other first class PETSc objects do.
Changed old VecScatterCreate() to VecScatterCreateWithData()

Completed basic refactorizati

Changed VecScatterCreate() to take a communicator and return scatter as other first class PETSc objects do.
Changed old VecScatterCreate() to VecScatterCreateWithData()

Completed basic refactorization of VecScatter code, ready for new implementations to be added.

(There is still some code duplication between the VecScatter MPI1 and MPI3 code but it can be fixed at anytime without requiring changes to the VecScatter API.

Part 5 of a refactorization of VecScatter to be treated as a first class PETSc object

Commit-type: style-fix

show more ...

Tao: Added alternative output for chwirut2 test.

Obtained on an optimized build on FreeBSD in single precision.

another alt file to account for Inf function on first iteration

alt files for archs with optimized BLAS impls

underscored removed from BLAS routine name strings

test problems outputs updated