// filename:c2011-F-10-3-4-ex.c // original examples and/or notes: // (c) ISO/IEC JTC1 SC22 WG14 N1570, April 12, 2011 // http://www.open-std.org/jtc1/sc22/wg14/www/docs/n1570.pdf // C2011 Annex F.10.3.4 The frexp functions // compile and output mechanism: // (c) Ogawa Kiyoshi, kaizen@gifu-u.ac.jp, January.04, 2014 // compile errors and/or wornings: // 1 (c) Apple LLVM version 4.2 (clang-425.0.27) (based on LLVM 3.2svn) // Target: x86_64-apple-darwin11.4.2 //Thread model: posix // (c) LLVM 2003-2009 University of Illinois at Urbana-Champaign. // 2 gcc-4.9 (GCC) 4.9.0 20131229 (experimental) // Copyright (C) 2013 Free Software Foundation, Inc. #include #include #define MAX 5 double frexp_a(double value, int *exp){ *exp = (value == 0) ? 0 : (int)(1 + logb(value)); printf("logb(%f)=%d, scalbn(value, -(%d))=%f \n", value, (int)(1 + logb(value)), *exp, scalbn(value, -(*exp))); return scalbn(value, -(*exp)); } int main(void) { int exp,exp_a; double value[]={3.14159,-2.5,0,FP_INFINITE,FP_NAN },ret,ret_a; for(int i=0;i