| line |
true |
false |
branch |
|
101
|
734 |
734 |
for(jj=0; jj
|
|
|
13 |
13 |
for(jj=0; jj
|
|
102
|
1573 |
734 |
for(i=0; i
|
|
|
39 |
13 |
for(i=0; i
|
|
104
|
2517 |
1573 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j) |
|
|
78 |
39 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j) |
|
108
|
76675 |
734 |
for(kk=0; kk
|
|
|
4147 |
13 |
for(kk=0; kk
|
|
109
|
160441 |
76675 |
for(i=0; i
|
|
|
12441 |
4147 |
for(i=0; i
|
|
111
|
251298 |
160441 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j){ |
|
|
24882 |
12441 |
for(j=__MAX__(jj, i); j<__MIN__(jj+bsize, m); ++j){ |
|
113
|
7986678 |
251298 |
for(k=kk; k<__MIN__(kk+bsize, n); ++k){ |
|
|
795678 |
24882 |
for(k=kk; k<__MIN__(kk+bsize, n); ++k){ |
|
124
|
1573 |
734 |
for(i=0; i
|
|
|
39 |
13 |
for(i=0; i
|
|
125
|
944 |
1573 |
for(j=0; j
|
|
|
39 |
39 |
for(j=0; j
|
|
151
|
40169 |
10080 |
for(j=0; j
|
|
|
21126 |
5312 |
for(j=0; j
|
|
154
|
40071 |
98 |
d=FABS(d); |
|
|
21022 |
104 |
d=FABS(d); |
|
155
|
74 |
40095 |
if(d
|
|
|
76 |
21050 |
if(d
|
|
166
|
422478 |
40169 |
for(i=0; i
|
|
|
145458 |
21126 |
for(i=0; i
|
|
189
|
0 |
0 |
for(j=0; j
|
|
|
0 |
0 |
for(j=0; j
|
|
192
|
0 |
0 |
d=FABS(d); |
|
|
0 |
0 |
d=FABS(d); |
|
193
|
0 |
0 |
if(d
|
|
|
0 |
0 |
if(d
|
|
205
|
0 |
0 |
for(i=0; i
|
|
|
0 |
0 |
for(i=0; i
|
|
267
|
0 |
2 |
if(!buf){ |
|
|
0 |
2 |
if(!buf){ |
|
283
|
4 |
2 |
for(j=0; j
|
|
|
4 |
2 |
for(j=0; j
|
|
284
|
4 |
0 |
temp=eps*FABS(p[j]); |
|
|
4 |
0 |
temp=eps*FABS(p[j]); |
|
285
|
0 |
4 |
if(temp==zero) temp=eps; |
|
|
0 |
4 |
if(temp==zero) temp=eps; |
|
295
|
20 |
2 |
for(i=0; i
|
|
|
20 |
2 |
for(i=0; i
|
|
298
|
4 |
2 |
for(j=0; j
|
|
|
4 |
2 |
for(j=0; j
|
|
299
|
4 |
0 |
temp=FABS(p[j]); |
|
|
4 |
0 |
temp=FABS(p[j]); |
|
300
|
0 |
4 |
if(temp==zero) temp=one; |
|
|
0 |
4 |
if(temp==zero) temp=one; |
|
302
|
40 |
4 |
for(i=0; i
|
|
|
40 |
4 |
for(i=0; i
|
|
306
|
20 |
2 |
for(i=0; i
|
|
|
20 |
2 |
for(i=0; i
|
|
308
|
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])) |
|
309
|
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])); |
|
311
|
6 |
14 |
if(temp>epsmch && temp
|
|
|
6 |
0 |
if(temp>epsmch && temp
|
|
|
6 |
14 |
if(temp>epsmch && temp
|
|
|
6 |
0 |
if(temp>epsmch && temp
|
|
313
|
0 |
20 |
if(temp>=eps) err[i]=zero; |
|
|
0 |
20 |
if(temp>=eps) err[i]=zero; |
|
442
|
0 |
143 |
if(!buf){ |
|
|
0 |
16 |
if(!buf){ |
|
453
|
770 |
143 |
for(i=0; i
|
|
|
129 |
16 |
for(i=0; i
|
|
456
|
322 |
143 |
for(i=0; i
|
|
|
43 |
16 |
for(i=0; i
|
|
458
|
770 |
322 |
for(j=0; j
|
|
|
129 |
43 |
for(j=0; j
|
|
459
|
422 |
348 |
if((tmp=FABS(a[i*m+j]))>max) |
|
|
497 |
273 |
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) |
|
461
|
0 |
322 |
if(max==0.0){ |
|
|
0 |
43 |
if(max==0.0){ |
|
470
|
322 |
143 |
for(j=0; j
|
|
|
43 |
16 |
for(j=0; j
|
|
471
|
224 |
322 |
for(i=0; i
|
|
|
43 |
43 |
for(i=0; i
|
|
473
|
54 |
224 |
for(k=0; k
|
|
|
21 |
43 |
for(k=0; k
|
|
478
|
546 |
322 |
for(i=j; i
|
|
|
86 |
43 |
for(i=j; i
|
|
480
|
278 |
546 |
for(k=0; k
|
|
|
64 |
86 |
for(k=0; k
|
|
483
|
368 |
178 |
if((tmp=work[i]*FABS(sum))>=max){ |
|
|
399 |
147 |
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){ |
|
488
|
76 |
246 |
if(j!=maxi){ |
|
|
15 |
28 |
if(j!=maxi){ |
|
489
|
185 |
76 |
for(k=0; k
|
|
|
41 |
15 |
for(k=0; k
|
|
497
|
0 |
322 |
if(a[j*m+j]==0.0) |
|
|
2 |
41 |
if(a[j*m+j]==0.0) |
|
499
|
179 |
143 |
if(j!=m-1){ |
|
|
27 |
16 |
if(j!=m-1){ |
|
501
|
224 |
179 |
for(i=j+1; i
|
|
|
43 |
27 |
for(i=j+1; i
|
|
509
|
322 |
143 |
for(l=0; l
|
|
|
43 |
16 |
for(l=0; l
|
|
510
|
770 |
322 |
for(i=0; i
|
|
|
129 |
43 |
for(i=0; i
|
|
513
|
770 |
322 |
for(i=k=0; i
|
|
|
129 |
43 |
for(i=k=0; i
|
|
517
|
224 |
546 |
if(k!=0) |
|
|
43 |
86 |
if(k!=0) |
|
518
|
278 |
224 |
for(j=k-1; j
|
|
|
64 |
43 |
for(j=k-1; j
|
|
521
|
322 |
224 |
if(sum!=0.0) |
|
|
43 |
43 |
if(sum!=0.0) |
|
526
|
770 |
322 |
for(i=m-1; i>=0; --i){ |
|
|
129 |
43 |
for(i=m-1; i>=0; --i){ |
|
528
|
610 |
770 |
for(j=i+1; j
|
|
|
149 |
129 |
for(j=i+1; j
|
|
533
|
770 |
322 |
for(i=0; i
|
|
|
129 |
43 |
for(i=0; i
|
|
579
|
0 |
143 |
if(!rnk) return 0; |
|
|
0 |
16 |
if(!rnk) return 0; |
|
585
|
770 |
143 |
for(i=0; i
|
|
|
129 |
16 |
for(i=0; i
|
|
624
|
0 |
0 |
if((hx=(LM_REAL *)malloc(n*sizeof(LM_REAL)))==NULL){ |
|
|
0 |
0 |
if((hx=(LM_REAL *)malloc(n*sizeof(LM_REAL)))==NULL){ |
|
632
|
0 |
0 |
for(i=0, tmp=0.0; i
|
|
|
0 |
0 |
for(i=0, tmp=0.0; i
|
|
636
|
0 |
0 |
for(i=0, SSerr=SStot=0.0; i
|
|
|
0 |
0 |
for(i=0, SSerr=SStot=0.0; i
|
|
654
|
11 |
0 |
if(!lb || !ub) return 1; |
|
|
0 |
11 |
if(!lb || !ub) return 1; |
|
|
5 |
0 |
if(!lb || !ub) return 1; |
|
|
0 |
5 |
if(!lb || !ub) return 1; |
|
656
|
42 |
11 |
for(i=0; i
|
|
|
20 |
5 |
for(i=0; i
|
|
657
|
0 |
42 |
if(lb[i]>ub[i]) return 0; |
|
|
0 |
20 |
if(lb[i]>ub[i]) return 0; |
|
726
|
159958 |
0 |
if(x){ |
|
|
72713 |
0 |
if(x){ |
|
727
|
445769 |
159958 |
for(i=blockn-1; i>0; i-=blocksize){ |
|
|
17934 |
72713 |
for(i=blockn-1; i>0; i-=blocksize){ |
|
745
|
158823 |
1135 |
if(i
|
|
|
72713 |
0 |
if(i
|
|
762
|
0 |
0 |
for(i=blockn-1; i>0; i-=blocksize){ |
|
|
0 |
0 |
for(i=blockn-1; i>0; i-=blocksize){ |
|
780
|
0 |
0 |
if(i
|
|
|
0 |
0 |
if(i
|