- Int_t *iPivot = new Int_t[kNn];
- Int_t *indxR = new Int_t[kNn];
- Int_t *indxC = new Int_t[kNn];
- for(Int_t i=0;i<kNn;i++) iPivot[i]=0;
- Int_t iCol=-1,iRow=-1;
- for(Int_t i=0;i<kNn;i++){
- Float_t big=0.;
- for(Int_t j=0;j<kNn;j++){
- if(iPivot[j]!=1){
- for(Int_t k=0;k<kNn;k++){
- if(iPivot[k]==0){
- if(TMath::Abs(mat[j][k])>=big){
- big=TMath::Abs(mat[j][k]);
- iRow=j;
- iCol=k;
- }
+ if(npts<fPolOrder+1){
+ if(fParam) delete [] fParam;
+ fParam=new Float_t[kNn];
+ for(Int_t i=0; i<kNn;i++)fParam[i]=0;
+ }else{
+ Int_t *iPivot = new Int_t[kNn];
+ Int_t *indxR = new Int_t[kNn];
+ Int_t *indxC = new Int_t[kNn];
+ for(Int_t i=0;i<kNn;i++) iPivot[i]=0;
+ Int_t iCol=-1,iRow=-1;
+ for(Int_t i=0;i<kNn;i++){
+ Float_t big=0.;
+ for(Int_t j=0;j<kNn;j++){
+ if(iPivot[j]!=1){
+ for(Int_t k=0;k<kNn;k++){
+ if(iPivot[k]==0){
+ if(TMath::Abs(mat[j][k])>=big){
+ big=TMath::Abs(mat[j][k]);
+ iRow=j;
+ iCol=k;
+ }
+ }