------- Comment #2 from eli dot friedman at gmail dot com 2010-07-08 03:08 ------- I'm afraid I mis-reduced the issue; try the following (which reproduces on both trunk r161941 and 4.5.0):
typedef struct LLSModel{ double variance[32]; int indep_count; }LLSModel; void av_solve_lls(LLSModel *m, double (*factor)[33], int count, int min_order){ int i,j,k; for(i=0; i<count; i++){ for(j=i; j<count; j++){ double sum= factor[i][j]; for(k=i-1; k>=0; k--) sum -= factor[i][k]; factor[j][i] = sum; } } for(j=count-1; j>=min_order; j--){ for(i=j; i>=0; i--) factor[j][i]= 1 / factor[i][i]; m->variance[j]= factor[0][0]; for(i=0; i<=j; i++) m->variance[j] += factor[j][i]; } } -- eli dot friedman at gmail dot com changed: What |Removed |Added ---------------------------------------------------------------------------- Status|RESOLVED |UNCONFIRMED Resolution|DUPLICATE | http://gcc.gnu.org/bugzilla/show_bug.cgi?id=44848