cvector_d.c

Go to the documentation of this file.

#include "cvector_d.h"
 
 
 
#if defined(CVEC_MALLOC) && defined(CVEC_FREE) && defined(CVEC_REALLOC)
/* ok */
#elif !defined(CVEC_MALLOC) && !defined(CVEC_FREE) && !defined(CVEC_REALLOC)
/* ok */
#else
#error "Must define all or none of CVEC_MALLOC, CVEC_FREE, and CVEC_REALLOC."
#endif
 
#ifndef CVEC_MALLOC
#define CVEC_MALLOC(sz)      malloc(sz)
#define CVEC_REALLOC(p, sz)  realloc(p, sz)
#define CVEC_FREE(p)         free(p)
#endif
 
#ifndef CVEC_MEMMOVE
#include <string.h>
#define CVEC_MEMMOVE(dst, src, sz)  memmove(dst, src, sz)
#endif
 
#ifndef CVEC_ASSERT
#include <assert.h>
#define CVEC_ASSERT(x)       assert(x)
#endif
 
cvec_sz CVEC_D_START_SZ = 50;
 
#define CVEC_D_ALLOCATOR(x) ((x+1) * 2)
 
cvector_d* cvec_d_heap(cvec_sz size, cvec_sz capacity)
{
    cvector_d* vec;
    
    if (!(vec = (cvector_d*)CVEC_MALLOC(sizeof(cvector_d)))) {
        CVEC_ASSERT(vec != NULL);
        return NULL;
    }
 
    vec->size = size;
    vec->capacity = (capacity > vec->size || (vec->size && capacity == vec->size)) ? capacity : vec->size + CVEC_D_START_SZ;
 
    if (!(vec->a = (double*)CVEC_MALLOC(vec->capacity*sizeof(double)))) {
        CVEC_ASSERT(vec->a != NULL);
        CVEC_FREE(vec);
        return NULL;
    }
 
    return vec;
}
 
cvector_d* cvec_init_d_heap(double* vals, cvec_sz num)
{
    cvector_d* vec;
    
    if (!(vec = (cvector_d*)CVEC_MALLOC(sizeof(cvector_d)))) {
        CVEC_ASSERT(vec != NULL);
        return NULL;
    }
 
    vec->capacity = num + CVEC_D_START_SZ;
    vec->size = num;
    if (!(vec->a = (double*)CVEC_MALLOC(vec->capacity*sizeof(double)))) {
        CVEC_ASSERT(vec->a != NULL);
        CVEC_FREE(vec);
        return NULL;
    }
 
    CVEC_MEMMOVE(vec->a, vals, sizeof(double)*num);
 
    return vec;
}
 
int cvec_d(cvector_d* vec, cvec_sz size, cvec_sz capacity)
{
    vec->size = size;
    vec->capacity = (capacity > vec->size || (vec->size && capacity == vec->size)) ? capacity : vec->size + CVEC_D_START_SZ;
 
    if (!(vec->a = (double*)CVEC_MALLOC(vec->capacity*sizeof(double)))) {
        CVEC_ASSERT(vec->a != NULL);
        vec->size = vec->capacity = 0;
        return 0;
    }
 
    return 1;
}
 
int cvec_init_d(cvector_d* vec, double* vals, cvec_sz num)
{
    vec->capacity = num + CVEC_D_START_SZ;
    vec->size = num;
    if (!(vec->a = (double*)CVEC_MALLOC(vec->capacity*sizeof(double)))) {
        CVEC_ASSERT(vec->a != NULL);
        vec->size = vec->capacity = 0;
        return 0;
    }
 
    CVEC_MEMMOVE(vec->a, vals, sizeof(double)*num);
 
    return 1;
}
 
int cvec_copyc_d(void* dest, void* src)
{
    cvector_d* vec1 = (cvector_d*)dest;
    cvector_d* vec2 = (cvector_d*)src;
 
    vec1->a = NULL;
    vec1->size = 0;
    vec1->capacity = 0;
 
    return cvec_copy_d(vec1, vec2);
}
 
int cvec_copy_d(cvector_d* dest, cvector_d* src)
{
    double* tmp = NULL;
    if (!(tmp = (double*)CVEC_REALLOC(dest->a, src->capacity*sizeof(double)))) {
        CVEC_ASSERT(tmp != NULL);
        return 0;
    }
    dest->a = tmp;
    
    CVEC_MEMMOVE(dest->a, src->a, src->size*sizeof(double));
    dest->size = src->size;
    dest->capacity = src->capacity;
    return 1;
}
 
 
int cvec_push_d(cvector_d* vec, double a)
{
    double* tmp;
    cvec_sz tmp_sz;
    if (vec->capacity == vec->size) {
        tmp_sz = CVEC_D_ALLOCATOR(vec->capacity);
        if (!(tmp = (double*)CVEC_REALLOC(vec->a, sizeof(double)*tmp_sz))) {
            CVEC_ASSERT(tmp != NULL);
            return 0;
        }
        vec->a = tmp;
        vec->capacity = tmp_sz;
    }
    vec->a[vec->size++] = a;
    return 1;
}
 
double cvec_pop_d(cvector_d* vec)
{
    return vec->a[--vec->size];
}
 
double* cvec_back_d(cvector_d* vec)
{
    return &vec->a[vec->size-1];
}
 
int cvec_extend_d(cvector_d* vec, cvec_sz num)
{
    double* tmp;
    cvec_sz tmp_sz;
    if (vec->capacity < vec->size + num) {
        tmp_sz = vec->capacity + num + CVEC_D_START_SZ;
        if (!(tmp = (double*)CVEC_REALLOC(vec->a, sizeof(double)*tmp_sz))) {
            CVEC_ASSERT(tmp != NULL);
            return 0;
        }
        vec->a = tmp;
        vec->capacity = tmp_sz;
    }
 
    vec->size += num;
    return 1;
}
 
int cvec_insert_d(cvector_d* vec, cvec_sz i, double a)
{
    double* tmp;
    cvec_sz tmp_sz;
    if (vec->capacity == vec->size) {
        tmp_sz = CVEC_D_ALLOCATOR(vec->capacity);
        if (!(tmp = (double*)CVEC_REALLOC(vec->a, sizeof(double)*tmp_sz))) {
            CVEC_ASSERT(tmp != NULL);
            return 0;
        }
        vec->a = tmp;
        vec->capacity = tmp_sz;
    }
    
    CVEC_MEMMOVE(&vec->a[i+1], &vec->a[i], (vec->size-i)*sizeof(double));
    vec->a[i] = a;
    vec->size++;
    return 1;
}
 
int cvec_insert_array_d(cvector_d* vec, cvec_sz i, double* a, cvec_sz num)
{
    double* tmp;
    cvec_sz tmp_sz;
    if (vec->capacity < vec->size + num) {
        tmp_sz = vec->capacity + num + CVEC_D_START_SZ;
        if (!(tmp = (double*)CVEC_REALLOC(vec->a, sizeof(double)*tmp_sz))) {
            CVEC_ASSERT(tmp != NULL);
            return 0;
        }
        vec->a = tmp;
        vec->capacity = tmp_sz;
    }
 
    CVEC_MEMMOVE(&vec->a[i+num], &vec->a[i], (vec->size-i)*sizeof(double));
    CVEC_MEMMOVE(&vec->a[i], a, num*sizeof(double));
    vec->size += num;
    return 1;
}
 
double cvec_replace_d(cvector_d* vec, cvec_sz i, double a)
{
    double tmp = vec->a[i];
    vec->a[i] = a;
    return tmp;
}
 
void cvec_erase_d(cvector_d* vec, cvec_sz start, cvec_sz end)
{
    cvec_sz d = end - start + 1;
    CVEC_MEMMOVE(&vec->a[start], &vec->a[end+1], (vec->size-1-end)*sizeof(double));
    vec->size -= d;
}
 
int cvec_reserve_d(cvector_d* vec, cvec_sz size)
{
    double* tmp;
    if (vec->capacity < size) {
        if (!(tmp = (double*)CVEC_REALLOC(vec->a, sizeof(double)*(size+CVEC_D_START_SZ)))) {
            CVEC_ASSERT(tmp != NULL);
            return 0;
        }
        vec->a = tmp;
        vec->capacity = size + CVEC_D_START_SZ;
    }
    return 1;
}
 
int cvec_set_cap_d(cvector_d* vec, cvec_sz size)
{
    double* tmp;
    if (size < vec->size)
        vec->size = size;
 
    if (!(tmp = (double*)CVEC_REALLOC(vec->a, sizeof(double)*size))) {
        CVEC_ASSERT(tmp != NULL);
        return 0;
    }
    vec->a = tmp;
    vec->capacity = size;
    return 1;
}
 
void cvec_set_val_sz_d(cvector_d* vec, double val)
{
    cvec_sz i;
    for(i=0; i<vec->size; i++) {
        vec->a[i] = val;
    }
}
 
void cvec_set_val_cap_d(cvector_d* vec, double val)
{
    cvec_sz i;
    for(i=0; i<vec->capacity; i++) {
        vec->a[i] = val;
    }
}
 
void cvec_clear_d(cvector_d* vec) { vec->size = 0; }
 
void cvec_free_d_heap(void* vec)
{
    cvector_d* tmp = (cvector_d*)vec;
    if (!tmp) return;
    CVEC_FREE(tmp->a);
    CVEC_FREE(tmp);
}
 
void cvec_free_d(void* vec)
{
    cvector_d* tmp = (cvector_d*)vec;
    CVEC_FREE(tmp->a); tmp->size = 0;
    tmp->capacity = 0;
}