.file "__vexp.S"
#include "libm.h"
RO_DATA
.align 16
TBL:
.word 0x3ff00000,0x00000000
.word 0x00000000,0x00000000
.word 0x3ff00b1a,0xfa5abcbf
.word 0xbc84f6b2,0xa7609f71
.word 0x3ff0163d,0xa9fb3335
.word 0x3c9b6129,0x9ab8cdb7
.word 0x3ff02168,0x143b0281
.word 0xbc82bf31,0x0fc54eb6
.word 0x3ff02c9a,0x3e778061
.word 0xbc719083,0x535b085d
.word 0x3ff037d4,0x2e11bbcc
.word 0x3c656811,0xeeade11a
.word 0x3ff04315,0xe86e7f85
.word 0xbc90a31c,0x1977c96e
.word 0x3ff04e5f,0x72f654b1
.word 0x3c84c379,0x3aa0d08c
.word 0x3ff059b0,0xd3158574
.word 0x3c8d73e2,0xa475b465
.word 0x3ff0650a,0x0e3c1f89
.word 0xbc95cb7b,0x5799c396
.word 0x3ff0706b,0x29ddf6de
.word 0xbc8c91df,0xe2b13c26
.word 0x3ff07bd4,0x2b72a836
.word 0x3c832334,0x54458700
.word 0x3ff08745,0x18759bc8
.word 0x3c6186be,0x4bb284ff
.word 0x3ff092bd,0xf66607e0
.word 0xbc968063,0x800a3fd1
.word 0x3ff09e3e,0xcac6f383
.word 0x3c914878,0x18316136
.word 0x3ff0a9c7,0x9b1f3919
.word 0x3c85d16c,0x873d1d38
.word 0x3ff0b558,0x6cf9890f
.word 0x3c98a62e,0x4adc610a
.word 0x3ff0c0f1,0x45e46c85
.word 0x3c94f989,0x06d21cef
.word 0x3ff0cc92,0x2b7247f7
.word 0x3c901edc,0x16e24f71
.word 0x3ff0d83b,0x23395dec
.word 0xbc9bc14d,0xe43f316a
.word 0x3ff0e3ec,0x32d3d1a2
.word 0x3c403a17,0x27c57b53
.word 0x3ff0efa5,0x5fdfa9c5
.word 0xbc949db9,0xbc54021b
.word 0x3ff0fb66,0xaffed31b
.word 0xbc6b9bed,0xc44ebd7b
.word 0x3ff10730,0x28d7233e
.word 0x3c8d46eb,0x1692fdd5
.word 0x3ff11301,0xd0125b51
.word 0xbc96c510,0x39449b3a
.word 0x3ff11edb,0xab5e2ab6
.word 0xbc9ca454,0xf703fb72
.word 0x3ff12abd,0xc06c31cc
.word 0xbc51b514,0xb36ca5c7
.word 0x3ff136a8,0x14f204ab
.word 0xbc67108f,0xba48dcf0
.word 0x3ff1429a,0xaea92de0
.word 0xbc932fbf,0x9af1369e
.word 0x3ff14e95,0x934f312e
.word 0xbc8b91e8,0x39bf44ab
.word 0x3ff15a98,0xc8a58e51
.word 0x3c82406a,0xb9eeab0a
.word 0x3ff166a4,0x5471c3c2
.word 0x3c58f23b,0x82ea1a32
.word 0x3ff172b8,0x3c7d517b
.word 0xbc819041,0xb9d78a76
.word 0x3ff17ed4,0x8695bbc0
.word 0x3c709e3f,0xe2ac5a64
.word 0x3ff18af9,0x388c8dea
.word 0xbc911023,0xd1970f6c
.word 0x3ff19726,0x58375d2f
.word 0x3c94aadd,0x85f17e08
.word 0x3ff1a35b,0xeb6fcb75
.word 0x3c8e5b4c,0x7b4968e4
.word 0x3ff1af99,0xf8138a1c
.word 0x3c97bf85,0xa4b69280
.word 0x3ff1bbe0,0x84045cd4
.word 0xbc995386,0x352ef607
.word 0x3ff1c82f,0x95281c6b
.word 0x3c900977,0x8010f8c9
.word 0x3ff1d487,0x3168b9aa
.word 0x3c9e016e,0x00a2643c
.word 0x3ff1e0e7,0x5eb44027
.word 0xbc96fdd8,0x088cb6de
.word 0x3ff1ed50,0x22fcd91d
.word 0xbc91df98,0x027bb78c
.word 0x3ff1f9c1,0x8438ce4d
.word 0xbc9bf524,0xa097af5c
.word 0x3ff2063b,0x88628cd6
.word 0x3c8dc775,0x814a8494
.word 0x3ff212be,0x3578a819
.word 0x3c93592d,0x2cfcaac9
.word 0x3ff21f49,0x917ddc96
.word 0x3c82a97e,0x9494a5ee
.word 0x3ff22bdd,0xa27912d1
.word 0x3c8d34fb,0x5577d69e
.word 0x3ff2387a,0x6e756238
.word 0x3c99b07e,0xb6c70573
.word 0x3ff2451f,0xfb82140a
.word 0x3c8acfcc,0x911ca996
.word 0x3ff251ce,0x4fb2a63f
.word 0x3c8ac155,0xbef4f4a4
.word 0x3ff25e85,0x711ece75
.word 0x3c93e1a2,0x4ac31b2c
.word 0x3ff26b45,0x65e27cdd
.word 0x3c82bd33,0x9940e9d9
.word 0x3ff2780e,0x341ddf29
.word 0x3c9e067c,0x05f9e76c
.word 0x3ff284df,0xe1f56381
.word 0xbc9a4c3a,0x8c3f0d7e
.word 0x3ff291ba,0x7591bb70
.word 0xbc82cc72,0x28401cbc
.word 0x3ff29e9d,0xf51fdee1
.word 0x3c8612e8,0xafad1255
.word 0x3ff2ab8a,0x66d10f13
.word 0xbc995743,0x191690a7
.word 0x3ff2b87f,0xd0dad990
.word 0xbc410adc,0xd6381aa4
.word 0x3ff2c57e,0x39771b2f
.word 0xbc950145,0xa6eb5124
.word 0x3ff2d285,0xa6e4030b
.word 0x3c900247,0x54db41d5
.word 0x3ff2df96,0x1f641589
.word 0x3c9d16cf,0xfbbce198
.word 0x3ff2ecaf,0xa93e2f56
.word 0x3c71ca0f,0x45d52383
.word 0x3ff2f9d2,0x4abd886b
.word 0xbc653c55,0x532bda93
.word 0x3ff306fe,0x0a31b715
.word 0x3c86f46a,0xd23182e4
.word 0x3ff31432,0xedeeb2fd
.word 0x3c8959a3,0xf3f3fcd0
.word 0x3ff32170,0xfc4cd831
.word 0x3c8a9ce7,0x8e18047c
.word 0x3ff32eb8,0x3ba8ea32
.word 0xbc9c45e8,0x3cb4f318
.word 0x3ff33c08,0xb26416ff
.word 0x3c932721,0x843659a6
.word 0x3ff34962,0x66e3fa2d
.word 0xbc835a75,0x930881a4
.word 0x3ff356c5,0x5f929ff1
.word 0xbc8b5cee,0x5c4e4628
.word 0x3ff36431,0xa2de883b
.word 0xbc8c3144,0xa06cb85e
.word 0x3ff371a7,0x373aa9cb
.word 0xbc963aea,0xbf42eae2
.word 0x3ff37f26,0x231e754a
.word 0xbc99f5ca,0x9eceb23c
.word 0x3ff38cae,0x6d05d866
.word 0xbc9e958d,0x3c9904bd
.word 0x3ff39a40,0x1b7140ef
.word 0xbc99a9a5,0xfc8e2934
.word 0x3ff3a7db,0x34e59ff7
.word 0xbc75e436,0xd661f5e3
.word 0x3ff3b57f,0xbfec6cf4
.word 0x3c954c66,0xe26fff18
.word 0x3ff3c32d,0xc313a8e5
.word 0xbc9efff8,0x375d29c3
.word 0x3ff3d0e5,0x44ede173
.word 0x3c7fe8d0,0x8c284c71
.word 0x3ff3dea6,0x4c123422
.word 0x3c8ada09,0x11f09ebc
.word 0x3ff3ec70,0xdf1c5175
.word 0xbc8af663,0x7b8c9bca
.word 0x3ff3fa45,0x04ac801c
.word 0xbc97d023,0xf956f9f3
.word 0x3ff40822,0xc367a024
.word 0x3c8bddf8,0xb6f4d048
.word 0x3ff4160a,0x21f72e2a
.word 0xbc5ef369,0x1c309278
.word 0x3ff423fb,0x2709468a
.word 0xbc98462d,0xc0b314dd
.word 0x3ff431f5,0xd950a897
.word 0xbc81c7dd,0xe35f7998
.word 0x3ff43ffa,0x3f84b9d4
.word 0x3c8880be,0x9704c002
.word 0x3ff44e08,0x6061892d
.word 0x3c489b7a,0x04ef80d0
.word 0x3ff45c20,0x42a7d232
.word 0xbc686419,0x82fb1f8e
.word 0x3ff46a41,0xed1d0057
.word 0x3c9c944b,0xd1648a76
.word 0x3ff4786d,0x668b3237
.word 0xbc9c20f0,0xed445733
.word 0x3ff486a2,0xb5c13cd0
.word 0x3c73c1a3,0xb69062f0
.word 0x3ff494e1,0xe192aed2
.word 0xbc83b289,0x5e499ea0
.word 0x3ff4a32a,0xf0d7d3de
.word 0x3c99cb62,0xf3d1be56
.word 0x3ff4b17d,0xea6db7d7
.word 0xbc8125b8,0x7f2897f0
.word 0x3ff4bfda,0xd5362a27
.word 0x3c7d4397,0xafec42e2
.word 0x3ff4ce41,0xb817c114
.word 0x3c905e29,0x690abd5d
.word 0x3ff4dcb2,0x99fddd0d
.word 0x3c98ecdb,0xbc6a7833
.word 0x3ff4eb2d,0x81d8abff
.word 0xbc95257d,0x2e5d7a52
.word 0x3ff4f9b2,0x769d2ca7
.word 0xbc94b309,0xd25957e3
.word 0x3ff50841,0x7f4531ee
.word 0x3c7a249b,0x49b7465f
.word 0x3ff516da,0xa2cf6642
.word 0xbc8f7685,0x69bd93ee
.word 0x3ff5257d,0xe83f4eef
.word 0xbc7c998d,0x43efef71
.word 0x3ff5342b,0x569d4f82
.word 0xbc807abe,0x1db13cac
.word 0x3ff542e2,0xf4f6ad27
.word 0x3c87926d,0x192d5f7e
.word 0x3ff551a4,0xca5d920f
.word 0xbc8d689c,0xefede59a
.word 0x3ff56070,0xdde910d2
.word 0xbc90fb6e,0x168eebf0
.word 0x3ff56f47,0x36b527da
.word 0x3c99bb2c,0x011d93ad
.word 0x3ff57e27,0xdbe2c4cf
.word 0xbc90b98c,0x8a57b9c4
.word 0x3ff58d12,0xd497c7fd
.word 0x3c8295e1,0x5b9a1de8
.word 0x3ff59c08,0x27ff07cc
.word 0xbc97e2ce,0xe467e60f
.word 0x3ff5ab07,0xdd485429
.word 0x3c96324c,0x054647ad
.word 0x3ff5ba11,0xfba87a03
.word 0xbc9b77a1,0x4c233e1a
.word 0x3ff5c926,0x8a5946b7
.word 0x3c3c4b1b,0x816986a2
.word 0x3ff5d845,0x90998b93
.word 0xbc9cd6a7,0xa8b45642
.word 0x3ff5e76f,0x15ad2148
.word 0x3c9ba6f9,0x3080e65e
.word 0x3ff5f6a3,0x20dceb71
.word 0xbc89eadd,0xe3cdcf92
.word 0x3ff605e1,0xb976dc09
.word 0xbc93e242,0x9b56de47
.word 0x3ff6152a,0xe6cdf6f4
.word 0x3c9e4b3e,0x4ab84c27
.word 0x3ff6247e,0xb03a5585
.word 0xbc9383c1,0x7e40b497
.word 0x3ff633dd,0x1d1929fd
.word 0x3c984710,0xbeb964e5
.word 0x3ff64346,0x34ccc320
.word 0xbc8c483c,0x759d8932
.word 0x3ff652b9,0xfebc8fb7
.word 0xbc9ae3d5,0xc9a73e08
.word 0x3ff66238,0x82552225
.word 0xbc9bb609,0x87591c34
.word 0x3ff671c1,0xc70833f6
.word 0xbc8e8732,0x586c6134
.word 0x3ff68155,0xd44ca973
.word 0x3c6038ae,0x44f73e65
.word 0x3ff690f4,0xb19e9538
.word 0x3c8804bd,0x9aeb445c
.word 0x3ff6a09e,0x667f3bcd
.word 0xbc9bdd34,0x13b26456
.word 0x3ff6b052,0xfa75173e
.word 0x3c7a38f5,0x2c9a9d0e
.word 0x3ff6c012,0x750bdabf
.word 0xbc728956,0x67ff0b0d
.word 0x3ff6cfdc,0xddd47645
.word 0x3c9c7aa9,0xb6f17309
.word 0x3ff6dfb2,0x3c651a2f
.word 0xbc6bbe3a,0x683c88ab
.word 0x3ff6ef92,0x98593ae5
.word 0xbc90b974,0x9e1ac8b2
.word 0x3ff6ff7d,0xf9519484
.word 0xbc883c0f,0x25860ef6
.word 0x3ff70f74,0x66f42e87
.word 0x3c59d644,0xd45aa65f
.word 0x3ff71f75,0xe8ec5f74
.word 0xbc816e47,0x86887a99
.word 0x3ff72f82,0x86ead08a
.word 0xbc920aa0,0x2cd62c72
.word 0x3ff73f9a,0x48a58174
.word 0xbc90a8d9,0x6c65d53c
.word 0x3ff74fbd,0x35d7cbfd
.word 0x3c9047fd,0x618a6e1c
.word 0x3ff75feb,0x564267c9
.word 0xbc902459,0x57316dd3
.word 0x3ff77024,0xb1ab6e09
.word 0x3c9b7877,0x169147f8
.word 0x3ff78069,0x4fde5d3f
.word 0x3c9866b8,0x0a02162c
.word 0x3ff790b9,0x38ac1cf6
.word 0x3c9349a8,0x62aadd3e
.word 0x3ff7a114,0x73eb0187
.word 0xbc841577,0xee04992f
.word 0x3ff7b17b,0x0976cfdb
.word 0xbc9bebb5,0x8468dc88
.word 0x3ff7c1ed,0x0130c132
.word 0x3c9f124c,0xd1164dd6
.word 0x3ff7d26a,0x62ff86f0
.word 0x3c91bddb,0xfb72b8b4
.word 0x3ff7e2f3,0x36cf4e62
.word 0x3c705d02,0xba15797e
.word 0x3ff7f387,0x8491c491
.word 0xbc807f11,0xcf9311ae
.word 0x3ff80427,0x543e1a12
.word 0xbc927c86,0x626d972b
.word 0x3ff814d2,0xadd106d9
.word 0x3c946437,0x0d151d4d
.word 0x3ff82589,0x994cce13
.word 0xbc9d4c1d,0xd41532d8
.word 0x3ff8364c,0x1eb941f7
.word 0x3c999b9a,0x31df2bd5
.word 0x3ff8471a,0x4623c7ad
.word 0xbc88d684,0xa341cdfb
.word 0x3ff857f4,0x179f5b21
.word 0xbc5ba748,0xf8b216d0
.word 0x3ff868d9,0x9b4492ec
.word 0x3ca01c83,0xb21584a3
.word 0x3ff879ca,0xd931a436
.word 0x3c85d2d7,0xd2db47bc
.word 0x3ff88ac7,0xd98a6699
.word 0x3c9994c2,0xf37cb53a
.word 0x3ff89bd0,0xa478580f
.word 0x3c9d5395,0x4475202a
.word 0x3ff8ace5,0x422aa0db
.word 0x3c96e9f1,0x56864b27
.word 0x3ff8be05,0xbad61778
.word 0x3c9ecb5e,0xfc43446e
.word 0x3ff8cf32,0x16b5448c
.word 0xbc70d55e,0x32e9e3aa
.word 0x3ff8e06a,0x5e0866d9
.word 0xbc97114a,0x6fc9b2e6
.word 0x3ff8f1ae,0x99157736
.word 0x3c85cc13,0xa2e3976c
.word 0x3ff902fe,0xd0282c8a
.word 0x3c9592ca,0x85fe3fd2
.word 0x3ff9145b,0x0b91ffc6
.word 0xbc9dd679,0x2e582524
.word 0x3ff925c3,0x53aa2fe2
.word 0xbc83455f,0xa639db7f
.word 0x3ff93737,0xb0cdc5e5
.word 0xbc675fc7,0x81b57ebc
.word 0x3ff948b8,0x2b5f98e5
.word 0xbc8dc3d6,0x797d2d99
.word 0x3ff95a44,0xcbc8520f
.word 0xbc764b7c,0x96a5f039
.word 0x3ff96bdd,0x9a7670b3
.word 0xbc5ba596,0x7f19c896
.word 0x3ff97d82,0x9fde4e50
.word 0xbc9d185b,0x7c1b85d0
.word 0x3ff98f33,0xe47a22a2
.word 0x3c7cabda,0xa24c78ed
.word 0x3ff9a0f1,0x70ca07ba
.word 0xbc9173bd,0x91cee632
.word 0x3ff9b2bb,0x4d53fe0d
.word 0xbc9dd84e,0x4df6d518
.word 0x3ff9c491,0x82a3f090
.word 0x3c7c7c46,0xb071f2be
.word 0x3ff9d674,0x194bb8d5
.word 0xbc9516be,0xa3dd8233
.word 0x3ff9e863,0x19e32323
.word 0x3c7824ca,0x78e64c6e
.word 0x3ff9fa5e,0x8d07f29e
.word 0xbc84a9ce,0xaaf1face
.word 0x3ffa0c66,0x7b5de565
.word 0xbc935949,0x5d1cd533
.word 0x3ffa1e7a,0xed8eb8bb
.word 0x3c9c6618,0xee8be70e
.word 0x3ffa309b,0xec4a2d33
.word 0x3c96305c,0x7ddc36ab
.word 0x3ffa42c9,0x80460ad8
.word 0xbc9aa780,0x589fb120
.word 0x3ffa5503,0xb23e255d
.word 0xbc9d2f6e,0xdb8d41e1
.word 0x3ffa674a,0x8af46052
.word 0x3c650f56,0x30670366
.word 0x3ffa799e,0x1330b358
.word 0x3c9bcb7e,0xcac563c6
.word 0x3ffa8bfe,0x53c12e59
.word 0xbc94f867,0xb2ba15a8
.word 0x3ffa9e6b,0x5579fdbf
.word 0x3c90fac9,0x0ef7fd31
.word 0x3ffab0e5,0x21356eba
.word 0x3c889c31,0xdae94544
.word 0x3ffac36b,0xbfd3f37a
.word 0xbc8f9234,0xcae76cd0
.word 0x3ffad5ff,0x3a3c2774
.word 0x3c97ef3b,0xb6b1b8e4
.word 0x3ffae89f,0x995ad3ad
.word 0x3c97a1cd,0x345dcc81
.word 0x3ffafb4c,0xe622f2ff
.word 0xbc94b2fc,0x0f315ecc
.word 0x3ffb0e07,0x298db666
.word 0xbc9bdef5,0x4c80e425
.word 0x3ffb20ce,0x6c9a8952
.word 0x3c94dd02,0x4a0756cc
.word 0x3ffb33a2,0xb84f15fb
.word 0xbc62805e,0x3084d708
.word 0x3ffb4684,0x15b749b1
.word 0xbc7f763d,0xe9df7c90
.word 0x3ffb5972,0x8de5593a
.word 0xbc9c71df,0xbbba6de3
.word 0x3ffb6c6e,0x29f1c52a
.word 0x3c92a8f3,0x52883f6e
.word 0x3ffb7f76,0xf2fb5e47
.word 0xbc75584f,0x7e54ac3b
.word 0x3ffb928c,0xf22749e4
.word 0xbc9b7216,0x54cb65c6
.word 0x3ffba5b0,0x30a1064a
.word 0xbc9efcd3,0x0e54292e
.word 0x3ffbb8e0,0xb79a6f1f
.word 0xbc3f52d1,0xc9696205
.word 0x3ffbcc1e,0x904bc1d2
.word 0x3c823dd0,0x7a2d9e84
.word 0x3ffbdf69,0xc3f3a207
.word 0xbc3c2623,0x60ea5b52
.word 0x3ffbf2c2,0x5bd71e09
.word 0xbc9efdca,0x3f6b9c73
.word 0x3ffc0628,0x6141b33d
.word 0xbc8d8a5a,0xa1fbca34
.word 0x3ffc199b,0xdd85529c
.word 0x3c811065,0x895048dd
.word 0x3ffc2d1c,0xd9fa652c
.word 0xbc96e516,0x17c8a5d7
.word 0x3ffc40ab,0x5fffd07a
.word 0x3c9b4537,0xe083c60a
.word 0x3ffc5447,0x78fafb22
.word 0x3c912f07,0x2493b5af
.word 0x3ffc67f1,0x2e57d14b
.word 0x3c92884d,0xff483cad
.word 0x3ffc7ba8,0x8988c933
.word 0xbc8e76bb,0xbe255559
.word 0x3ffc8f6d,0x9406e7b5
.word 0x3c71acbc,0x48805c44
.word 0x3ffca340,0x5751c4db
.word 0xbc87f2be,0xd10d08f4
.word 0x3ffcb720,0xdcef9069
.word 0x3c7503cb,0xd1e949db
.word 0x3ffccb0f,0x2e6d1675
.word 0xbc7d220f,0x86009093
.word 0x3ffcdf0b,0x555dc3fa
.word 0xbc8dd83b,0x53829d72
.word 0x3ffcf315,0x5b5bab74
.word 0xbc9a08e9,0xb86dff57
.word 0x3ffd072d,0x4a07897c
.word 0xbc9cbc37,0x43797a9c
.word 0x3ffd1b53,0x2b08c968
.word 0x3c955636,0x219a36ee
.word 0x3ffd2f87,0x080d89f2
.word 0xbc9d487b,0x719d8578
.word 0x3ffd43c8,0xeacaa1d6
.word 0x3c93db53,0xbf5a1614
.word 0x3ffd5818,0xdcfba487
.word 0x3c82ed02,0xd75b3706
.word 0x3ffd6c76,0xe862e6d3
.word 0x3c5fe87a,0x4a8165a0
.word 0x3ffd80e3,0x16c98398
.word 0xbc911ec1,0x8beddfe8
.word 0x3ffd955d,0x71ff6075
.word 0x3c9a052d,0xbb9af6be
.word 0x3ffda9e6,0x03db3285
.word 0x3c9c2300,0x696db532
.word 0x3ffdbe7c,0xd63a8315
.word 0xbc9b76f1,0x926b8be4
.word 0x3ffdd321,0xf301b460
.word 0x3c92da57,0x78f018c2
.word 0x3ffde7d5,0x641c0658
.word 0xbc9ca552,0x8e79ba8f
.word 0x3ffdfc97,0x337b9b5f
.word 0xbc91a5cd,0x4f184b5c
.word 0x3ffe1167,0x6b197d17
.word 0xbc72b529,0xbd5c7f44
.word 0x3ffe2646,0x14f5a129
.word 0xbc97b627,0x817a1496
.word 0x3ffe3b33,0x3b16ee12
.word 0xbc99f4a4,0x31fdc68a
.word 0x3ffe502e,0xe78b3ff6
.word 0x3c839e89,0x80a9cc8f
.word 0x3ffe6539,0x24676d76
.word 0xbc863ff8,0x7522b734
.word 0x3ffe7a51,0xfbc74c83
.word 0x3c92d522,0xca0c8de2
.word 0x3ffe8f79,0x77cdb740
.word 0xbc910894,0x80b054b1
.word 0x3ffea4af,0xa2a490da
.word 0xbc9e9c23,0x179c2893
.word 0x3ffeb9f4,0x867cca6e
.word 0x3c94832f,0x2293e4f2
.word 0x3ffecf48,0x2d8e67f1
.word 0xbc9c93f3,0xb411ad8c
.word 0x3ffee4aa,0xa2188510
.word 0x3c91c68d,0xa487568d
.word 0x3ffefa1b,0xee615a27
.word 0x3c9dc7f4,0x86a4b6b0
.word 0x3fff0f9c,0x1cb6412a
.word 0xbc932200,0x65181d45
.word 0x3fff252b,0x376bba97
.word 0x3c93a1a5,0xbf0d8e43
.word 0x3fff3ac9,0x48dd7274
.word 0xbc795a5a,0x3ed837de
.word 0x3fff5076,0x5b6e4540
.word 0x3c99d3e1,0x2dd8a18b
.word 0x3fff6632,0x798844f8
.word 0x3c9fa37b,0x3539343e
.word 0x3fff7bfd,0xad9cbe14
.word 0xbc9dbb12,0xd006350a
.word 0x3fff91d8,0x02243c89
.word 0xbc612ea8,0xa779f689
.word 0x3fffa7c1,0x819e90d8
.word 0x3c874853,0xf3a5931e
.word 0x3fffbdba,0x3692d514
.word 0xbc796773,0x15098eb6
.word 0x3fffd3c2,0x2b8f71f1
.word 0x3c62eb74,0x966579e7
.word 0x3fffe9d9,0x6b2a23d9
.word 0x3c74a603,0x7442fde3
.align 16
constants:
.word 0x3ef00000,0x00000000
.word 0x40862e42,0xfefa39ef
.word 0x01000000,0x00000000
.word 0x7f000000,0x00000000
.word 0x80000000,0x00000000
.word 0x43f00000,0x00000000 ! scaling 2^12 two96
.word 0xfff00000,0x00000000
.word 0x3ff00000,0x00000000
.word 0x3fdfffff,0xfffffff6
.word 0x3fc55555,0x721a1d14
.word 0x3fa55555,0x6e0896af
.word 0x41371547,0x652b82fe ! scaling 2^12 invln2_256
.word 0x3ea62e42,0xfee00000 ! scaling 2^(-12) ln2_256h
.word 0x3caa39ef,0x35793c76 ! scaling 2^(-12) ln2_256l
! base set w/o scaling
! .word 0x43300000,0x00000000 ! scaling two96
! .word 0x40771547,0x652b82fe ! scaling invln2_256
! .word 0x3f662e42,0xfee00000 ! scaling ln2_256h
! .word 0x3d6a39ef,0x35793c76 ! scaling ln2_256l
#define ox3ef 0x0
#define thresh 0x8
#define tiny 0x10
#define huge 0x18
#define signbit 0x20
#define two96 0x28
#define neginf 0x30
#define one 0x38
#define B1OFF 0x40
#define B2OFF 0x48
#define B3OFF 0x50
#define invln2_256 0x58
#define ln2_256h 0x60
#define ln2_256l 0x68
! local storage indices
#define m2 STACK_BIAS-0x4
#define m1 STACK_BIAS-0x8
#define m0 STACK_BIAS-0xc
#define jnk STACK_BIAS-0x20
! sizeof temp storage - must be a multiple of 16 for V9
#define tmps 0x20
! register use
! i0 n
! i1 x
! i2 stridex
! i3 y
! i4 stridey
! i5 0x80000000
! g1 TBL
! l0 m0
! l1 m1
! l2 m2
! l3 j0,oy0
! l4 j1,oy1
! l5 j2,oy2
! l6 0x3e300000
! l7 0x40862e41
! o0 py0
! o1 py1
! o2 py2
! o3 scratch
! o4 scratch
! o5 0x40874910
! o7 0x7ff00000
! f0 x0
! f2
! f4
! f6
! f8
! f10 x1
! f12
! f14
! f16
! f18
! f20 x2
! f22
! f24
! f26
! f28
! f30
! f32
! f34
! f36 0x3ef0...
! f38 thresh
! f40 tiny
! f42 huge
! f44 signbit
! f46 two96
! f48 neginf
! f50 one
! f52 B1
! f54 B2
! f56 B3
! f58 invln2_256
! f60 ln2_256h
! f62 ln2_256l
#define BOUNDRY %f36
#define THRESH %f38
#define TINY %f40
#define HUGE %f42
#define SIGNBIT %f44
#define TWO96 %f46
#define NEGINF %f48
#define ONE %f50
#define B1 %f52
#define B2 %f54
#define B3 %f56
#define INVLN2_256 %f58
#define LN2_256H %f60
#define LN2_256L %f62
ENTRY(__vexp)
save %sp,-SA(MINFRAME)-tmps,%sp
PIC_SETUP(l7)
PIC_SET(l7,constants,o3)
PIC_SET(l7,TBL,o0)
mov %o0,%g1
wr %g0,0x82,%asi ! set %asi for non-faulting loads
sethi %hi(0x80000000),%i5
sethi %hi(0x3e300000),%l6
sethi %hi(0x40862e41),%l7
or %l7,%lo(0x40862e41),%l7
sethi %hi(0x40874910),%o5
or %o5,%lo(0x40874910),%o5
sethi %hi(0x7ff00000),%o7
ldd [%o3+ox3ef],BOUNDRY
ldd [%o3+thresh],THRESH
ldd [%o3+tiny],TINY
ldd [%o3+huge],HUGE
ldd [%o3+signbit],SIGNBIT
ldd [%o3+two96],TWO96
ldd [%o3+neginf],NEGINF
ldd [%o3+one],ONE
ldd [%o3+B1OFF],B1
ldd [%o3+B2OFF],B2
ldd [%o3+B3OFF],B3
ldd [%o3+invln2_256],INVLN2_256
ldd [%o3+ln2_256h],LN2_256H
ldd [%o3+ln2_256l],LN2_256L
sll %i2,3,%i2 ! scale strides
sll %i4,3,%i4
add %fp,jnk,%l3 ! precondition loop
add %fp,jnk,%l4
add %fp,jnk,%l5
ld [%i1],%l0 ! hx = *x
ld [%i1],%f0
ld [%i1+4],%f1
andn %l0,%i5,%l0 ! hx &= ~0x80000000
ba .loop0
add %i1,%i2,%i1 ! x += stridex
.align 16
! -- 16 byte aligned
.loop0:
lda [%i1]%asi,%l1 ! preload next argument
sub %l0,%l6,%o3
sub %l7,%l0,%o4
fand %f0,SIGNBIT,%f2 ! get sign bit
lda [%i1]%asi,%f10
orcc %o3,%o4,%g0
mov %i3,%o0 ! py0 = y
bl,pn %icc,.range0 ! if hx < 0x3e300000 or > 0x40862e41
! delay slot
lda [%i1+4]%asi,%f11
addcc %i0,-1,%i0
add %i3,%i4,%i3 ! y += stridey
ble,pn %icc,.endloop1
! delay slot
andn %l1,%i5,%l1
add %i1,%i2,%i1 ! x += stridex
for %f2,TWO96,%f2 ! used to strip least sig bits
fmuld %f0,INVLN2_256,%f4 ! x/ (ln2/256) , creating k
.loop1:
lda [%i1]%asi,%l2 ! preload next argument
sub %l1,%l6,%o3
sub %l7,%l1,%o4
fand %f10,SIGNBIT,%f12
lda [%i1]%asi,%f20
orcc %o3,%o4,%g0
mov %i3,%o1 ! py1 = y
bl,pn %icc,.range1 ! if hx < 0x3e300000 or > 0x40862e41
! delay slot
lda [%i1+4]%asi,%f21
addcc %i0,-1,%i0
add %i3,%i4,%i3 ! y += stridey
ble,pn %icc,.endloop2
! delay slot
andn %l2,%i5,%l2
add %i1,%i2,%i1 ! x += stridex
for %f12,TWO96,%f12
fmuld %f10,INVLN2_256,%f14
.loop2:
sub %l2,%l6,%o3
sub %l7,%l2,%o4
fand %f20,SIGNBIT,%f22
fmuld %f20,INVLN2_256,%f24 ! okay to put this here; for alignment
orcc %o3,%o4,%g0
bl,pn %icc,.range2 ! if hx < 0x3e300000 or > 0x40862e41
! delay slot
for %f22,TWO96,%f22
faddd %f4,%f2,%f4 ! creating k+j/256, sra to zero bits
.cont:
faddd %f14,%f12,%f14
mov %i3,%o2 ! py2 = y
faddd %f24,%f22,%f24
add %i3,%i4,%i3 ! y += stridey
! BUBBLE USIII
fsubd %f4,%f2,%f8 ! creating k+j/256: sll
st %f6,[%l3] ! store previous loop x0
fsubd %f14,%f12,%f18
st %f7,[%l3+4] ! store previous loop x0
fsubd %f24,%f22,%f28
st %f16,[%l4]
! BUBBLE USIII
fmuld %f8,LN2_256H,%f2 ! closest LN2_256 to x
st %f17,[%l4+4]
fmuld %f18,LN2_256H,%f12
st %f26,[%l5]
fmuld %f28,LN2_256H,%f22
st %f27,[%l5+4]
! BUBBLE USIII
fsubd %f0,%f2,%f0 ! r = x - p*LN2_256H
fmuld %f8,LN2_256L,%f4 ! closest LN2_256 to x , added prec
fsubd %f10,%f12,%f10
fmuld %f18,LN2_256L,%f14
fsubd %f20,%f22,%f20
fmuld %f28,LN2_256L,%f24
! BUBBLE USIII
fsubd %f0,%f4,%f0 ! r -= p*LN2_256L
fsubd %f10,%f14,%f10
fsubd %f20,%f24,%f20
!!!!!!!!!!!!!!!!!!! New polynomial reorder starts here
! Alternate polynomial grouping allowing non-sequential calc of p
! OLD : p = r * ( 1 + r * ( B1 + r * ( B2 + r * B3) ) )
! NEW : p = r * [ (1+r*B1) + (r*r) * ( B2 + r * B3) ) ]
!
! let SLi Ri SRi be accumulators
fmuld %f0,B3,%f2 ! SR1 = r1 * B3
fdtoi %f8,%f8 ! convert k+j/256 to int
st %f8,[%fp+m0] ! store k, to shift return/use
fmuld %f10,B3,%f12 ! SR2 = r2 * B3
fdtoi %f18,%f18 ! convert k+j/256 to int
st %f18,[%fp+m1] ! store k, to shift return/use
fmuld %f20,B3,%f22 ! SR3 = r3 * B3
fdtoi %f28,%f28 ! convert k+j/256 to int
st %f28,[%fp+m2] ! store k, to shift return/use
fmuld %f0,%f0,%f4 ! R1 = r1 * r1
fmuld %f10,%f10,%f14 ! R2 = r2 * r2
faddd %f2,B2,%f2 ! SR1 += B2
fmuld %f20,%f20,%f24 ! R3 = r3 * r3
faddd %f12,B2,%f12 ! SR2 += B2
faddd %f22,B2,%f22 ! SR3 += B2
fmuld %f0,B1,%f6 ! SL1 = r1 * B1
fmuld %f10,B1,%f32 ! SL2 = r2 * B1
fand %f8,NEGINF,%f8
! best here for RAW BYPASS
ld [%fp+m0],%l0 ! get nonshifted k into intreg
fmuld %f20,B1,%f34 ! SL3 = r3 * B1
fand %f18,NEGINF,%f18
ld [%fp+m1],%l1 ! get nonshifted k into intreg
fmuld %f4,%f2,%f4 ! R1 = R1 * SR1
fand %f28,NEGINF,%f28
ld [%fp+m2],%l2 ! get nonshifted k into intreg
fmuld %f14,%f12,%f14 ! R2 = R2 * SR2
faddd %f6,ONE,%f6 ! SL1 += 1
fmuld %f24,%f22,%f24 ! R3 = R3 * SR3
faddd %f32,ONE,%f32 ! SL2 += 1
sra %l0,8,%l3 ! shift k tobe offset 256-8byte
faddd %f34,ONE,%f34 ! SL3 += 1
sra %l1,8,%l4 ! shift k tobe offset 256-8byte
sra %l2,8,%l5 ! shift k tobe offset 256-8byte
! BUBBLE in USIII
and %l3,0xff0,%l3
and %l4,0xff0,%l4
faddd %f6,%f4,%f6 ! R1 = SL1 + R1
ldd [%g1+%l3],%f4 ! tbl[j]
add %l3,8,%l3 ! inc j
and %l5,0xff0,%l5
faddd %f32,%f14,%f32 ! R2 = SL2 + R2
ldd [%g1+%l4],%f14 ! tbl[j]
add %l4,8,%l4 ! inc j
sra %l0,20,%o3
faddd %f34,%f24,%f34 ! R3 = SL3 + R3
ldd [%g1+%l5],%f24 ! tbl[j]
add %l5,8,%l5 ! inc j
sra %l1,20,%l1
! BUBBLE in USIII
ldd [%g1+%l4],%f16 ! tbl[j+1]
add %o3,1021,%o3 ! inc j
fmuld %f0,%f6,%f0 ! p1 = r1 * R1
ldd [%g1+%l3],%f6 ! tbl[j+1]
add %l1,1021,%l1 ! inc j
sra %l2,20,%l2
fmuld %f10,%f32,%f10 ! p2 = r2 * R2
ldd [%g1+%l5],%f26 ! tbl[j+1]
add %l2,1021,%l2 ! inc j
fmuld %f20,%f34,%f20 ! p3 = r3 * R3
!!!!!!!!!!!!!!!!!!! poly-reorder - ends here
fmuld %f0,%f4,%f0 ! start exp(x) = exp(r) * tbl[j]
mov %o0,%l3
fmuld %f10,%f14,%f10
mov %o1,%l4
fmuld %f20,%f24,%f20
mov %o2,%l5
faddd %f0,%f6,%f6 ! cont exp(x) : apply tbl[j] high bits
lda [%i1]%asi,%l0 ! preload next argument
faddd %f10,%f16,%f16
lda [%i1]%asi,%f0
faddd %f20,%f26,%f26
lda [%i1+4]%asi,%f1
faddd %f6,%f4,%f6 ! cont exp(x) : apply tbl[j+1] low bits
add %i1,%i2,%i1 ! x += stridex
faddd %f16,%f14,%f16
andn %l0,%i5,%l0
or %o3,%l1,%o4
! -- 16 byte aligned
orcc %o4,%l2,%o4
bl,pn %icc,.small
! delay slot
faddd %f26,%f24,%f26
fpadd32 %f6,%f8,%f6 ! done exp(x) : apply 2^k
fpadd32 %f16,%f18,%f16
addcc %i0,-1,%i0
bg,pn %icc,.loop0
! delay slot
fpadd32 %f26,%f28,%f26
ba,pt %icc,.endloop0
! delay slot
nop
.align 16
.small:
tst %o3
bge,pt %icc,1f
! delay slot
fpadd32 %f6,%f8,%f6
fpadd32 %f6,BOUNDRY,%f6
fmuld %f6,TINY,%f6
1:
tst %l1
bge,pt %icc,1f
! delay slot
fpadd32 %f16,%f18,%f16
fpadd32 %f16,BOUNDRY,%f16
fmuld %f16,TINY,%f16
1:
tst %l2
bge,pt %icc,1f
! delay slot
fpadd32 %f26,%f28,%f26
fpadd32 %f26,BOUNDRY,%f26
fmuld %f26,TINY,%f26
1:
addcc %i0,-1,%i0
bg,pn %icc,.loop0
! delay slot
nop
ba,pt %icc,.endloop0
! delay slot
nop
.endloop2:
for %f12,TWO96,%f12
fmuld %f10,INVLN2_256,%f14
faddd %f14,%f12,%f14
fsubd %f14,%f12,%f18
fmuld %f18,LN2_256H,%f12
fsubd %f10,%f12,%f10
fmuld %f18,LN2_256L,%f14
fsubd %f10,%f14,%f10
fmuld %f10,B3,%f12
fdtoi %f18,%f18
st %f18,[%fp+m1]
fmuld %f10,%f10,%f14
faddd %f12,B2,%f12
fmuld %f10,B1,%f32
fand %f18,NEGINF,%f18
ld [%fp+m1],%l1
fmuld %f14,%f12,%f14
faddd %f32,ONE,%f32
sra %l1,8,%o4
and %o4,0xff0,%o4
faddd %f32,%f14,%f32
ldd [%g1+%o4],%f14
add %o4,8,%o4
sra %l1,20,%l1
ldd [%g1+%o4],%f30
addcc %l1,1021,%l1
fmuld %f10,%f32,%f10
fmuld %f10,%f14,%f10
faddd %f10,%f30,%f30
faddd %f30,%f14,%f30
bge,pt %icc,1f
! delay slot
fpadd32 %f30,%f18,%f30
fpadd32 %f30,BOUNDRY,%f30
fmuld %f30,TINY,%f30
1:
st %f30,[%o1]
st %f31,[%o1+4]
.endloop1:
for %f2,TWO96,%f2
fmuld %f0,INVLN2_256,%f4
faddd %f4,%f2,%f4
fsubd %f4,%f2,%f8
fmuld %f8,LN2_256H,%f2
fsubd %f0,%f2,%f0
fmuld %f8,LN2_256L,%f4
fsubd %f0,%f4,%f0
fmuld %f0,B3,%f2
fdtoi %f8,%f8
st %f8,[%fp+m0]
fmuld %f0,%f0,%f4
faddd %f2,B2,%f2
fmuld %f0,B1,%f32
fand %f8,NEGINF,%f8
ld [%fp+m0],%l0
fmuld %f4,%f2,%f4
faddd %f32,ONE,%f32
sra %l0,8,%o4
and %o4,0xff0,%o4
faddd %f32,%f4,%f32
ldd [%g1+%o4],%f4
add %o4,8,%o4
sra %l0,20,%o3
ldd [%g1+%o4],%f30
addcc %o3,1021,%o3
fmuld %f0,%f32,%f0
fmuld %f0,%f4,%f0
faddd %f0,%f30,%f30
faddd %f30,%f4,%f30
bge,pt %icc,1f
! delay slot
fpadd32 %f30,%f8,%f30
fpadd32 %f30,BOUNDRY,%f30
fmuld %f30,TINY,%f30
1:
st %f30,[%o0]
st %f31,[%o0+4]
.endloop0:
st %f6,[%l3]
st %f7,[%l3+4]
st %f16,[%l4]
st %f17,[%l4+4]
st %f26,[%l5]
st %f27,[%l5+4]
ret
restore
.range0:
cmp %l0,%l6
bl,a,pt %icc,3f ! if x is tiny
! delay slot, annulled if branch not taken
faddd %f0,ONE,%f4
cmp %l0,%o5
bg,pt %icc,1f ! if x is huge, inf, nan
! delay slot
nop
fcmpd %fcc0,%f0,THRESH
fbg,a,pt %fcc0,3f ! if x is huge and positive
! delay slot, annulled if branch not taken
fmuld HUGE,HUGE,%f4
! x is near the extremes but within range; return to the loop
addcc %i0,-1,%i0
add %i3,%i4,%i3 ! y += stridey
ble,pn %icc,.endloop1
! delay slot
andn %l1,%i5,%l1
add %i1,%i2,%i1 ! x += stridex
for %f2,TWO96,%f2
ba,pt %icc,.loop1
! delay slot
fmuld %f0,INVLN2_256,%f4
1:
cmp %l0,%o7
bl,pn %icc,2f ! if x is finite
! delay slot
nop
fzero %f4
fcmpd %fcc0,%f0,NEGINF
fmovdne %fcc0,%f0,%f4
ba,pt %icc,3f
fmuld %f4,%f4,%f4 ! x*x or zero*zero
2:
fmovd HUGE,%f4
fcmpd %fcc0,%f0,ONE
fmovdl %fcc0,TINY,%f4
fmuld %f4,%f4,%f4 ! huge*huge or tiny*tiny
3:
st %f4,[%o0]
andn %l1,%i5,%l0
add %i1,%i2,%i1 ! x += stridex
fmovd %f10,%f0
st %f5,[%o0+4]
addcc %i0,-1,%i0
bg,pt %icc,.loop0
! delay slot
add %i3,%i4,%i3 ! y += stridey
ba,pt %icc,.endloop0
! delay slot
nop
.range1:
cmp %l1,%l6
bl,a,pt %icc,3f ! if x is tiny
! delay slot, annulled if branch not taken
faddd %f10,ONE,%f14
cmp %l1,%o5
bg,pt %icc,1f ! if x is huge, inf, nan
! delay slot
nop
fcmpd %fcc0,%f10,THRESH
fbg,a,pt %fcc0,3f ! if x is huge and positive
! delay slot, annulled if branch not taken
fmuld HUGE,HUGE,%f14
! x is near the extremes but within range; return to the loop
addcc %i0,-1,%i0
add %i3,%i4,%i3 ! y += stridey
ble,pn %icc,.endloop2
! delay slot
andn %l2,%i5,%l2
add %i1,%i2,%i1 ! x += stridex
for %f12,TWO96,%f12
ba,pt %icc,.loop2
! delay slot
fmuld %f10,INVLN2_256,%f14
1:
cmp %l1,%o7
bl,pn %icc,2f ! if x is finite
! delay slot
nop
fzero %f14
fcmpd %fcc0,%f10,NEGINF
fmovdne %fcc0,%f10,%f14
ba,pt %icc,3f
fmuld %f14,%f14,%f14 ! x*x or zero*zero
2:
fmovd HUGE,%f14
fcmpd %fcc0,%f10,ONE
fmovdl %fcc0,TINY,%f14
fmuld %f14,%f14,%f14 ! huge*huge or tiny*tiny
3:
st %f14,[%o1]
andn %l2,%i5,%l1
add %i1,%i2,%i1 ! x += stridex
fmovd %f20,%f10
st %f15,[%o1+4]
addcc %i0,-1,%i0
bg,pt %icc,.loop1
! delay slot
add %i3,%i4,%i3 ! y += stridey
ba,pt %icc,.endloop1
! delay slot
nop
.range2:
cmp %l2,%l6
bl,a,pt %icc,3f ! if x is tiny
! delay slot, annulled if branch not taken
faddd %f20,ONE,%f24
cmp %l2,%o5
bg,pt %icc,1f ! if x is huge, inf, nan
! delay slot
nop
fcmpd %fcc0,%f20,THRESH
fbg,a,pt %fcc0,3f ! if x is huge and positive
! delay slot, annulled if branch not taken
fmuld HUGE,HUGE,%f24
! x is near the extremes but within range; return to the loop
ba,pt %icc,.cont
! delay slot
faddd %f4,%f2,%f4
1:
cmp %l2,%o7
bl,pn %icc,2f ! if x is finite
! delay slot
nop
fzero %f24
fcmpd %fcc0,%f20,NEGINF
fmovdne %fcc0,%f20,%f24
ba,pt %icc,3f
fmuld %f24,%f24,%f24 ! x*x or zero*zero
2:
fmovd HUGE,%f24
fcmpd %fcc0,%f20,ONE
fmovdl %fcc0,TINY,%f24
fmuld %f24,%f24,%f24 ! huge*huge or tiny*tiny
3:
st %f24,[%i3]
st %f25,[%i3+4]
lda [%i1]%asi,%l2 ! preload next argument
lda [%i1]%asi,%f20
lda [%i1+4]%asi,%f21
andn %l2,%i5,%l2
add %i1,%i2,%i1 ! x += stridex
addcc %i0,-1,%i0
bg,pt %icc,.loop2
! delay slot
add %i3,%i4,%i3 ! y += stridey
ba,pt %icc,.endloop2
! delay slot
nop
SET_SIZE(__vexp)
