@@ -30,24 +30,15 @@ var EPS = require( '@stdlib/constants/float64/eps' );
3030var SMALLEST_NORMAL = require ( '@stdlib/constants/float64/smallest-normal' ) ;
3131var PI = require ( '@stdlib/constants/float64/pi' ) ;
3232var kernelBetainc = require ( './../lib' ) ;
33- var fixtures = require ( './fixtures/cpp/output.json' ) ;
34- var smallAB = require ( './fixtures/cpp/small_ab.json' ) ;
35- var mediumAB = require ( './fixtures/cpp/medium_ab.json' ) ;
36- var largeRangeAB = require ( './fixtures/cpp/large_range_ab.json' ) ;
37- var intAB = require ( './fixtures/cpp/int_ab.json' ) ;
38- var largeAsym = require ( './fixtures/cpp/large_asym.json' ) ;
3933
4034
4135// FIXTURES //
4236
43- var lowerRegularized = fixtures . lower_regularized ;
44- var lowerUnregularized = fixtures . lower_unregularized ;
45- var upperRegularized = fixtures . upper_regularized ;
46- var upperUnregularized = fixtures . upper_unregularized ;
47- var derivativeRegularized = fixtures . derivative_regularized ;
48- var x = fixtures . x ;
49- var a = fixtures . a ;
50- var b = fixtures . b ;
37+ var smallRange = require ( './fixtures/cpp/small_range.json' ) ;
38+ var mediumRange = require ( './fixtures/cpp/medium_range.json' ) ;
39+ var largeRangeRange = require ( './fixtures/cpp/large_range.json' ) ;
40+ var intValues = require ( './fixtures/cpp/int_values.json' ) ;
41+ var largeAsym = require ( './fixtures/cpp/large_asym.json' ) ;
5142
5243
5344// TESTS //
@@ -549,132 +540,6 @@ tape( 'the function handles tiny parameter values', function test( t ) {
549540 t . end ( ) ;
550541} ) ;
551542
552- tape ( 'the function evaluates the lower regularized incomplete beta function' , function test ( t ) {
553- var expected ;
554- var delta ;
555- var tol ;
556- var i ;
557- var y ;
558-
559- expected = lowerRegularized ;
560- for ( i = 0 ; i < x . length ; i ++ ) {
561- y = kernelBetainc ( x [ i ] , a [ i ] , b [ i ] , true , false ) [ 0 ] ;
562- if ( y === expected [ i ] ) {
563- t . strictEqual ( y , expected [ i ] , 'x: ' + x [ i ] + ', y: ' + y + '. a: ' + a [ i ] + '. b: ' + b [ i ] + ', expected: ' + expected [ i ] ) ;
564- } else {
565- delta = abs ( y - expected [ i ] ) ;
566- tol = 50.0 * EPS * abs ( expected [ i ] ) ;
567- t . ok ( delta <= tol , 'within tolerance. x: ' + x [ i ] + '. a: ' + a [ i ] + '. b: ' + b [ i ] + '. y: ' + y + '. E: ' + expected [ i ] + '. Δ: ' + delta + '. tol: ' + tol ) ;
568- }
569- }
570- t . end ( ) ;
571- } ) ;
572-
573- tape ( 'the function evaluates the lower unregularized incomplete beta function' , function test ( t ) {
574- var expected ;
575- var delta ;
576- var tol ;
577- var i ;
578- var y ;
579-
580- expected = lowerUnregularized ;
581- for ( i = 0 ; i < x . length ; i ++ ) {
582- y = kernelBetainc ( x [ i ] , a [ i ] , b [ i ] , false , false ) [ 0 ] ;
583- if ( y === expected [ i ] ) {
584- t . strictEqual ( y , expected [ i ] , 'x: ' + x [ i ] + ', y: ' + y + '. a: ' + a [ i ] + '. b: ' + b [ i ] + ', expected: ' + expected [ i ] ) ;
585- } else {
586- delta = abs ( y - expected [ i ] ) ;
587- tol = 90.0 * EPS * abs ( expected [ i ] ) ;
588- t . ok ( delta <= tol , 'within tolerance. x: ' + x [ i ] + '. a: ' + a [ i ] + '. b: ' + b [ i ] + '. y: ' + y + '. E: ' + expected [ i ] + '. Δ: ' + delta + '. tol: ' + tol ) ;
589- }
590- }
591- t . end ( ) ;
592- } ) ;
593-
594- tape ( 'the function evaluates the upper regularized incomplete beta function' , function test ( t ) {
595- var expected ;
596- var delta ;
597- var tol ;
598- var i ;
599- var y ;
600-
601- expected = upperRegularized ;
602- for ( i = 0 ; i < x . length ; i ++ ) {
603- y = kernelBetainc ( x [ i ] , a [ i ] , b [ i ] , true , true ) [ 0 ] ;
604- if ( y === expected [ i ] ) {
605- t . strictEqual ( y , expected [ i ] , 'x: ' + x [ i ] + ', y: ' + y + '. a: ' + a [ i ] + '. b: ' + b [ i ] + ', expected: ' + expected [ i ] ) ;
606- } else {
607- delta = abs ( y - expected [ i ] ) ;
608- tol = 55.0 * EPS * abs ( expected [ i ] ) ;
609- t . ok ( delta <= tol , 'within tolerance. x: ' + x [ i ] + '. a: ' + a [ i ] + '. b: ' + b [ i ] + '. y: ' + y + '. E: ' + expected [ i ] + '. Δ: ' + delta + '. tol: ' + tol ) ;
610- }
611- }
612- t . end ( ) ;
613- } ) ;
614-
615- tape ( 'the function evaluates the upper unregularized incomplete beta function' , function test ( t ) {
616- var expected ;
617- var delta ;
618- var tol ;
619- var i ;
620- var y ;
621-
622- expected = upperUnregularized ;
623- for ( i = 0 ; i < x . length ; i ++ ) {
624- y = kernelBetainc ( x [ i ] , a [ i ] , b [ i ] , false , true ) [ 0 ] ;
625- if ( y === expected [ i ] ) {
626- t . strictEqual ( y , expected [ i ] , 'x: ' + x [ i ] + ', y: ' + y + '. a: ' + a [ i ] + '. b: ' + b [ i ] + ', expected: ' + expected [ i ] ) ;
627- } else {
628- delta = abs ( y - expected [ i ] ) ;
629- tol = 70.0 * EPS * abs ( expected [ i ] ) ;
630- t . ok ( delta <= tol , 'within tolerance. x: ' + x [ i ] + '. a: ' + a [ i ] + '. b: ' + b [ i ] + '. y: ' + y + '. E: ' + expected [ i ] + '. Δ: ' + delta + '. tol: ' + tol ) ;
631- }
632- }
633- t . end ( ) ;
634- } ) ;
635-
636- tape ( 'the function evaluates the lower derivative regularized incomplete beta function' , function test ( t ) {
637- var expected ;
638- var delta ;
639- var tol ;
640- var i ;
641- var y ;
642-
643- expected = derivativeRegularized ;
644- for ( i = 0 ; i < x . length ; i ++ ) {
645- y = kernelBetainc ( x [ i ] , a [ i ] , b [ i ] , true , false ) [ 1 ] ;
646- if ( y === expected [ i ] ) {
647- t . strictEqual ( y , expected [ i ] , 'x: ' + x [ i ] + ', y: ' + y + '. a: ' + a [ i ] + '. b: ' + b [ i ] + ', expected: ' + expected [ i ] ) ;
648- } else {
649- delta = abs ( y - expected [ i ] ) ;
650- tol = 650.0 * EPS * abs ( expected [ i ] ) ;
651- t . ok ( delta <= tol , 'within tolerance. x: ' + x [ i ] + '. a: ' + a [ i ] + '. b: ' + b [ i ] + '. y: ' + y + '. E: ' + expected [ i ] + '. Δ: ' + delta + '. tol: ' + tol ) ;
652- }
653- }
654- t . end ( ) ;
655- } ) ;
656-
657- tape ( 'the function evaluates the upper derivative regularized incomplete beta function' , function test ( t ) {
658- var expected ;
659- var delta ;
660- var tol ;
661- var i ;
662- var y ;
663-
664- expected = derivativeRegularized ;
665- for ( i = 0 ; i < x . length ; i ++ ) {
666- y = kernelBetainc ( x [ i ] , a [ i ] , b [ i ] , true , true ) [ 1 ] ;
667- if ( y === expected [ i ] ) {
668- t . strictEqual ( y , expected [ i ] , 'x: ' + x [ i ] + ', y: ' + y + '. a: ' + a [ i ] + '. b: ' + b [ i ] + ', expected: ' + expected [ i ] ) ;
669- } else {
670- delta = abs ( y - expected [ i ] ) ;
671- tol = 650.0 * EPS * abs ( expected [ i ] ) ;
672- t . ok ( delta <= tol , 'within tolerance. x: ' + x [ i ] + '. a: ' + a [ i ] + '. b: ' + b [ i ] + '. y: ' + y + '. E: ' + expected [ i ] + '. Δ: ' + delta + '. tol: ' + tol ) ;
673- }
674- }
675- t . end ( ) ;
676- } ) ;
677-
678543tape ( 'the function evaluates the lower regularized incomplete beta function for small a,b in [0,10]' , function test ( t ) {
679544 var expected ;
680545 var delta ;
@@ -685,10 +550,10 @@ tape( 'the function evaluates the lower regularized incomplete beta function for
685550 var i ;
686551 var y ;
687552
688- expected = smallAB . lower_regularized ;
689- sv = smallAB . x ;
690- sa = smallAB . a ;
691- sb = smallAB . b ;
553+ expected = smallRange . lower_regularized ;
554+ sv = smallRange . x ;
555+ sa = smallRange . a ;
556+ sb = smallRange . b ;
692557 for ( i = 0 ; i < sv . length ; i ++ ) {
693558 y = kernelBetainc ( sv [ i ] , sa [ i ] , sb [ i ] , true , false ) [ 0 ] ;
694559 if ( y === expected [ i ] ) {
@@ -712,10 +577,10 @@ tape( 'the function evaluates the upper regularized incomplete beta function for
712577 var i ;
713578 var y ;
714579
715- expected = smallAB . upper_regularized ;
716- sv = smallAB . x ;
717- sa = smallAB . a ;
718- sb = smallAB . b ;
580+ expected = smallRange . upper_regularized ;
581+ sv = smallRange . x ;
582+ sa = smallRange . a ;
583+ sb = smallRange . b ;
719584 for ( i = 0 ; i < sv . length ; i ++ ) {
720585 y = kernelBetainc ( sv [ i ] , sa [ i ] , sb [ i ] , true , true ) [ 0 ] ;
721586 if ( y === expected [ i ] ) {
@@ -739,10 +604,10 @@ tape( 'the function evaluates the lower regularized incomplete beta function for
739604 var i ;
740605 var y ;
741606
742- expected = mediumAB . lower_regularized ;
743- mv = mediumAB . x ;
744- ma = mediumAB . a ;
745- mb = mediumAB . b ;
607+ expected = mediumRange . lower_regularized ;
608+ mv = mediumRange . x ;
609+ ma = mediumRange . a ;
610+ mb = mediumRange . b ;
746611 for ( i = 0 ; i < mv . length ; i ++ ) {
747612 y = kernelBetainc ( mv [ i ] , ma [ i ] , mb [ i ] , true , false ) [ 0 ] ;
748613 if ( y === expected [ i ] ) {
@@ -766,10 +631,10 @@ tape( 'the function evaluates the upper regularized incomplete beta function for
766631 var i ;
767632 var y ;
768633
769- expected = mediumAB . upper_regularized ;
770- mv = mediumAB . x ;
771- ma = mediumAB . a ;
772- mb = mediumAB . b ;
634+ expected = mediumRange . upper_regularized ;
635+ mv = mediumRange . x ;
636+ ma = mediumRange . a ;
637+ mb = mediumRange . b ;
773638 for ( i = 0 ; i < mv . length ; i ++ ) {
774639 y = kernelBetainc ( mv [ i ] , ma [ i ] , mb [ i ] , true , true ) [ 0 ] ;
775640 if ( y === expected [ i ] ) {
@@ -793,10 +658,10 @@ tape( 'the function evaluates the lower regularized incomplete beta function for
793658 var i ;
794659 var y ;
795660
796- expected = largeRangeAB . lower_regularized ;
797- lv = largeRangeAB . x ;
798- la = largeRangeAB . a ;
799- lb = largeRangeAB . b ;
661+ expected = largeRangeRange . lower_regularized ;
662+ lv = largeRangeRange . x ;
663+ la = largeRangeRange . a ;
664+ lb = largeRangeRange . b ;
800665 for ( i = 0 ; i < lv . length ; i ++ ) {
801666 y = kernelBetainc ( lv [ i ] , la [ i ] , lb [ i ] , true , false ) [ 0 ] ;
802667 if ( y === expected [ i ] ) {
@@ -820,10 +685,10 @@ tape( 'the function evaluates the upper regularized incomplete beta function for
820685 var i ;
821686 var y ;
822687
823- expected = largeRangeAB . upper_regularized ;
824- lv = largeRangeAB . x ;
825- la = largeRangeAB . a ;
826- lb = largeRangeAB . b ;
688+ expected = largeRangeRange . upper_regularized ;
689+ lv = largeRangeRange . x ;
690+ la = largeRangeRange . a ;
691+ lb = largeRangeRange . b ;
827692 for ( i = 0 ; i < lv . length ; i ++ ) {
828693 y = kernelBetainc ( lv [ i ] , la [ i ] , lb [ i ] , true , true ) [ 0 ] ;
829694 if ( y === expected [ i ] ) {
@@ -847,10 +712,10 @@ tape( 'the function evaluates the lower regularized incomplete beta function for
847712 var i ;
848713 var y ;
849714
850- expected = intAB . lower_regularized ;
851- iv = intAB . x ;
852- ia = intAB . a ;
853- ib = intAB . b ;
715+ expected = intValues . lower_regularized ;
716+ iv = intValues . x ;
717+ ia = intValues . a ;
718+ ib = intValues . b ;
854719 for ( i = 0 ; i < iv . length ; i ++ ) {
855720 y = kernelBetainc ( iv [ i ] , ia [ i ] , ib [ i ] , true , false ) [ 0 ] ;
856721 if ( y === expected [ i ] ) {
@@ -874,10 +739,10 @@ tape( 'the function evaluates the upper regularized incomplete beta function for
874739 var i ;
875740 var y ;
876741
877- expected = intAB . upper_regularized ;
878- iv = intAB . x ;
879- ia = intAB . a ;
880- ib = intAB . b ;
742+ expected = intValues . upper_regularized ;
743+ iv = intValues . x ;
744+ ia = intValues . a ;
745+ ib = intValues . b ;
881746 for ( i = 0 ; i < iv . length ; i ++ ) {
882747 y = kernelBetainc ( iv [ i ] , ia [ i ] , ib [ i ] , true , true ) [ 0 ] ;
883748 if ( y === expected [ i ] ) {
0 commit comments