13c827367444ee418f129b2c238299f49d3264554Jarkko Poyry/*------------------------------------------------------------------------- 23c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * drawElements Base Portability Library 33c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * ------------------------------------- 43c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * 53c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * Copyright 2014 The Android Open Source Project 63c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * 73c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * Licensed under the Apache License, Version 2.0 (the "License"); 83c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * you may not use this file except in compliance with the License. 93c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * You may obtain a copy of the License at 103c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * 113c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * http://www.apache.org/licenses/LICENSE-2.0 123c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * 133c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * Unless required by applicable law or agreed to in writing, software 143c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * distributed under the License is distributed on an "AS IS" BASIS, 153c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 163c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * See the License for the specific language governing permissions and 173c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * limitations under the License. 183c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * 193c827367444ee418f129b2c238299f49d3264554Jarkko Poyry *//*! 203c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * \file 213c827367444ee418f129b2c238299f49d3264554Jarkko Poyry * \brief Basic mathematical operations. 223c827367444ee418f129b2c238299f49d3264554Jarkko Poyry *//*--------------------------------------------------------------------*/ 233c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 243c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#include "deMath.h" 253c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 263c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#if (DE_COMPILER == DE_COMPILER_MSC) 273c827367444ee418f129b2c238299f49d3264554Jarkko Poyry# include <float.h> 283c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#endif 293c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 303c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#if (DE_COMPILER == DE_COMPILER_GCC) || (DE_COMPILER == DE_COMPILER_CLANG) 313c827367444ee418f129b2c238299f49d3264554Jarkko Poyry# include <fenv.h> 323c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#endif 333c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 343c827367444ee418f129b2c238299f49d3264554Jarkko PoyrydeRoundingMode deGetRoundingMode (void) 353c827367444ee418f129b2c238299f49d3264554Jarkko Poyry{ 363c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#if (DE_COMPILER == DE_COMPILER_MSC) 373c827367444ee418f129b2c238299f49d3264554Jarkko Poyry unsigned int status = 0; 383c827367444ee418f129b2c238299f49d3264554Jarkko Poyry int ret; 393c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 403c827367444ee418f129b2c238299f49d3264554Jarkko Poyry ret = _controlfp_s(&status, 0, 0); 413c827367444ee418f129b2c238299f49d3264554Jarkko Poyry DE_ASSERT(ret == 0); 423c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 433c827367444ee418f129b2c238299f49d3264554Jarkko Poyry switch (status & _MCW_RC) 443c827367444ee418f129b2c238299f49d3264554Jarkko Poyry { 453c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case _RC_CHOP: return DE_ROUNDINGMODE_TO_ZERO; 463c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case _RC_UP: return DE_ROUNDINGMODE_TO_POSITIVE_INF; 473c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case _RC_DOWN: return DE_ROUNDINGMODE_TO_NEGATIVE_INF; 483c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case _RC_NEAR: return DE_ROUNDINGMODE_TO_NEAREST; 493c827367444ee418f129b2c238299f49d3264554Jarkko Poyry default: return DE_ROUNDINGMODE_LAST; 503c827367444ee418f129b2c238299f49d3264554Jarkko Poyry } 513c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#elif (DE_COMPILER == DE_COMPILER_GCC) || (DE_COMPILER == DE_COMPILER_CLANG) 523c827367444ee418f129b2c238299f49d3264554Jarkko Poyry int mode = fegetround(); 533c827367444ee418f129b2c238299f49d3264554Jarkko Poyry switch (mode) 543c827367444ee418f129b2c238299f49d3264554Jarkko Poyry { 553c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case FE_TOWARDZERO: return DE_ROUNDINGMODE_TO_ZERO; 563c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case FE_UPWARD: return DE_ROUNDINGMODE_TO_POSITIVE_INF; 573c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case FE_DOWNWARD: return DE_ROUNDINGMODE_TO_NEGATIVE_INF; 583c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case FE_TONEAREST: return DE_ROUNDINGMODE_TO_NEAREST; 593c827367444ee418f129b2c238299f49d3264554Jarkko Poyry default: return DE_ROUNDINGMODE_LAST; 603c827367444ee418f129b2c238299f49d3264554Jarkko Poyry } 613c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#else 623c827367444ee418f129b2c238299f49d3264554Jarkko Poyry# error Implement deGetRoundingMode(). 633c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#endif 643c827367444ee418f129b2c238299f49d3264554Jarkko Poyry} 653c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 663c827367444ee418f129b2c238299f49d3264554Jarkko PoyrydeBool deSetRoundingMode (deRoundingMode mode) 673c827367444ee418f129b2c238299f49d3264554Jarkko Poyry{ 683c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#if (DE_COMPILER == DE_COMPILER_MSC) 693c827367444ee418f129b2c238299f49d3264554Jarkko Poyry unsigned int flag = 0; 703c827367444ee418f129b2c238299f49d3264554Jarkko Poyry unsigned int oldState; 713c827367444ee418f129b2c238299f49d3264554Jarkko Poyry int ret; 723c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 733c827367444ee418f129b2c238299f49d3264554Jarkko Poyry switch (mode) 743c827367444ee418f129b2c238299f49d3264554Jarkko Poyry { 753c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_ZERO: flag = _RC_CHOP; break; 763c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_POSITIVE_INF: flag = _RC_UP; break; 773c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_NEGATIVE_INF: flag = _RC_DOWN; break; 783c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_NEAREST: flag = _RC_NEAR; break; 793c827367444ee418f129b2c238299f49d3264554Jarkko Poyry default: 803c827367444ee418f129b2c238299f49d3264554Jarkko Poyry DE_ASSERT(DE_FALSE); 813c827367444ee418f129b2c238299f49d3264554Jarkko Poyry } 823c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 833c827367444ee418f129b2c238299f49d3264554Jarkko Poyry ret = _controlfp_s(&oldState, flag, _MCW_RC); 843c827367444ee418f129b2c238299f49d3264554Jarkko Poyry return ret == 0; 853c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#elif (DE_COMPILER == DE_COMPILER_GCC) || (DE_COMPILER == DE_COMPILER_CLANG) 863c827367444ee418f129b2c238299f49d3264554Jarkko Poyry int flag = 0; 873c827367444ee418f129b2c238299f49d3264554Jarkko Poyry int ret; 883c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 893c827367444ee418f129b2c238299f49d3264554Jarkko Poyry switch (mode) 903c827367444ee418f129b2c238299f49d3264554Jarkko Poyry { 913c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_ZERO: flag = FE_TOWARDZERO; break; 923c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_POSITIVE_INF: flag = FE_UPWARD; break; 933c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_NEGATIVE_INF: flag = FE_DOWNWARD; break; 943c827367444ee418f129b2c238299f49d3264554Jarkko Poyry case DE_ROUNDINGMODE_TO_NEAREST: flag = FE_TONEAREST; break; 953c827367444ee418f129b2c238299f49d3264554Jarkko Poyry default: 963c827367444ee418f129b2c238299f49d3264554Jarkko Poyry DE_ASSERT(DE_FALSE); 973c827367444ee418f129b2c238299f49d3264554Jarkko Poyry } 983c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 993c827367444ee418f129b2c238299f49d3264554Jarkko Poyry ret = fesetround(flag); 1003c827367444ee418f129b2c238299f49d3264554Jarkko Poyry return ret == 0; 1013c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#else 1023c827367444ee418f129b2c238299f49d3264554Jarkko Poyry# error Implement deSetRoundingMode(). 1033c827367444ee418f129b2c238299f49d3264554Jarkko Poyry#endif 1043c827367444ee418f129b2c238299f49d3264554Jarkko Poyry} 1053c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 1063c827367444ee418f129b2c238299f49d3264554Jarkko Poyrydouble deFractExp (double x, int* exponent) 1073c827367444ee418f129b2c238299f49d3264554Jarkko Poyry{ 1083c827367444ee418f129b2c238299f49d3264554Jarkko Poyry if (deIsInf(x)) 1093c827367444ee418f129b2c238299f49d3264554Jarkko Poyry { 1103c827367444ee418f129b2c238299f49d3264554Jarkko Poyry *exponent = 0; 1113c827367444ee418f129b2c238299f49d3264554Jarkko Poyry return x; 1123c827367444ee418f129b2c238299f49d3264554Jarkko Poyry } 1133c827367444ee418f129b2c238299f49d3264554Jarkko Poyry else 1143c827367444ee418f129b2c238299f49d3264554Jarkko Poyry { 1153c827367444ee418f129b2c238299f49d3264554Jarkko Poyry int tmpExp = 0; 1163c827367444ee418f129b2c238299f49d3264554Jarkko Poyry double fract = frexp(x, &tmpExp); 1173c827367444ee418f129b2c238299f49d3264554Jarkko Poyry *exponent = tmpExp - 1; 1183c827367444ee418f129b2c238299f49d3264554Jarkko Poyry return fract * 2.0; 1193c827367444ee418f129b2c238299f49d3264554Jarkko Poyry } 1203c827367444ee418f129b2c238299f49d3264554Jarkko Poyry} 1213c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 1223c827367444ee418f129b2c238299f49d3264554Jarkko Poyry/* We could use frexpf, if available. */ 1233c827367444ee418f129b2c238299f49d3264554Jarkko Poyryfloat deFloatFractExp (float x, int* exponent) 1243c827367444ee418f129b2c238299f49d3264554Jarkko Poyry{ 1253c827367444ee418f129b2c238299f49d3264554Jarkko Poyry return (float)deFractExp(x, exponent); 1263c827367444ee418f129b2c238299f49d3264554Jarkko Poyry} 1273c827367444ee418f129b2c238299f49d3264554Jarkko Poyry 1283c827367444ee418f129b2c238299f49d3264554Jarkko Poyrydouble deRoundEven (double a) 1293c827367444ee418f129b2c238299f49d3264554Jarkko Poyry{ 1303c827367444ee418f129b2c238299f49d3264554Jarkko Poyry double integer; 1313c827367444ee418f129b2c238299f49d3264554Jarkko Poyry double fract = modf(a, &integer); 1323c827367444ee418f129b2c238299f49d3264554Jarkko Poyry if (fabs(fract) == 0.5) 1333c827367444ee418f129b2c238299f49d3264554Jarkko Poyry return 2.0 * deRound(a / 2.0); 1343c827367444ee418f129b2c238299f49d3264554Jarkko Poyry return deRound(a); 1353c827367444ee418f129b2c238299f49d3264554Jarkko Poyry} 136