@@ -236,8 +236,8 @@ tape( 'the function computes a complex quotient (difficult cases)', function tes
236
236
* See section 3.6 in https://arxiv.org/pdf/1210.4539.pdf.
237
237
*
238
238
* ```text
239
- * q[0] : 0011111111100011001100110011001100110011001100110011001100110100
240
- * 0.6: 0011111111100011001100110011001100110011001100110011001100110011
239
+ * real(q) : 0011111111100011001100110011001100110011001100110011001100110100
240
+ * 0.6: 0011111111100011001100110011001100110011001100110011001100110011
241
241
* ```
242
242
*
243
243
* If we add
@@ -246,7 +246,7 @@ tape( 'the function computes a complex quotient (difficult cases)', function tes
246
246
* 0000000000000000000000000000000000000000000000000000000000000001
247
247
* ```
248
248
*
249
- * to `0.6`, we get `q[0] `; thus, the result is 1 bit off.
249
+ * to `0.6`, we get `real( q ) `; thus, the result is 1 bit off.
250
250
*/
251
251
idx = bitdiff ( real ( q ) , 0.6 ) ;
252
252
t . strictEqual ( idx , 61 , 'real component has expected binary representation' ) ;
@@ -322,14 +322,14 @@ tape( 'the function computes a complex quotient (tested against fixtures)', func
322
322
} else {
323
323
delta = abs ( real ( q ) - qre [ i ] ) ;
324
324
tol = EPS * abs ( qre [ i ] ) ;
325
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
325
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
326
326
}
327
327
if ( imag ( q ) === qim [ i ] ) {
328
328
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
329
329
} else {
330
330
delta = abs ( imag ( q ) - qim [ i ] ) ;
331
331
tol = EPS * abs ( qim [ i ] ) ;
332
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
332
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
333
333
}
334
334
}
335
335
t . end ( ) ;
@@ -366,14 +366,14 @@ tape( 'the function computes a complex quotient (different component scales)', f
366
366
} else {
367
367
delta = abs ( real ( q ) - qre [ i ] ) ;
368
368
tol = EPS * abs ( qre [ i ] ) ;
369
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
369
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
370
370
}
371
371
if ( imag ( q ) === qim [ i ] ) {
372
372
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
373
373
} else {
374
374
delta = abs ( imag ( q ) - qim [ i ] ) ;
375
375
tol = EPS * abs ( qim [ i ] ) ;
376
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
376
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
377
377
}
378
378
}
379
379
t . end ( ) ;
@@ -410,14 +410,14 @@ tape( 'the function computes a complex quotient (different component scales)', f
410
410
} else {
411
411
delta = abs ( real ( q ) - qre [ i ] ) ;
412
412
tol = EPS * abs ( qre [ i ] ) ;
413
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
413
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
414
414
}
415
415
if ( imag ( q ) === qim [ i ] ) {
416
416
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
417
417
} else {
418
418
delta = abs ( imag ( q ) - qim [ i ] ) ;
419
419
tol = EPS * abs ( qim [ i ] ) ;
420
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
420
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
421
421
}
422
422
}
423
423
t . end ( ) ;
@@ -454,14 +454,14 @@ tape( 'the function computes a complex quotient (different imaginary component s
454
454
} else {
455
455
delta = abs ( real ( q ) - qre [ i ] ) ;
456
456
tol = EPS * abs ( qre [ i ] ) ;
457
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
457
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
458
458
}
459
459
if ( imag ( q ) === qim [ i ] ) {
460
460
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
461
461
} else {
462
462
delta = abs ( imag ( q ) - qim [ i ] ) ;
463
463
tol = EPS * abs ( qim [ i ] ) ;
464
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
464
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
465
465
}
466
466
}
467
467
t . end ( ) ;
@@ -498,14 +498,14 @@ tape( 'the function computes a complex quotient (real imaginary component scales
498
498
} else {
499
499
delta = abs ( real ( q ) - qre [ i ] ) ;
500
500
tol = EPS * abs ( qre [ i ] ) ;
501
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
501
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
502
502
}
503
503
if ( imag ( q ) === qim [ i ] ) {
504
504
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
505
505
} else {
506
506
delta = abs ( imag ( q ) - qim [ i ] ) ;
507
507
tol = EPS * abs ( qim [ i ] ) ;
508
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
508
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
509
509
}
510
510
}
511
511
t . end ( ) ;
@@ -542,14 +542,14 @@ tape( 'the function computes a complex quotient (large negative imaginary compon
542
542
} else {
543
543
delta = abs ( real ( q ) - qre [ i ] ) ;
544
544
tol = EPS * abs ( qre [ i ] ) ;
545
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
545
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
546
546
}
547
547
if ( imag ( q ) === qim [ i ] ) {
548
548
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
549
549
} else {
550
550
delta = abs ( imag ( q ) - qim [ i ] ) ;
551
551
tol = EPS * abs ( qim [ i ] ) ;
552
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
552
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
553
553
}
554
554
}
555
555
t . end ( ) ;
@@ -586,14 +586,14 @@ tape( 'the function computes a complex quotient (large negative real components)
586
586
} else {
587
587
delta = abs ( real ( q ) - qre [ i ] ) ;
588
588
tol = EPS * abs ( qre [ i ] ) ;
589
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
589
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
590
590
}
591
591
if ( imag ( q ) === qim [ i ] ) {
592
592
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
593
593
} else {
594
594
delta = abs ( imag ( q ) - qim [ i ] ) ;
595
595
tol = EPS * abs ( qim [ i ] ) ;
596
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
596
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
597
597
}
598
598
}
599
599
t . end ( ) ;
@@ -630,14 +630,14 @@ tape( 'the function computes a complex quotient (large positive imaginary compon
630
630
} else {
631
631
delta = abs ( real ( q ) - qre [ i ] ) ;
632
632
tol = EPS * abs ( qre [ i ] ) ;
633
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
633
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
634
634
}
635
635
if ( imag ( q ) === qim [ i ] ) {
636
636
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
637
637
} else {
638
638
delta = abs ( imag ( q ) - qim [ i ] ) ;
639
639
tol = EPS * abs ( qim [ i ] ) ;
640
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
640
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
641
641
}
642
642
}
643
643
t . end ( ) ;
@@ -674,14 +674,14 @@ tape( 'the function computes a complex quotient (large positive real components)
674
674
} else {
675
675
delta = abs ( real ( q ) - qre [ i ] ) ;
676
676
tol = EPS * abs ( qre [ i ] ) ;
677
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
677
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
678
678
}
679
679
if ( imag ( q ) === qim [ i ] ) {
680
680
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
681
681
} else {
682
682
delta = abs ( imag ( q ) - qim [ i ] ) ;
683
683
tol = EPS * abs ( qim [ i ] ) ;
684
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
684
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
685
685
}
686
686
}
687
687
t . end ( ) ;
@@ -718,14 +718,14 @@ tape( 'the function computes a complex quotient (tiny negative imaginary compone
718
718
} else {
719
719
delta = abs ( real ( q ) - qre [ i ] ) ;
720
720
tol = EPS * abs ( qre [ i ] ) ;
721
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
721
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
722
722
}
723
723
if ( imag ( q ) === qim [ i ] ) {
724
724
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
725
725
} else {
726
726
delta = abs ( imag ( q ) - qim [ i ] ) ;
727
727
tol = EPS * abs ( qim [ i ] ) ;
728
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
728
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
729
729
}
730
730
}
731
731
t . end ( ) ;
@@ -762,14 +762,14 @@ tape( 'the function computes a complex quotient (tiny negative real components)'
762
762
} else {
763
763
delta = abs ( real ( q ) - qre [ i ] ) ;
764
764
tol = EPS * abs ( qre [ i ] ) ;
765
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
765
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
766
766
}
767
767
if ( imag ( q ) === qim [ i ] ) {
768
768
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
769
769
} else {
770
770
delta = abs ( imag ( q ) - qim [ i ] ) ;
771
771
tol = EPS * abs ( qim [ i ] ) ;
772
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
772
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
773
773
}
774
774
}
775
775
t . end ( ) ;
@@ -806,14 +806,14 @@ tape( 'the function computes a complex quotient (tiny positive imaginary compone
806
806
} else {
807
807
delta = abs ( real ( q ) - qre [ i ] ) ;
808
808
tol = EPS * abs ( qre [ i ] ) ;
809
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
809
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
810
810
}
811
811
if ( imag ( q ) === qim [ i ] ) {
812
812
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
813
813
} else {
814
814
delta = abs ( imag ( q ) - qim [ i ] ) ;
815
815
tol = EPS * abs ( qim [ i ] ) ;
816
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
816
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
817
817
}
818
818
}
819
819
t . end ( ) ;
@@ -850,14 +850,14 @@ tape( 'the function computes a complex quotient (tiny positive real components)'
850
850
} else {
851
851
delta = abs ( real ( q ) - qre [ i ] ) ;
852
852
tol = EPS * abs ( qre [ i ] ) ;
853
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + q [ 0 ] + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
853
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. real: ' + real ( q ) + '. expected: ' + qre [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
854
854
}
855
855
if ( imag ( q ) === qim [ i ] ) {
856
856
t . strictEqual ( imag ( q ) , qim [ i ] , 'returns expected imaginary component' ) ;
857
857
} else {
858
858
delta = abs ( imag ( q ) - qim [ i ] ) ;
859
859
tol = EPS * abs ( qim [ i ] ) ;
860
- t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + q [ 1 ] + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
860
+ t . ok ( delta <= tol , 'within tolerance. x: ' + re1 [ i ] + ' + ' + im1 [ i ] + 'i. y: ' + re2 [ i ] + ' + ' + im2 [ i ] + 'i. imag: ' + imag ( q ) + '. expected: ' + qim [ i ] + '. delta: ' + delta + '. tol: ' + tol + '.' ) ;
861
861
}
862
862
}
863
863
t . end ( ) ;
0 commit comments