Ignore:
Timestamp:
05/13/11 22:54:42 (9 years ago)
Author:
aivar
Message:

Changed GCI_gauss_jordan to be GCI_solve and GCI_invert. The current default CGI_solve uses Gaussian Elimination. Added inversion method. Added Lower/Upper Decomposition code for both solve and invert, but it is commented out.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/projects/slim-curve/src/main/c/EcfSingle.c

    r7675 r7706  
    614614 
    615615 
    616 int GCI_gauss_jordan(float **a, int n, float *b); 
    617  
    618616int GCI_marquardt_step_instr(float xincr, float y[], 
    619617                                        int ndata, int fit_start, int fit_end, 
     
    680678        } 
    681679 
    682         /* Matrix solution; GCI_gauss_jordan solves Ax=b rather than AX=B */ 
    683         if (GCI_gauss_jordan(covar, *pmfit, dparam) != 0) { 
     680        /* Matrix solution; GCI_solve solves Ax=b rather than AX=B */ 
     681        if (GCI_solve(covar, *pmfit, dparam) != 0) { 
    684682                return -1; 
    685683        } 
    686          
     684 
     685        //TODO need to make sure covar gets inverted.  Previously the Gauss 
     686        // Jordan solution would invert covar as a side effect. 
    687687 
    688688        /* Once converged, evaluate covariance matrix */ 
Note: See TracChangeset for help on using the changeset viewer.