Prolog "sincos"
/**
* Warranty & Liability
* To the extent permitted by applicable law and unless explicitly
* otherwise agreed upon, XLOG Technologies AG makes no warranties
* regarding the provided information. XLOG Technologies AG assumes
* no liability that any problems might be solved with the information
* provided by XLOG Technologies AG.
*
* Rights & License
* All industrial property rights regarding the information - copyright
* and patent rights in particular - are the sole property of XLOG
* Technologies AG. If the company was not the originator of some
* excerpts, XLOG Technologies AG has at least obtained the right to
* reproduce, change and translate the information.
*
* Reproduction is restricted to the whole unaltered document. Reproduction
* of the information is only allowed for non-commercial uses. Selling,
* giving away or letting of the execution of the library is prohibited.
* The library can be distributed as part of your applications and libraries
* for execution provided this comment remains unchanged.
*
* Restrictions
* Only to be distributed with programs that add significant and primary
* functionality to the library. Not to be distributed with additional
* software intended to replace any components of the library.
*
* Trademarks
* Jekejeke is a registered trademark of XLOG Technologies AG.
*/
/**
* Source of test cases are the following standards:
* - Prolog General Core ISO/IEC 13211-1
* - Draft Technical Corrigendum 2, WG17, Ulrich Neumerkel
* <a href="https://www.complang.tuwien.ac.at/ulrich/iso-prolog/dtc2">www.complang.tuwien.ac.at/ulrich/iso-prolog/dtc2</a>
*/
:- multifile runner_case/5.
:- discontiguous(runner_case/5).
:- multifile(runner_pred/5).
:- discontiguous(runner_pred/5).
:- multifile(runner_file/3).
runner_file(arithmetic, sincos, 'ISO 9.3 trigo').
/****************************************************************/
/* Trigonometric Operations */
/****************************************************************/
/* X ** Y */
runner_pred(**,3, arithmetic, sincos, 'ISO 9.3.1.4').
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, ISO 1') :-
125.0 is 5 ** 3.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, ISO 2') :-
-125.0 is (-5) ** 3.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, ISO 3') :-
0.2 is 5 ** -1.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, ISO 4') :-
catch(_ is 77 ** _, error(E, _), true),
E == instantiation_error.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, ISO 6') :-
125.0 is 5 ** 3.0.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, ISO 7') :-
1.0 is 0.0 ** 0.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, XLOG 1') :-
0.0025252525252525255 =:= 396** -1.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, XLOG 2') :-
1.0250068783051207E95 =:= 1889**29.
runner_case(**,3, arithmetic, sincos, 'ISO 9.3.1.4, XLOG 3') :-
1.8678023704280373E-14 =:= -2705** -4.
/* sin(X) */
runner_pred(sin,2, arithmetic, sincos, 'ISO 9.3.2.4').
runner_case(sin,2, arithmetic, sincos, 'ISO 9.3.2.4, ISO 1') :-
0.0 is sin(0.0).
runner_case(sin,2, arithmetic, sincos, 'ISO 9.3.2.4, ISO 3') :-
0.0 is sin(0).
runner_case(sin,2, arithmetic, sincos, 'ISO 9.3.2.4, ISO 4') :-
1.0 is sin(pi/2.0).
runner_case(sin,2, arithmetic, sincos, 'ISO 9.3.2.4, XLOG 1') :-
0.7457052121767203 is sin(2.3).
runner_case(sin,2, arithmetic, sincos, 'ISO 9.3.2.4, XLOG 2') :-
catch(_ is sin(inf), error(E,_), true),
E == evaluation_error(undefined).
runner_case(sin,2, arithmetic, sincos, 'ISO 9.3.2.4, ISO 5') :-
catch(_ is sin(foobar), error(E, _), true),
E == type_error(evaluable, foobar/0).
runner_case(sin,2, arithmetic, sincos, 'ISO 9.3.2.4, XLOG 3') :-
catch(sin(0.1,_), error(E,_), true),
E == existence_error(procedure, sin/2).
/* cos(X) */
runner_pred(cos,2, arithmetic, sincos, 'ISO 9.3.3.4').
runner_case(cos,2, arithmetic, sincos, 'ISO 9.3.3.4, ISO 1') :-
1.0 is cos(0.0).
runner_case(cos,2, arithmetic, sincos, 'ISO 9.3.3.4, ISO 2') :-
catch(_ is cos(_), error(E, _), true),
E == instantiation_error.
runner_case(cos,2, arithmetic, sincos, 'ISO 9.3.3.4, ISO 3') :-
1.0 is cos(0).
runner_case(cos,2, arithmetic, sincos, 'ISO 9.3.3.4, ISO 5') :-
6.123233995736766e-17 is cos(pi/2.0).
runner_case(cos,2, arithmetic, sincos, 'ISO 9.3.3.4, XLOG 1') :-
X is cos(2.3), (X == -0.6662760212798241; X = -0.666276021279824).
/* tan(X) */
runner_pred(tan,2, arithmetic, sincos, 'Corr.2 9.3.14.4').
runner_case(tan,2, arithmetic, sincos, 'Corr.2 9.3.14.4, XLOG 1') :-
X is tan(pi/6), 1/3 =:= X*X.
runner_case(tan,2, arithmetic, sincos, 'Corr.2 9.3.14.4, XLOG 2') :-
0.9999999999999999 is tan(pi/4).
runner_case(tan,2, arithmetic, sincos, 'Corr.2 9.3.14.4, XLOG 3') :-
catch(_ is tan(foobar), error(E, _), true),
E == type_error(evaluable, foobar/0).
runner_case(tan,2, arithmetic, sincos, 'Corr.2 9.3.14.4, XLOG 4') :-
-1.1192136417341325 is tan(2.3).
/* asin(X) */
runner_pred(asin,2, arithmetic, sincos, 'Corr.2 9.3.11.4').
runner_case(asin,2, arithmetic, sincos, 'Corr.2 9.3.11.4, XLOG 1') :-
0.0 is asin(0).
runner_case(asin,2, arithmetic, sincos, 'Corr.2 9.3.11.4, XLOG 2') :-
pi =:= asin(sqrt(3/4))*3.
runner_case(asin,2, arithmetic, sincos, 'Corr.2 9.3.11.4, XLOG 3') :-
pi/2 =:= asin(1).
runner_case(asin,2, arithmetic, sincos, 'Corr.2 9.3.11.4, XLOG 4') :-
catch(_ is asin(_), error(E, _), true),
E == instantiation_error.
runner_case(asin,2, arithmetic, sincos, 'Corr.2 9.3.11.4, ISO 3') :-
catch(_ is asin(2), error(E, _), true),
E == evaluation_error(undefined).
/* acos(X) */
runner_pred(acos,2, arithmetic, sincos, 'Corr.2 9.3.12.4').
runner_case(acos,2, arithmetic, sincos, 'Corr.2 9.3.12.4, XLOG 1') :-
pi/2 =:= acos(0).
runner_case(acos,2, arithmetic, sincos, 'Corr.2 9.3.12.4, XLOG 2') :-
3.1415926535897936 is acos(sqrt(3/4))*6.
runner_case(acos,2, arithmetic, sincos, 'Corr.2 9.3.12.4, XLOG 3') :-
0.0 is acos(1).
runner_case(acos,2, arithmetic, sincos, 'Corr.2 9.3.12.4, XLOG 4') :-
catch(_ is acos(foobar), error(E, _), true),
E == type_error(evaluable, foobar/0).
runner_case(acos,2, arithmetic, sincos, 'Corr.2 9.3.12.4, ISO 3') :-
catch(_ is acos(1.5), error(E, _), true),
E == evaluation_error(undefined).
/* atan(X) */
runner_pred(atan,2, arithmetic, sincos, 'ISO 9.3.4.4').
runner_case(atan,2, arithmetic, sincos, 'ISO 9.3.4.4, ISO 1') :-
0.0 is atan(0.0).
runner_case(atan,2, arithmetic, sincos, 'ISO 9.3.4.4, ISO 2') :-
pi =:= atan(1.0)*4.
runner_case(atan,2, arithmetic, sincos, 'ISO 9.3.4.4, ISO 3') :-
catch(_ is atan(_), error(E, _), true),
E == instantiation_error.
runner_case(atan,2, arithmetic, sincos, 'ISO 9.3.4.4, ISO 4') :-
0.0 is atan(0).
/* exp(X) */
runner_pred(exp,2, arithmetic, sincos, 'ISO 9.3.5.4').
runner_case(exp,2, arithmetic, sincos, 'ISO 9.3.5.4, ISO 1') :-
1.0 is exp(0.0).
runner_case(exp,2, arithmetic, sincos, 'ISO 9.3.5.4, ISO 2') :-
(2.7182818284590455 is exp(1.0); 2.718281828459045 is exp(1.0)).
runner_case(exp,2, arithmetic, sincos, 'ISO 9.3.5.4, ISO 4') :-
1.0 is exp(0).
runner_case(exp,2, arithmetic, sincos, 'ISO 9.3.5.4, ISO 5') :-
catch(_ is exp(foobar), error(E, _), true),
E == type_error(evaluable, foobar/0).
/* log(X) */
runner_pred(log,2, arithmetic, sincos, 'ISO 9.3.6.4').
runner_case(log,2, arithmetic, sincos, 'ISO 9.3.6.4, ISO 1') :-
0.0 is log(1.0).
runner_case(log,2, arithmetic, sincos, 'ISO 9.3.6.4, ISO 2') :-
1.0 is log(e).
runner_case(log,2, arithmetic, sincos, 'ISO 9.3.6.4, ISO 3') :-
catch(_ is log(_), error(E, _), true),
E == instantiation_error.
runner_case(log,2, arithmetic, sincos, 'ISO 9.3.6.4, ISO 4') :-
catch(_ is log(0), error(E, _), true),
nonvar(E), E=evaluation_error(_).
runner_case(log,2, arithmetic, sincos, 'ISO 9.3.6.4, ISO 6') :-
catch(_ is log(0.0), error(E, _), true),
nonvar(E), E=evaluation_error(_).
runner_case(log,2, arithmetic, sincos, 'ISO 9.3.6.4, ISO 7') :-
catch(_ is log(-1), error(E, _), true),
nonvar(E), E=evaluation_error(undefined).
/* sqrt(X) */
runner_pred(sqrt,2, arithmetic, sincos, 'ISO 9.3.7.4').
runner_case(sqrt,2, arithmetic, sincos, 'ISO 9.3.7.4, ISO 1') :-
0.0 is sqrt(0.0).
runner_case(sqrt,2, arithmetic, sincos, 'ISO 9.3.7.4, ISO 2') :-
1.0 is sqrt(1).
runner_case(sqrt,2, arithmetic, sincos, 'ISO 9.3.7.4, ISO 3') :-
1.1 is sqrt(1.21).
runner_case(sqrt,2, arithmetic, sincos, 'ISO 9.3.7.4, ISO 5') :-
catch(_ is sqrt(-1.0), error(E, _), true),
E == evaluation_error(undefined).
runner_case(sqrt,2, arithmetic, sincos, 'ISO 9.3.7.4, ISO 6') :-
catch(_ is sqrt(foobar), error(E, _), true),
E == type_error(evaluable, foobar/0).
/* pi */
runner_pred(pi,1, arithmetic, sincos, 'Corr.2 9.3.15.4').
runner_case(pi,1, arithmetic, sincos, 'Corr.2 9.3.15.4, ISO 1') :-
3.141592653589793 is pi.
/* atan2 */
runner_pred(atan2,3, arithmetic, sincos, 'Corr.2 9.3.13').
runner_case(atan2,3, arithmetic, sincos, 'Corr.2 9.3.13, ISO 1') :-
0.0 is atan2(1,0)-pi/2.
runner_case(atan2,3, arithmetic, sincos, 'Corr.2 9.3.13, ISO 2') :-
0.0 is atan2(0,-1)-pi.
runner_case(atan2,3, arithmetic, sincos, 'Corr.2 9.3.13, ISO 3') :-
catch(_ is atan2(0,0), error(E, _), true),
E == evaluation_error(undefined).
/* epsilon */
runner_pred(epsilon,1, arithmetic, sincos, 'N208 9.7.3').
runner_case(epsilon,1, arithmetic, sincos, 'N208 9.7.3, XLOG 1') :-
\+ 1+epsilon =:= 1.
runner_case(epsilon,1, arithmetic, sincos, 'N208 9.7.3, XLOG 2') :-
2+epsilon =:= 2.
/* e */
runner_pred(e,1, arithmetic, sincos, 'N208 9.7.2').
runner_case(e,1, arithmetic, sincos, 'N208 9.7.2, XLOG 1') :-
2.718281828459045 is e.