| line |
true |
false |
branch |
|
103
|
720 |
720 |
for(jj=0; jj
|
|
|
13 |
13 |
for(jj=0; jj
|
|
104
|
1548 |
720 |
for(i=0; i
|
|
|
39 |
13 |
for(i=0; i
|
|
106
|
2484 |
1548 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j) |
|
|
78 |
39 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j) |
|
110
|
71450 |
720 |
for(kk=0; kk
|
|
|
4147 |
13 |
for(kk=0; kk
|
|
111
|
150087 |
71450 |
for(i=0; i
|
|
|
12441 |
4147 |
for(i=0; i
|
|
113
|
235911 |
150087 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j){ |
|
|
24882 |
12441 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j){ |
|
115
|
7494678 |
235911 |
for(k=kk; k<__MIN__(kk+bsize, n); ++k){ |
|
|
795678 |
24882 |
for(k=kk; k<__MIN__(kk+bsize, n); ++k){ |
|
126
|
1548 |
720 |
for(i=0; i
|
|
|
39 |
13 |
for(i=0; i
|
|
127
|
936 |
1548 |
for(j=0; j
|
|
|
39 |
39 |
for(j=0; j
|
|
153
|
40167 |
10079 |
for(j=0; j
|
|
|
21126 |
5312 |
for(j=0; j
|
|
156
|
40070 |
97 |
d=FABS(d); |
|
|
21022 |
104 |
d=FABS(d); |
|
157
|
74 |
40093 |
if(d
|
|
|
76 |
21050 |
if(d
|
|
168
|
402478 |
40167 |
for(i=0; i
|
|
|
145458 |
21126 |
for(i=0; i
|
|
191
|
0 |
0 |
for(j=0; j
|
|
|
0 |
0 |
for(j=0; j
|
|
194
|
0 |
0 |
d=FABS(d); |
|
|
0 |
0 |
d=FABS(d); |
|
195
|
0 |
0 |
if(d
|
|
|
0 |
0 |
if(d
|
|
207
|
0 |
0 |
for(i=0; i
|
|
|
0 |
0 |
for(i=0; i
|
|
269
|
0 |
2 |
if(!buf){ |
|
|
0 |
2 |
if(!buf){ |
|
285
|
4 |
2 |
for(j=0; j
|
|
|
4 |
2 |
for(j=0; j
|
|
286
|
4 |
0 |
temp=eps*FABS(p[j]); |
|
|
4 |
0 |
temp=eps*FABS(p[j]); |
|
287
|
0 |
4 |
if(temp==zero) temp=eps; |
|
|
0 |
4 |
if(temp==zero) temp=eps; |
|
297
|
20 |
2 |
for(i=0; i
|
|
|
20 |
2 |
for(i=0; i
|
|
300
|
4 |
2 |
for(j=0; j
|
|
|
4 |
2 |
for(j=0; j
|
|
301
|
4 |
0 |
temp=FABS(p[j]); |
|
|
4 |
0 |
temp=FABS(p[j]); |
|
302
|
0 |
4 |
if(temp==zero) temp=one; |
|
|
0 |
4 |
if(temp==zero) temp=one; |
|
304
|
40 |
4 |
for(i=0; i
|
|
|
40 |
4 |
for(i=0; i
|
|
308
|
20 |
2 |
for(i=0; i
|
|
|
20 |
2 |
for(i=0; i
|
|
310
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
6 |
14 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
6 |
14 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
|
20 |
0 |
if(fvec[i]!=zero && fvecp[i]!=zero && FABS(fvecp[i]-fvec[i])>=epsf*FABS(fvec[i])) |
|
311
|
14 |
6 |
temp=eps*FABS((fvecp[i]-fvec[i])/eps - err[i])/(FABS(fvec[i])+FABS(fvecp[i])); |
|
|
20 |
0 |
temp=eps*FABS((fvecp[i]-fvec[i])/eps - err[i])/(FABS(fvec[i])+FABS(fvecp[i])); |
|
|
20 |
0 |
temp=eps*FABS((fvecp[i]-fvec[i])/eps - err[i])/(FABS(fvec[i])+FABS(fvecp[i])); |
|
|
10 |
10 |
temp=eps*FABS((fvecp[i]-fvec[i])/eps - err[i])/(FABS(fvec[i])+FABS(fvecp[i])); |
|
|
20 |
0 |
temp=eps*FABS((fvecp[i]-fvec[i])/eps - err[i])/(FABS(fvec[i])+FABS(fvecp[i])); |
|
|
20 |
0 |
temp=eps*FABS((fvecp[i]-fvec[i])/eps - err[i])/(FABS(fvec[i])+FABS(fvecp[i])); |
|
313
|
6 |
14 |
if(temp>epsmch && temp
|
|
|
6 |
0 |
if(temp>epsmch && temp
|
|
|
6 |
14 |
if(temp>epsmch && temp
|
|
|
6 |
0 |
if(temp>epsmch && temp
|
|
315
|
0 |
20 |
if(temp>=eps) err[i]=zero; |
|
|
0 |
20 |
if(temp>=eps) err[i]=zero; |
|
444
|
0 |
144 |
if(!buf){ |
|
|
0 |
16 |
if(!buf){ |
|
455
|
779 |
144 |
for(i=0; i
|
|
|
129 |
16 |
for(i=0; i
|
|
458
|
325 |
144 |
for(i=0; i
|
|
|
43 |
16 |
for(i=0; i
|
|
460
|
779 |
325 |
for(j=0; j
|
|
|
129 |
43 |
for(j=0; j
|
|
461
|
427 |
352 |
if((tmp=FABS(a[i*m+j]))>max) |
|
|
503 |
276 |
if((tmp=FABS(a[i*m+j]))>max) |
|
|
81 |
48 |
if((tmp=FABS(a[i*m+j]))>max) |
|
|
59 |
70 |
if((tmp=FABS(a[i*m+j]))>max) |
|
463
|
0 |
325 |
if(max==0.0){ |
|
|
0 |
43 |
if(max==0.0){ |
|
472
|
325 |
144 |
for(j=0; j
|
|
|
43 |
16 |
for(j=0; j
|
|
473
|
227 |
325 |
for(i=0; i
|
|
|
43 |
43 |
for(i=0; i
|
|
475
|
55 |
227 |
for(k=0; k
|
|
|
21 |
43 |
for(k=0; k
|
|
480
|
552 |
325 |
for(i=j; i
|
|
|
86 |
43 |
for(i=j; i
|
|
482
|
282 |
552 |
for(k=0; k
|
|
|
64 |
86 |
for(k=0; k
|
|
485
|
371 |
181 |
if((tmp=work[i]*FABS(sum))>=max){ |
|
|
403 |
149 |
if((tmp=work[i]*FABS(sum))>=max){ |
|
|
61 |
25 |
if((tmp=work[i]*FABS(sum))>=max){ |
|
|
59 |
27 |
if((tmp=work[i]*FABS(sum))>=max){ |
|
490
|
77 |
248 |
if(j!=maxi){ |
|
|
16 |
27 |
if(j!=maxi){ |
|
491
|
188 |
77 |
for(k=0; k
|
|
|
44 |
16 |
for(k=0; k
|
|
499
|
0 |
325 |
if(a[j*m+j]==0.0) |
|
|
2 |
41 |
if(a[j*m+j]==0.0) |
|
501
|
181 |
144 |
if(j!=m-1){ |
|
|
27 |
16 |
if(j!=m-1){ |
|
503
|
227 |
181 |
for(i=j+1; i
|
|
|
43 |
27 |
for(i=j+1; i
|
|
511
|
325 |
144 |
for(l=0; l
|
|
|
43 |
16 |
for(l=0; l
|
|
512
|
779 |
325 |
for(i=0; i
|
|
|
129 |
43 |
for(i=0; i
|
|
515
|
779 |
325 |
for(i=k=0; i
|
|
|
129 |
43 |
for(i=k=0; i
|
|
519
|
227 |
552 |
if(k!=0) |
|
|
43 |
86 |
if(k!=0) |
|
520
|
282 |
227 |
for(j=k-1; j
|
|
|
64 |
43 |
for(j=k-1; j
|
|
523
|
325 |
227 |
if(sum!=0.0) |
|
|
43 |
43 |
if(sum!=0.0) |
|
528
|
779 |
325 |
for(i=m-1; i>=0; --i){ |
|
|
129 |
43 |
for(i=m-1; i>=0; --i){ |
|
530
|
619 |
779 |
for(j=i+1; j
|
|
|
149 |
129 |
for(j=i+1; j
|
|
535
|
779 |
325 |
for(i=0; i
|
|
|
129 |
43 |
for(i=0; i
|
|
581
|
0 |
144 |
if(!rnk) return 0; |
|
|
0 |
16 |
if(!rnk) return 0; |
|
587
|
779 |
144 |
for(i=0; i
|
|
|
129 |
16 |
for(i=0; i
|
|
626
|
0 |
0 |
if((hx=(LM_REAL *)malloc(n*sizeof(LM_REAL)))==NULL){ |
|
|
0 |
0 |
if((hx=(LM_REAL *)malloc(n*sizeof(LM_REAL)))==NULL){ |
|
634
|
0 |
0 |
for(i=n, tmp=0.0; i-->0; ) |
|
|
0 |
0 |
for(i=n, tmp=0.0; i-->0; ) |
|
638
|
0 |
0 |
if(x) |
|
|
0 |
0 |
if(x) |
|
639
|
0 |
0 |
for(i=n, SSerr=SStot=0.0; i-->0; ){ |
|
|
0 |
0 |
for(i=n, SSerr=SStot=0.0; i-->0; ){ |
|
647
|
0 |
0 |
for(i=n, SSerr=SStot=0.0; i-->0; ){ |
|
|
0 |
0 |
for(i=n, SSerr=SStot=0.0; i-->0; ){ |
|
665
|
12 |
0 |
if(!lb || !ub) return 1; |
|
|
0 |
12 |
if(!lb || !ub) return 1; |
|
|
5 |
0 |
if(!lb || !ub) return 1; |
|
|
0 |
5 |
if(!lb || !ub) return 1; |
|
667
|
45 |
12 |
for(i=0; i
|
|
|
20 |
5 |
for(i=0; i
|
|
668
|
0 |
45 |
if(lb[i]>ub[i]) return 0; |
|
|
0 |
20 |
if(lb[i]>ub[i]) return 0; |