123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349 |
- /*
- * Copyright (c) 2008-2020 Stefan Krah. All rights reserved.
- *
- * Redistribution and use in source and binary forms, with or without
- * modification, are permitted provided that the following conditions
- * are met:
- *
- * 1. Redistributions of source code must retain the above copyright
- * notice, this list of conditions and the following disclaimer.
- *
- * 2. Redistributions in binary form must reproduce the above copyright
- * notice, this list of conditions and the following disclaimer in the
- * documentation and/or other materials provided with the distribution.
- *
- * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS" AND
- * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
- * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
- * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
- * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
- * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
- * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
- * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
- * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
- * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
- * SUCH DAMAGE.
- */
- #include "mpdecimal.h"
- #include <assert.h>
- #include <stdio.h>
- #include <stdlib.h>
- #include <string.h>
- #include "mpalloc.h"
- #include "typearith.h"
- #if defined(_MSC_VER)
- #pragma warning(disable : 4232)
- #endif
- /* Guaranteed minimum allocation for a coefficient. May be changed once
- at program start using mpd_setminalloc(). */
- mpd_ssize_t MPD_MINALLOC = MPD_MINALLOC_MIN;
- /* Custom allocation and free functions */
- void *(* mpd_mallocfunc)(size_t size) = malloc;
- void *(* mpd_reallocfunc)(void *ptr, size_t size) = realloc;
- void *(* mpd_callocfunc)(size_t nmemb, size_t size) = calloc;
- void (* mpd_free)(void *ptr) = free;
- /* emulate calloc if it is not available */
- void *
- mpd_callocfunc_em(size_t nmemb, size_t size)
- {
- void *ptr;
- size_t req;
- mpd_size_t overflow;
- req = mul_size_t_overflow((mpd_size_t)nmemb, (mpd_size_t)size,
- &overflow);
- if (overflow) {
- return NULL;
- }
- ptr = mpd_mallocfunc(req);
- if (ptr == NULL) {
- return NULL;
- }
- /* used on uint32_t or uint64_t */
- memset(ptr, 0, req);
- return ptr;
- }
- /* malloc with overflow checking */
- void *
- mpd_alloc(mpd_size_t nmemb, mpd_size_t size)
- {
- mpd_size_t req, overflow;
- req = mul_size_t_overflow(nmemb, size, &overflow);
- if (overflow) {
- return NULL;
- }
- return mpd_mallocfunc(req);
- }
- /* calloc with overflow checking */
- void *
- mpd_calloc(mpd_size_t nmemb, mpd_size_t size)
- {
- mpd_size_t overflow;
- (void)mul_size_t_overflow(nmemb, size, &overflow);
- if (overflow) {
- return NULL;
- }
- return mpd_callocfunc(nmemb, size);
- }
- /* realloc with overflow checking */
- void *
- mpd_realloc(void *ptr, mpd_size_t nmemb, mpd_size_t size, uint8_t *err)
- {
- void *new;
- mpd_size_t req, overflow;
- req = mul_size_t_overflow(nmemb, size, &overflow);
- if (overflow) {
- *err = 1;
- return ptr;
- }
- new = mpd_reallocfunc(ptr, req);
- if (new == NULL) {
- *err = 1;
- return ptr;
- }
- return new;
- }
- /* struct hack malloc with overflow checking */
- void *
- mpd_sh_alloc(mpd_size_t struct_size, mpd_size_t nmemb, mpd_size_t size)
- {
- mpd_size_t req, overflow;
- req = mul_size_t_overflow(nmemb, size, &overflow);
- if (overflow) {
- return NULL;
- }
- req = add_size_t_overflow(req, struct_size, &overflow);
- if (overflow) {
- return NULL;
- }
- return mpd_mallocfunc(req);
- }
- /* Allocate a new decimal with a coefficient of length 'nwords'. In case
- of an error the return value is NULL. */
- mpd_t *
- mpd_qnew_size(mpd_ssize_t nwords)
- {
- mpd_t *result;
- nwords = (nwords < MPD_MINALLOC) ? MPD_MINALLOC : nwords;
- result = mpd_alloc(1, sizeof *result);
- if (result == NULL) {
- return NULL;
- }
- result->data = mpd_alloc(nwords, sizeof *result->data);
- if (result->data == NULL) {
- mpd_free(result);
- return NULL;
- }
- result->flags = 0;
- result->exp = 0;
- result->digits = 0;
- result->len = 0;
- result->alloc = nwords;
- return result;
- }
- /* Allocate a new decimal with a coefficient of length MPD_MINALLOC.
- In case of an error the return value is NULL. */
- mpd_t *
- mpd_qnew(void)
- {
- return mpd_qnew_size(MPD_MINALLOC);
- }
- /* Allocate new decimal. Caller can check for NULL or MPD_Malloc_error.
- Raises on error. */
- mpd_t *
- mpd_new(mpd_context_t *ctx)
- {
- mpd_t *result;
- result = mpd_qnew();
- if (result == NULL) {
- mpd_addstatus_raise(ctx, MPD_Malloc_error);
- }
- return result;
- }
- /*
- * Input: 'result' is a static mpd_t with a static coefficient.
- * Assumption: 'nwords' >= result->alloc.
- *
- * Resize the static coefficient to a larger dynamic one and copy the
- * existing data. If successful, the value of 'result' is unchanged.
- * Otherwise, set 'result' to NaN and update 'status' with MPD_Malloc_error.
- */
- int
- mpd_switch_to_dyn(mpd_t *result, mpd_ssize_t nwords, uint32_t *status)
- {
- mpd_uint_t *p = result->data;
- assert(nwords >= result->alloc);
- result->data = mpd_alloc(nwords, sizeof *result->data);
- if (result->data == NULL) {
- result->data = p;
- mpd_set_qnan(result);
- mpd_set_positive(result);
- result->exp = result->digits = result->len = 0;
- *status |= MPD_Malloc_error;
- return 0;
- }
- memcpy(result->data, p, result->alloc * (sizeof *result->data));
- result->alloc = nwords;
- mpd_set_dynamic_data(result);
- return 1;
- }
- /*
- * Input: 'result' is a static mpd_t with a static coefficient.
- *
- * Convert the coefficient to a dynamic one that is initialized to zero. If
- * malloc fails, set 'result' to NaN and update 'status' with MPD_Malloc_error.
- */
- int
- mpd_switch_to_dyn_zero(mpd_t *result, mpd_ssize_t nwords, uint32_t *status)
- {
- mpd_uint_t *p = result->data;
- result->data = mpd_calloc(nwords, sizeof *result->data);
- if (result->data == NULL) {
- result->data = p;
- mpd_set_qnan(result);
- mpd_set_positive(result);
- result->exp = result->digits = result->len = 0;
- *status |= MPD_Malloc_error;
- return 0;
- }
- result->alloc = nwords;
- mpd_set_dynamic_data(result);
- return 1;
- }
- /*
- * Input: 'result' is a static or a dynamic mpd_t with a dynamic coefficient.
- * Resize the coefficient to length 'nwords':
- * Case nwords > result->alloc:
- * If realloc is successful:
- * 'result' has a larger coefficient but the same value. Return 1.
- * Otherwise:
- * Set 'result' to NaN, update status with MPD_Malloc_error and return 0.
- * Case nwords < result->alloc:
- * If realloc is successful:
- * 'result' has a smaller coefficient. result->len is undefined. Return 1.
- * Otherwise (unlikely):
- * 'result' is unchanged. Reuse the now oversized coefficient. Return 1.
- */
- int
- mpd_realloc_dyn(mpd_t *result, mpd_ssize_t nwords, uint32_t *status)
- {
- uint8_t err = 0;
- result->data = mpd_realloc(result->data, nwords, sizeof *result->data, &err);
- if (!err) {
- result->alloc = nwords;
- }
- else if (nwords > result->alloc) {
- mpd_set_qnan(result);
- mpd_set_positive(result);
- result->exp = result->digits = result->len = 0;
- *status |= MPD_Malloc_error;
- return 0;
- }
- return 1;
- }
- /*
- * Input: 'result' is a static mpd_t with a static coefficient.
- * Assumption: 'nwords' >= result->alloc.
- *
- * Resize the static coefficient to a larger dynamic one and copy the
- * existing data.
- *
- * On failure the value of 'result' is unchanged.
- */
- int
- mpd_switch_to_dyn_cxx(mpd_t *result, mpd_ssize_t nwords)
- {
- assert(nwords >= result->alloc);
- mpd_uint_t *data = mpd_alloc(nwords, sizeof *result->data);
- if (data == NULL) {
- return 0;
- }
- memcpy(data, result->data, result->alloc * (sizeof *result->data));
- result->data = data;
- result->alloc = nwords;
- mpd_set_dynamic_data(result);
- return 1;
- }
- /*
- * Input: 'result' is a static or a dynamic mpd_t with a dynamic coefficient.
- * Resize the coefficient to length 'nwords':
- * Case nwords > result->alloc:
- * If realloc is successful:
- * 'result' has a larger coefficient but the same value. Return 1.
- * Otherwise:
- * 'result' has a the same coefficient. Return 0.
- * Case nwords < result->alloc:
- * If realloc is successful:
- * 'result' has a smaller coefficient. result->len is undefined. Return 1.
- * Otherwise (unlikely):
- * 'result' is unchanged. Reuse the now oversized coefficient. Return 1.
- */
- int
- mpd_realloc_dyn_cxx(mpd_t *result, mpd_ssize_t nwords)
- {
- uint8_t err = 0;
- mpd_uint_t *p = mpd_realloc(result->data, nwords, sizeof *result->data, &err);
- if (!err) {
- result->data = p;
- result->alloc = nwords;
- }
- else if (nwords > result->alloc) {
- return 0;
- }
- return 1;
- }
|