mpsignal.c 22 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967
  1. /*
  2. * Copyright (c) 2008-2020 Stefan Krah. All rights reserved.
  3. *
  4. * Redistribution and use in source and binary forms, with or without
  5. * modification, are permitted provided that the following conditions
  6. * are met:
  7. *
  8. * 1. Redistributions of source code must retain the above copyright
  9. * notice, this list of conditions and the following disclaimer.
  10. *
  11. * 2. Redistributions in binary form must reproduce the above copyright
  12. * notice, this list of conditions and the following disclaimer in the
  13. * documentation and/or other materials provided with the distribution.
  14. *
  15. * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS" AND
  16. * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  17. * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  18. * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
  19. * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  20. * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  21. * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  22. * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  23. * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  24. * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  25. * SUCH DAMAGE.
  26. */
  27. #include "mpdecimal.h"
  28. #include <stddef.h>
  29. #include <stdint.h>
  30. /* Signaling wrappers for the quiet functions in mpdecimal.c. */
  31. char *
  32. mpd_format(const mpd_t *dec, const char *fmt, mpd_context_t *ctx)
  33. {
  34. char *ret;
  35. uint32_t status = 0;
  36. ret = mpd_qformat(dec, fmt, ctx, &status);
  37. mpd_addstatus_raise(ctx, status);
  38. return ret;
  39. }
  40. void
  41. mpd_import_u16(mpd_t *result, const uint16_t *srcdata, size_t srclen,
  42. uint8_t srcsign, uint32_t base, mpd_context_t *ctx)
  43. {
  44. uint32_t status = 0;
  45. mpd_qimport_u16(result, srcdata, srclen, srcsign, base, ctx, &status);
  46. mpd_addstatus_raise(ctx, status);
  47. }
  48. void
  49. mpd_import_u32(mpd_t *result, const uint32_t *srcdata, size_t srclen,
  50. uint8_t srcsign, uint32_t base, mpd_context_t *ctx)
  51. {
  52. uint32_t status = 0;
  53. mpd_qimport_u32(result, srcdata, srclen, srcsign, base, ctx, &status);
  54. mpd_addstatus_raise(ctx, status);
  55. }
  56. size_t
  57. mpd_export_u16(uint16_t **rdata, size_t rlen, uint32_t base, const mpd_t *src,
  58. mpd_context_t *ctx)
  59. {
  60. size_t n;
  61. uint32_t status = 0;
  62. n = mpd_qexport_u16(rdata, rlen, base, src, &status);
  63. mpd_addstatus_raise(ctx, status);
  64. return n;
  65. }
  66. size_t
  67. mpd_export_u32(uint32_t **rdata, size_t rlen, uint32_t base, const mpd_t *src,
  68. mpd_context_t *ctx)
  69. {
  70. size_t n;
  71. uint32_t status = 0;
  72. n = mpd_qexport_u32(rdata, rlen, base, src, &status);
  73. mpd_addstatus_raise(ctx, status);
  74. return n;
  75. }
  76. void
  77. mpd_finalize(mpd_t *result, mpd_context_t *ctx)
  78. {
  79. uint32_t status = 0;
  80. mpd_qfinalize(result, ctx, &status);
  81. mpd_addstatus_raise(ctx, status);
  82. }
  83. int
  84. mpd_check_nan(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  85. {
  86. uint32_t status = 0;
  87. if (mpd_qcheck_nan(result, a, ctx, &status)) {
  88. mpd_addstatus_raise(ctx, status);
  89. return 1;
  90. }
  91. return 0;
  92. }
  93. int
  94. mpd_check_nans(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  95. {
  96. uint32_t status = 0;
  97. if (mpd_qcheck_nans(result, a, b, ctx, &status)) {
  98. mpd_addstatus_raise(ctx, status);
  99. return 1;
  100. }
  101. return 0;
  102. }
  103. void
  104. mpd_set_string(mpd_t *result, const char *s, mpd_context_t *ctx)
  105. {
  106. uint32_t status = 0;
  107. mpd_qset_string(result, s, ctx, &status);
  108. mpd_addstatus_raise(ctx, status);
  109. }
  110. void
  111. mpd_maxcoeff(mpd_t *result, mpd_context_t *ctx)
  112. {
  113. uint32_t status = 0;
  114. mpd_qmaxcoeff(result, ctx, &status);
  115. mpd_addstatus_raise(ctx, status);
  116. }
  117. /* set static mpd from signed integer */
  118. void
  119. mpd_sset_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx)
  120. {
  121. uint32_t status = 0;
  122. mpd_qsset_ssize(result, a, ctx, &status);
  123. mpd_addstatus_raise(ctx, status);
  124. }
  125. void
  126. mpd_sset_i32(mpd_t *result, int32_t a, mpd_context_t *ctx)
  127. {
  128. uint32_t status = 0;
  129. mpd_qsset_i32(result, a, ctx, &status);
  130. mpd_addstatus_raise(ctx, status);
  131. }
  132. #ifdef CONFIG_64
  133. void
  134. mpd_sset_i64(mpd_t *result, int64_t a, mpd_context_t *ctx)
  135. {
  136. uint32_t status = 0;
  137. mpd_qsset_i64(result, a, ctx, &status);
  138. mpd_addstatus_raise(ctx, status);
  139. }
  140. #endif
  141. /* set static mpd from unsigned integer */
  142. void
  143. mpd_sset_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx)
  144. {
  145. uint32_t status = 0;
  146. mpd_qsset_uint(result, a, ctx, &status);
  147. mpd_addstatus_raise(ctx, status);
  148. }
  149. void
  150. mpd_sset_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx)
  151. {
  152. uint32_t status = 0;
  153. mpd_qsset_u32(result, a, ctx, &status);
  154. mpd_addstatus_raise(ctx, status);
  155. }
  156. #ifdef CONFIG_64
  157. void
  158. mpd_sset_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx)
  159. {
  160. uint32_t status = 0;
  161. mpd_qsset_u64(result, a, ctx, &status);
  162. mpd_addstatus_raise(ctx, status);
  163. }
  164. #endif
  165. /* set mpd from signed integer */
  166. void
  167. mpd_set_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx)
  168. {
  169. uint32_t status = 0;
  170. mpd_qset_ssize(result, a, ctx, &status);
  171. mpd_addstatus_raise(ctx, status);
  172. }
  173. void
  174. mpd_set_i32(mpd_t *result, int32_t a, mpd_context_t *ctx)
  175. {
  176. uint32_t status = 0;
  177. mpd_qset_i32(result, a, ctx, &status);
  178. mpd_addstatus_raise(ctx, status);
  179. }
  180. #ifndef LEGACY_COMPILER
  181. void
  182. mpd_set_i64(mpd_t *result, int64_t a, mpd_context_t *ctx)
  183. {
  184. uint32_t status = 0;
  185. mpd_qset_i64(result, a, ctx, &status);
  186. mpd_addstatus_raise(ctx, status);
  187. }
  188. #endif
  189. /* set mpd from unsigned integer */
  190. void
  191. mpd_set_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx)
  192. {
  193. uint32_t status = 0;
  194. mpd_qset_uint(result, a, ctx, &status);
  195. mpd_addstatus_raise(ctx, status);
  196. }
  197. void
  198. mpd_set_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx)
  199. {
  200. uint32_t status = 0;
  201. mpd_qset_u32(result, a, ctx, &status);
  202. mpd_addstatus_raise(ctx, status);
  203. }
  204. #ifndef LEGACY_COMPILER
  205. void
  206. mpd_set_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx)
  207. {
  208. uint32_t status = 0;
  209. mpd_qset_u64(result, a, ctx, &status);
  210. mpd_addstatus_raise(ctx, status);
  211. }
  212. #endif
  213. /* convert mpd to signed integer */
  214. mpd_ssize_t
  215. mpd_get_ssize(const mpd_t *a, mpd_context_t *ctx)
  216. {
  217. uint32_t status = 0;
  218. mpd_ssize_t ret;
  219. ret = mpd_qget_ssize(a, &status);
  220. mpd_addstatus_raise(ctx, status);
  221. return ret;
  222. }
  223. int32_t
  224. mpd_get_i32(const mpd_t *a, mpd_context_t *ctx)
  225. {
  226. uint32_t status = 0;
  227. int32_t ret;
  228. ret = mpd_qget_i32(a, &status);
  229. mpd_addstatus_raise(ctx, status);
  230. return ret;
  231. }
  232. #ifndef LEGACY_COMPILER
  233. int64_t
  234. mpd_get_i64(const mpd_t *a, mpd_context_t *ctx)
  235. {
  236. uint32_t status = 0;
  237. int64_t ret;
  238. ret = mpd_qget_i64(a, &status);
  239. mpd_addstatus_raise(ctx, status);
  240. return ret;
  241. }
  242. #endif
  243. mpd_uint_t
  244. mpd_get_uint(const mpd_t *a, mpd_context_t *ctx)
  245. {
  246. uint32_t status = 0;
  247. mpd_uint_t ret;
  248. ret = mpd_qget_uint(a, &status);
  249. mpd_addstatus_raise(ctx, status);
  250. return ret;
  251. }
  252. mpd_uint_t
  253. mpd_abs_uint(const mpd_t *a, mpd_context_t *ctx)
  254. {
  255. uint32_t status = 0;
  256. mpd_uint_t ret;
  257. ret = mpd_qabs_uint(a, &status);
  258. mpd_addstatus_raise(ctx, status);
  259. return ret;
  260. }
  261. uint32_t
  262. mpd_get_u32(const mpd_t *a, mpd_context_t *ctx)
  263. {
  264. uint32_t status = 0;
  265. uint32_t ret;
  266. ret = mpd_qget_u32(a, &status);
  267. mpd_addstatus_raise(ctx, status);
  268. return ret;
  269. }
  270. #ifndef LEGACY_COMPILER
  271. uint64_t
  272. mpd_get_u64(const mpd_t *a, mpd_context_t *ctx)
  273. {
  274. uint32_t status = 0;
  275. uint64_t ret;
  276. ret = mpd_qget_u64(a, &status);
  277. mpd_addstatus_raise(ctx, status);
  278. return ret;
  279. }
  280. #endif
  281. void
  282. mpd_and(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  283. {
  284. uint32_t status = 0;
  285. mpd_qand(result, a, b, ctx, &status);
  286. mpd_addstatus_raise(ctx, status);
  287. }
  288. void
  289. mpd_copy(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  290. {
  291. uint32_t status = 0;
  292. if (!mpd_qcopy(result, a, &status)) {
  293. mpd_addstatus_raise(ctx, status);
  294. }
  295. }
  296. void
  297. mpd_canonical(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  298. {
  299. mpd_copy(result, a, ctx);
  300. }
  301. void
  302. mpd_copy_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  303. {
  304. uint32_t status = 0;
  305. if (!mpd_qcopy_abs(result, a, &status)) {
  306. mpd_addstatus_raise(ctx, status);
  307. }
  308. }
  309. void
  310. mpd_copy_negate(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  311. {
  312. uint32_t status = 0;
  313. if (!mpd_qcopy_negate(result, a, &status)) {
  314. mpd_addstatus_raise(ctx, status);
  315. }
  316. }
  317. void
  318. mpd_copy_sign(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  319. {
  320. uint32_t status = 0;
  321. if (!mpd_qcopy_sign(result, a, b, &status)) {
  322. mpd_addstatus_raise(ctx, status);
  323. }
  324. }
  325. void
  326. mpd_invert(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  327. {
  328. uint32_t status = 0;
  329. mpd_qinvert(result, a, ctx, &status);
  330. mpd_addstatus_raise(ctx, status);
  331. }
  332. void
  333. mpd_logb(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  334. {
  335. uint32_t status = 0;
  336. mpd_qlogb(result, a, ctx, &status);
  337. mpd_addstatus_raise(ctx, status);
  338. }
  339. void
  340. mpd_or(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  341. {
  342. uint32_t status = 0;
  343. mpd_qor(result, a, b, ctx, &status);
  344. mpd_addstatus_raise(ctx, status);
  345. }
  346. void
  347. mpd_rotate(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  348. {
  349. uint32_t status = 0;
  350. mpd_qrotate(result, a, b, ctx, &status);
  351. mpd_addstatus_raise(ctx, status);
  352. }
  353. void
  354. mpd_scaleb(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  355. {
  356. uint32_t status = 0;
  357. mpd_qscaleb(result, a, b, ctx, &status);
  358. mpd_addstatus_raise(ctx, status);
  359. }
  360. void
  361. mpd_shiftl(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx)
  362. {
  363. uint32_t status = 0;
  364. mpd_qshiftl(result, a, n, &status);
  365. mpd_addstatus_raise(ctx, status);
  366. }
  367. mpd_uint_t
  368. mpd_shiftr(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx)
  369. {
  370. uint32_t status = 0;
  371. mpd_uint_t rnd;
  372. rnd = mpd_qshiftr(result, a, n, &status);
  373. mpd_addstatus_raise(ctx, status);
  374. return rnd;
  375. }
  376. void
  377. mpd_shiftn(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx)
  378. {
  379. uint32_t status = 0;
  380. mpd_qshiftn(result, a, n, ctx, &status);
  381. mpd_addstatus_raise(ctx, status);
  382. }
  383. void
  384. mpd_shift(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  385. {
  386. uint32_t status = 0;
  387. mpd_qshift(result, a, b, ctx, &status);
  388. mpd_addstatus_raise(ctx, status);
  389. }
  390. void
  391. mpd_xor(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  392. {
  393. uint32_t status = 0;
  394. mpd_qxor(result, a, b, ctx, &status);
  395. mpd_addstatus_raise(ctx, status);
  396. }
  397. void
  398. mpd_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  399. {
  400. uint32_t status = 0;
  401. mpd_qabs(result, a, ctx, &status);
  402. mpd_addstatus_raise(ctx, status);
  403. }
  404. int
  405. mpd_cmp(const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  406. {
  407. uint32_t status = 0;
  408. int c;
  409. c = mpd_qcmp(a, b, &status);
  410. mpd_addstatus_raise(ctx, status);
  411. return c;
  412. }
  413. int
  414. mpd_compare(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  415. {
  416. uint32_t status = 0;
  417. int c;
  418. c = mpd_qcompare(result, a, b, ctx, &status);
  419. mpd_addstatus_raise(ctx, status);
  420. return c;
  421. }
  422. int
  423. mpd_compare_signal(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  424. {
  425. uint32_t status = 0;
  426. int c;
  427. c = mpd_qcompare_signal(result, a, b, ctx, &status);
  428. mpd_addstatus_raise(ctx, status);
  429. return c;
  430. }
  431. void
  432. mpd_add(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  433. {
  434. uint32_t status = 0;
  435. mpd_qadd(result, a, b, ctx, &status);
  436. mpd_addstatus_raise(ctx, status);
  437. }
  438. void
  439. mpd_sub(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  440. {
  441. uint32_t status = 0;
  442. mpd_qsub(result, a, b, ctx, &status);
  443. mpd_addstatus_raise(ctx, status);
  444. }
  445. void
  446. mpd_add_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx)
  447. {
  448. uint32_t status = 0;
  449. mpd_qadd_ssize(result, a, b, ctx, &status);
  450. mpd_addstatus_raise(ctx, status);
  451. }
  452. void
  453. mpd_add_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx)
  454. {
  455. uint32_t status = 0;
  456. mpd_qadd_i32(result, a, b, ctx, &status);
  457. mpd_addstatus_raise(ctx, status);
  458. }
  459. #ifndef LEGACY_COMPILER
  460. void
  461. mpd_add_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx)
  462. {
  463. uint32_t status = 0;
  464. mpd_qadd_i64(result, a, b, ctx, &status);
  465. mpd_addstatus_raise(ctx, status);
  466. }
  467. #endif
  468. void
  469. mpd_add_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx)
  470. {
  471. uint32_t status = 0;
  472. mpd_qadd_uint(result, a, b, ctx, &status);
  473. mpd_addstatus_raise(ctx, status);
  474. }
  475. void
  476. mpd_add_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx)
  477. {
  478. uint32_t status = 0;
  479. mpd_qadd_u32(result, a, b, ctx, &status);
  480. mpd_addstatus_raise(ctx, status);
  481. }
  482. #ifndef LEGACY_COMPILER
  483. void
  484. mpd_add_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx)
  485. {
  486. uint32_t status = 0;
  487. mpd_qadd_u64(result, a, b, ctx, &status);
  488. mpd_addstatus_raise(ctx, status);
  489. }
  490. #endif
  491. void
  492. mpd_sub_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx)
  493. {
  494. uint32_t status = 0;
  495. mpd_qsub_ssize(result, a, b, ctx, &status);
  496. mpd_addstatus_raise(ctx, status);
  497. }
  498. void
  499. mpd_sub_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx)
  500. {
  501. uint32_t status = 0;
  502. mpd_qsub_i32(result, a, b, ctx, &status);
  503. mpd_addstatus_raise(ctx, status);
  504. }
  505. #ifndef LEGACY_COMPILER
  506. void
  507. mpd_sub_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx)
  508. {
  509. uint32_t status = 0;
  510. mpd_qsub_i64(result, a, b, ctx, &status);
  511. mpd_addstatus_raise(ctx, status);
  512. }
  513. #endif
  514. void
  515. mpd_sub_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx)
  516. {
  517. uint32_t status = 0;
  518. mpd_qsub_uint(result, a, b, ctx, &status);
  519. mpd_addstatus_raise(ctx, status);
  520. }
  521. void
  522. mpd_sub_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx)
  523. {
  524. uint32_t status = 0;
  525. mpd_qsub_u32(result, a, b, ctx, &status);
  526. mpd_addstatus_raise(ctx, status);
  527. }
  528. #ifndef LEGACY_COMPILER
  529. void
  530. mpd_sub_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx)
  531. {
  532. uint32_t status = 0;
  533. mpd_qsub_u64(result, a, b, ctx, &status);
  534. mpd_addstatus_raise(ctx, status);
  535. }
  536. #endif
  537. void
  538. mpd_div(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  539. {
  540. uint32_t status = 0;
  541. mpd_qdiv(q, a, b, ctx, &status);
  542. mpd_addstatus_raise(ctx, status);
  543. }
  544. void
  545. mpd_div_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx)
  546. {
  547. uint32_t status = 0;
  548. mpd_qdiv_ssize(result, a, b, ctx, &status);
  549. mpd_addstatus_raise(ctx, status);
  550. }
  551. void
  552. mpd_div_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx)
  553. {
  554. uint32_t status = 0;
  555. mpd_qdiv_i32(result, a, b, ctx, &status);
  556. mpd_addstatus_raise(ctx, status);
  557. }
  558. #ifndef LEGACY_COMPILER
  559. void
  560. mpd_div_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx)
  561. {
  562. uint32_t status = 0;
  563. mpd_qdiv_i64(result, a, b, ctx, &status);
  564. mpd_addstatus_raise(ctx, status);
  565. }
  566. #endif
  567. void
  568. mpd_div_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx)
  569. {
  570. uint32_t status = 0;
  571. mpd_qdiv_uint(result, a, b, ctx, &status);
  572. mpd_addstatus_raise(ctx, status);
  573. }
  574. void
  575. mpd_div_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx)
  576. {
  577. uint32_t status = 0;
  578. mpd_qdiv_u32(result, a, b, ctx, &status);
  579. mpd_addstatus_raise(ctx, status);
  580. }
  581. #ifndef LEGACY_COMPILER
  582. void
  583. mpd_div_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx)
  584. {
  585. uint32_t status = 0;
  586. mpd_qdiv_u64(result, a, b, ctx, &status);
  587. mpd_addstatus_raise(ctx, status);
  588. }
  589. #endif
  590. void
  591. mpd_divmod(mpd_t *q, mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  592. {
  593. uint32_t status = 0;
  594. mpd_qdivmod(q, r, a, b, ctx, &status);
  595. mpd_addstatus_raise(ctx, status);
  596. }
  597. void
  598. mpd_divint(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  599. {
  600. uint32_t status = 0;
  601. mpd_qdivint(q, a, b, ctx, &status);
  602. mpd_addstatus_raise(ctx, status);
  603. }
  604. void
  605. mpd_exp(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  606. {
  607. uint32_t status = 0;
  608. mpd_qexp(result, a, ctx, &status);
  609. mpd_addstatus_raise(ctx, status);
  610. }
  611. void
  612. mpd_fma(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_t *c,
  613. mpd_context_t *ctx)
  614. {
  615. uint32_t status = 0;
  616. mpd_qfma(result, a, b, c, ctx, &status);
  617. mpd_addstatus_raise(ctx, status);
  618. }
  619. void
  620. mpd_ln(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  621. {
  622. uint32_t status = 0;
  623. mpd_qln(result, a, ctx, &status);
  624. mpd_addstatus_raise(ctx, status);
  625. }
  626. void
  627. mpd_log10(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  628. {
  629. uint32_t status = 0;
  630. mpd_qlog10(result, a, ctx, &status);
  631. mpd_addstatus_raise(ctx, status);
  632. }
  633. void
  634. mpd_max(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  635. {
  636. uint32_t status = 0;
  637. mpd_qmax(result, a, b, ctx, &status);
  638. mpd_addstatus_raise(ctx, status);
  639. }
  640. void
  641. mpd_max_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  642. {
  643. uint32_t status = 0;
  644. mpd_qmax_mag(result, a, b, ctx, &status);
  645. mpd_addstatus_raise(ctx, status);
  646. }
  647. void
  648. mpd_min(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  649. {
  650. uint32_t status = 0;
  651. mpd_qmin(result, a, b, ctx, &status);
  652. mpd_addstatus_raise(ctx, status);
  653. }
  654. void
  655. mpd_min_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  656. {
  657. uint32_t status = 0;
  658. mpd_qmin_mag(result, a, b, ctx, &status);
  659. mpd_addstatus_raise(ctx, status);
  660. }
  661. void
  662. mpd_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  663. {
  664. uint32_t status = 0;
  665. mpd_qminus(result, a, ctx, &status);
  666. mpd_addstatus_raise(ctx, status);
  667. }
  668. void
  669. mpd_mul(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  670. {
  671. uint32_t status = 0;
  672. mpd_qmul(result, a, b, ctx, &status);
  673. mpd_addstatus_raise(ctx, status);
  674. }
  675. void
  676. mpd_mul_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx)
  677. {
  678. uint32_t status = 0;
  679. mpd_qmul_ssize(result, a, b, ctx, &status);
  680. mpd_addstatus_raise(ctx, status);
  681. }
  682. void
  683. mpd_mul_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx)
  684. {
  685. uint32_t status = 0;
  686. mpd_qmul_i32(result, a, b, ctx, &status);
  687. mpd_addstatus_raise(ctx, status);
  688. }
  689. #ifndef LEGACY_COMPILER
  690. void
  691. mpd_mul_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx)
  692. {
  693. uint32_t status = 0;
  694. mpd_qmul_i64(result, a, b, ctx, &status);
  695. mpd_addstatus_raise(ctx, status);
  696. }
  697. #endif
  698. void
  699. mpd_mul_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx)
  700. {
  701. uint32_t status = 0;
  702. mpd_qmul_uint(result, a, b, ctx, &status);
  703. mpd_addstatus_raise(ctx, status);
  704. }
  705. void
  706. mpd_mul_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx)
  707. {
  708. uint32_t status = 0;
  709. mpd_qmul_u32(result, a, b, ctx, &status);
  710. mpd_addstatus_raise(ctx, status);
  711. }
  712. #ifndef LEGACY_COMPILER
  713. void
  714. mpd_mul_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx)
  715. {
  716. uint32_t status = 0;
  717. mpd_qmul_u64(result, a, b, ctx, &status);
  718. mpd_addstatus_raise(ctx, status);
  719. }
  720. #endif
  721. void
  722. mpd_next_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  723. {
  724. uint32_t status = 0;
  725. mpd_qnext_minus(result, a, ctx, &status);
  726. mpd_addstatus_raise(ctx, status);
  727. }
  728. void
  729. mpd_next_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  730. {
  731. uint32_t status = 0;
  732. mpd_qnext_plus(result, a, ctx, &status);
  733. mpd_addstatus_raise(ctx, status);
  734. }
  735. void
  736. mpd_next_toward(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  737. {
  738. uint32_t status = 0;
  739. mpd_qnext_toward(result, a, b, ctx, &status);
  740. mpd_addstatus_raise(ctx, status);
  741. }
  742. void
  743. mpd_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  744. {
  745. uint32_t status = 0;
  746. mpd_qplus(result, a, ctx, &status);
  747. mpd_addstatus_raise(ctx, status);
  748. }
  749. void
  750. mpd_pow(mpd_t *result, const mpd_t *base, const mpd_t *exp, mpd_context_t *ctx)
  751. {
  752. uint32_t status = 0;
  753. mpd_qpow(result, base, exp, ctx, &status);
  754. mpd_addstatus_raise(ctx, status);
  755. }
  756. void
  757. mpd_powmod(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_t *mod,
  758. mpd_context_t *ctx)
  759. {
  760. uint32_t status = 0;
  761. mpd_qpowmod(result, base, exp, mod, ctx, &status);
  762. mpd_addstatus_raise(ctx, status);
  763. }
  764. void
  765. mpd_quantize(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  766. {
  767. uint32_t status = 0;
  768. mpd_qquantize(result, a, b, ctx, &status);
  769. mpd_addstatus_raise(ctx, status);
  770. }
  771. void
  772. mpd_rescale(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, mpd_context_t *ctx)
  773. {
  774. uint32_t status = 0;
  775. mpd_qrescale(result, a, exp, ctx, &status);
  776. mpd_addstatus_raise(ctx, status);
  777. }
  778. void
  779. mpd_reduce(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  780. {
  781. uint32_t status = 0;
  782. mpd_qreduce(result, a, ctx, &status);
  783. mpd_addstatus_raise(ctx, status);
  784. }
  785. void
  786. mpd_rem(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  787. {
  788. uint32_t status = 0;
  789. mpd_qrem(r, a, b, ctx, &status);
  790. mpd_addstatus_raise(ctx, status);
  791. }
  792. void
  793. mpd_rem_near(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx)
  794. {
  795. uint32_t status = 0;
  796. mpd_qrem_near(r, a, b, ctx, &status);
  797. mpd_addstatus_raise(ctx, status);
  798. }
  799. void
  800. mpd_round_to_intx(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  801. {
  802. uint32_t status = 0;
  803. mpd_qround_to_intx(result, a, ctx, &status);
  804. mpd_addstatus_raise(ctx, status);
  805. }
  806. void
  807. mpd_round_to_int(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  808. {
  809. uint32_t status = 0;
  810. mpd_qround_to_int(result, a, ctx, &status);
  811. mpd_addstatus_raise(ctx, status);
  812. }
  813. void
  814. mpd_trunc(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  815. {
  816. uint32_t status = 0;
  817. mpd_qtrunc(result, a, ctx, &status);
  818. mpd_addstatus_raise(ctx, status);
  819. }
  820. void
  821. mpd_floor(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  822. {
  823. uint32_t status = 0;
  824. mpd_qfloor(result, a, ctx, &status);
  825. mpd_addstatus_raise(ctx, status);
  826. }
  827. void
  828. mpd_ceil(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  829. {
  830. uint32_t status = 0;
  831. mpd_qceil(result, a, ctx, &status);
  832. mpd_addstatus_raise(ctx, status);
  833. }
  834. void
  835. mpd_sqrt(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  836. {
  837. uint32_t status = 0;
  838. mpd_qsqrt(result, a, ctx, &status);
  839. mpd_addstatus_raise(ctx, status);
  840. }
  841. void
  842. mpd_invroot(mpd_t *result, const mpd_t *a, mpd_context_t *ctx)
  843. {
  844. uint32_t status = 0;
  845. mpd_qinvroot(result, a, ctx, &status);
  846. mpd_addstatus_raise(ctx, status);
  847. }