mirror of
https://git.savannah.gnu.org/git/guile.git
synced 2025-05-01 04:10:18 +02:00
* srfi-4.h, srfi-4.c (scm_uniform_vector_p, scm_uniform_vector_ref, scm_uniform_vector_set_x, scm_uniform_vector_to_list, scm_is_uniform_vector, scm_c_uniform_vector_lengths, scm_c_uniform_vector_size, scm_uniform_vector_elements, scm_uniform_vector_element_size, scm_uniform_vector_release): New. (scm_i_uniform_vector_prototype, scm_i_uniform_vector_tag): New. (scm_uniform_element_size, scm_uniform_vector_length): Moved here from unif.h, unif.c and extended to handle both the old and new uniform vectors. * unif.h, unif.c (scm_uniform_vector_ref, scm_array_ref): Renamed the former to the latter. (scm_uniform_vector_length, scm_uniform_element_size): Moved to srfi-4.h, srfi-4.c. (scm_make_uve): Call scm_make_s8vector for #\nul prototype. (scm_array_p, scm_array_rank, scm_array_dimensions, scm_transpose_array, scm_enclose_array, scm_array_ref, scm_cvref, scm_array_set_x, scm_array_contents, scm_uniform_array_read_x, scm_array_to_list, scm_array_prototype): Handle srfi-4 uniform vectors. Removed code for scm_tc7_byvect. (scm_dimensions_to_uniform_array): Fill array with 0 when prototype is #\nul. (scm_i_print_array_dimension, scm_i_legacy_tag, scm_i_print_array): New. (scm_raprin1): Call scm_i_print_array for arrays. Removed code for scm_tc7_byvect.
735 lines
19 KiB
C
735 lines
19 KiB
C
/* srfi-4.c --- Homogeneous numeric vector datatypes.
|
|
*
|
|
* Copyright (C) 2001, 2004 Free Software Foundation, Inc.
|
|
*
|
|
* This library is free software; you can redistribute it and/or
|
|
* modify it under the terms of the GNU Lesser General Public
|
|
* License as published by the Free Software Foundation; either
|
|
* version 2.1 of the License, or (at your option) any later version.
|
|
*
|
|
* This library is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
* Lesser General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU Lesser General Public
|
|
* License along with this library; if not, write to the Free Software
|
|
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
|
|
*/
|
|
|
|
#include <libguile.h>
|
|
#include <string.h>
|
|
#include <stdio.h>
|
|
|
|
#include "libguile/srfi-4.h"
|
|
#include "libguile/error.h"
|
|
#include "libguile/read.h"
|
|
#include "libguile/ports.h"
|
|
#include "libguile/chars.h"
|
|
|
|
/* Smob type code for homogeneous numeric vectors. */
|
|
int scm_tc16_uvec = 0;
|
|
|
|
|
|
/* Accessor macros for the three components of a homogeneous numeric
|
|
vector:
|
|
- The type tag (one of the symbolic constants below).
|
|
- The vector's length (counted in elements).
|
|
- The address of the data area (holding the elements of the
|
|
vector). */
|
|
#define SCM_UVEC_TYPE(u) (SCM_CELL_WORD_1(u))
|
|
#define SCM_UVEC_LENGTH(u) ((size_t)SCM_CELL_WORD_2(u))
|
|
#define SCM_UVEC_BASE(u) ((void *)SCM_CELL_WORD_3(u))
|
|
|
|
|
|
/* Symbolic constants encoding the various types of homogeneous
|
|
numeric vectors. */
|
|
#define SCM_UVEC_U8 0
|
|
#define SCM_UVEC_S8 1
|
|
#define SCM_UVEC_U16 2
|
|
#define SCM_UVEC_S16 3
|
|
#define SCM_UVEC_U32 4
|
|
#define SCM_UVEC_S32 5
|
|
#define SCM_UVEC_U64 6
|
|
#define SCM_UVEC_S64 7
|
|
#define SCM_UVEC_F32 8
|
|
#define SCM_UVEC_F64 9
|
|
|
|
|
|
/* This array maps type tags to the size of the elements. */
|
|
static const int uvec_sizes[10] = {
|
|
1, 1,
|
|
2, 2,
|
|
4, 4,
|
|
8, 8,
|
|
sizeof(float), sizeof(double)
|
|
};
|
|
|
|
static const char *uvec_tags[10] = {
|
|
"u8", "s8",
|
|
"u16", "s16",
|
|
"u32", "s32",
|
|
"u64", "s64",
|
|
"f32", "f64"
|
|
};
|
|
|
|
static const char *uvec_names[10] = {
|
|
"u8vector", "s8vector",
|
|
"u16vector", "s16vector",
|
|
"u32vector", "s32vector",
|
|
"u64vector", "s64vector",
|
|
"f32vector", "f64vector"
|
|
};
|
|
|
|
/* ================================================================ */
|
|
/* SMOB procedures. */
|
|
/* ================================================================ */
|
|
|
|
|
|
/* Smob print hook for homogeneous vectors. */
|
|
static int
|
|
uvec_print (SCM uvec, SCM port, scm_print_state *pstate)
|
|
{
|
|
union {
|
|
scm_t_uint8 *u8;
|
|
scm_t_int8 *s8;
|
|
scm_t_uint16 *u16;
|
|
scm_t_int16 *s16;
|
|
scm_t_uint32 *u32;
|
|
scm_t_int32 *s32;
|
|
#if SCM_HAVE_T_INT64
|
|
scm_t_uint64 *u64;
|
|
scm_t_int64 *s64;
|
|
#endif
|
|
float *f32;
|
|
double *f64;
|
|
} np;
|
|
|
|
size_t i = 0;
|
|
const size_t uvlen = SCM_UVEC_LENGTH (uvec);
|
|
void *uptr = SCM_UVEC_BASE (uvec);
|
|
|
|
switch (SCM_UVEC_TYPE (uvec))
|
|
{
|
|
case SCM_UVEC_U8: np.u8 = (scm_t_uint8 *) uptr; break;
|
|
case SCM_UVEC_S8: np.s8 = (scm_t_int8 *) uptr; break;
|
|
case SCM_UVEC_U16: np.u16 = (scm_t_uint16 *) uptr; break;
|
|
case SCM_UVEC_S16: np.s16 = (scm_t_int16 *) uptr; break;
|
|
case SCM_UVEC_U32: np.u32 = (scm_t_uint32 *) uptr; break;
|
|
case SCM_UVEC_S32: np.s32 = (scm_t_int32 *) uptr; break;
|
|
#if SCM_HAVE_T_INT64
|
|
case SCM_UVEC_U64: np.u64 = (scm_t_uint64 *) uptr; break;
|
|
case SCM_UVEC_S64: np.s64 = (scm_t_int64 *) uptr; break;
|
|
#endif
|
|
case SCM_UVEC_F32: np.f32 = (float *) uptr; break;
|
|
case SCM_UVEC_F64: np.f64 = (double *) uptr; break;
|
|
default:
|
|
abort (); /* Sanity check. */
|
|
break;
|
|
}
|
|
|
|
scm_putc ('#', port);
|
|
scm_puts (uvec_tags [SCM_UVEC_TYPE (uvec)], port);
|
|
scm_putc ('(', port);
|
|
|
|
while (i < uvlen)
|
|
{
|
|
if (i != 0) scm_puts (" ", port);
|
|
switch (SCM_UVEC_TYPE (uvec))
|
|
{
|
|
case SCM_UVEC_U8: scm_uintprint (*np.u8, 10, port); np.u8++; break;
|
|
case SCM_UVEC_S8: scm_intprint (*np.s8, 10, port); np.s8++; break;
|
|
case SCM_UVEC_U16: scm_uintprint (*np.u16, 10, port); np.u16++; break;
|
|
case SCM_UVEC_S16: scm_intprint (*np.s16, 10, port); np.s16++; break;
|
|
case SCM_UVEC_U32: scm_uintprint (*np.u32, 10, port); np.u32++; break;
|
|
case SCM_UVEC_S32: scm_intprint (*np.s32, 10, port); np.s32++; break;
|
|
#if SCM_HAVE_T_INT64
|
|
case SCM_UVEC_U64: scm_uintprint (*np.u64, 10, port); np.u64++; break;
|
|
case SCM_UVEC_S64: scm_intprint (*np.s64, 10, port); np.s64++; break;
|
|
#endif
|
|
case SCM_UVEC_F32: scm_i_print_double (*np.f32, port); np.f32++; break;
|
|
case SCM_UVEC_F64: scm_i_print_double (*np.f64, port); np.f64++; break;
|
|
default:
|
|
abort (); /* Sanity check. */
|
|
break;
|
|
}
|
|
i++;
|
|
}
|
|
scm_remember_upto_here_1 (uvec);
|
|
scm_puts (")", port);
|
|
return 1;
|
|
}
|
|
|
|
const char *
|
|
scm_i_uniform_vector_tag (SCM uvec)
|
|
{
|
|
return uvec_tags[SCM_UVEC_TYPE (uvec)];
|
|
}
|
|
|
|
static SCM
|
|
uvec_equalp (SCM a, SCM b)
|
|
{
|
|
SCM result = SCM_BOOL_T;
|
|
if (SCM_UVEC_TYPE (a) != SCM_UVEC_TYPE (b))
|
|
result = SCM_BOOL_F;
|
|
else if (SCM_UVEC_LENGTH (a) != SCM_UVEC_LENGTH (b))
|
|
result = SCM_BOOL_F;
|
|
else if (memcmp (SCM_UVEC_BASE (a), SCM_UVEC_BASE (b),
|
|
SCM_UVEC_LENGTH (a) * uvec_sizes[SCM_UVEC_TYPE(a)]) != 0)
|
|
result = SCM_BOOL_F;
|
|
|
|
scm_remember_upto_here_2 (a, b);
|
|
return result;
|
|
}
|
|
|
|
/* Smob free hook for homogeneous numeric vectors. */
|
|
static size_t
|
|
uvec_free (SCM uvec)
|
|
{
|
|
int type = SCM_UVEC_TYPE (uvec);
|
|
scm_gc_free (SCM_UVEC_BASE (uvec),
|
|
SCM_UVEC_LENGTH (uvec) * uvec_sizes[type],
|
|
uvec_names[type]);
|
|
return 0;
|
|
}
|
|
|
|
/* ================================================================ */
|
|
/* Utility procedures. */
|
|
/* ================================================================ */
|
|
|
|
static SCM_C_INLINE int
|
|
is_uvec (int type, SCM obj)
|
|
{
|
|
return (SCM_SMOB_PREDICATE (scm_tc16_uvec, obj)
|
|
&& SCM_UVEC_TYPE (obj) == type);
|
|
}
|
|
|
|
static SCM_C_INLINE SCM
|
|
uvec_p (int type, SCM obj)
|
|
{
|
|
return scm_from_bool (is_uvec (type, obj));
|
|
}
|
|
|
|
static SCM_C_INLINE void
|
|
uvec_assert (int type, SCM obj)
|
|
{
|
|
if (!is_uvec (type, obj))
|
|
scm_wrong_type_arg_msg (NULL, 0, obj, uvec_names[type]);
|
|
}
|
|
|
|
/* Create a new, uninitialized homogeneous numeric vector of type TYPE
|
|
with space for LEN elements. */
|
|
static SCM
|
|
alloc_uvec (int type, size_t c_len)
|
|
{
|
|
void *base = scm_gc_malloc (c_len * uvec_sizes[type], uvec_names[type]);
|
|
SCM_RETURN_NEWSMOB3 (scm_tc16_uvec, type, c_len, (scm_t_bits) base);
|
|
}
|
|
|
|
/* GCC doesn't seem to want to optimize unused switch clauses away,
|
|
so we use a big 'if' in the next two functions.
|
|
*/
|
|
|
|
static SCM_C_INLINE SCM
|
|
uvec_fast_ref (int type, void *base, size_t c_idx)
|
|
{
|
|
if (type == SCM_UVEC_U8)
|
|
return scm_from_uint8 (((scm_t_uint8*)base)[c_idx]);
|
|
else if (type == SCM_UVEC_S8)
|
|
return scm_from_int8 (((scm_t_int8*)base)[c_idx]);
|
|
else if (type == SCM_UVEC_U16)
|
|
return scm_from_uint16 (((scm_t_uint16*)base)[c_idx]);
|
|
else if (type == SCM_UVEC_S16)
|
|
return scm_from_int16 (((scm_t_int16*)base)[c_idx]);
|
|
else if (type == SCM_UVEC_U32)
|
|
return scm_from_uint32 (((scm_t_uint32*)base)[c_idx]);
|
|
else if (type == SCM_UVEC_S32)
|
|
return scm_from_int32 (((scm_t_int32*)base)[c_idx]);
|
|
#if SCM_HAVE_T_INT64
|
|
else if (type == SCM_UVEC_U64)
|
|
return scm_from_uint64 (((scm_t_uint64*)base)[c_idx]);
|
|
else if (type == SCM_UVEC_S64)
|
|
return scm_from_int64 (((scm_t_int64*)base)[c_idx]);
|
|
#endif
|
|
else if (type == SCM_UVEC_F32)
|
|
return scm_from_double (((float*)base)[c_idx]);
|
|
else if (type == SCM_UVEC_F64)
|
|
return scm_from_double (((double*)base)[c_idx]);
|
|
}
|
|
|
|
static SCM_C_INLINE void
|
|
uvec_fast_set_x (int type, void *base, size_t c_idx, SCM val)
|
|
{
|
|
if (type == SCM_UVEC_U8)
|
|
(((scm_t_uint8*)base)[c_idx]) = scm_to_uint8 (val);
|
|
else if (type == SCM_UVEC_S8)
|
|
(((scm_t_int8*)base)[c_idx]) = scm_to_int8 (val);
|
|
else if (type == SCM_UVEC_U16)
|
|
(((scm_t_uint16*)base)[c_idx]) = scm_to_uint16 (val);
|
|
else if (type == SCM_UVEC_S16)
|
|
(((scm_t_int16*)base)[c_idx]) = scm_to_int16 (val);
|
|
else if (type == SCM_UVEC_U32)
|
|
(((scm_t_uint32*)base)[c_idx]) = scm_to_uint32 (val);
|
|
else if (type == SCM_UVEC_S32)
|
|
(((scm_t_int32*)base)[c_idx]) = scm_to_int32 (val);
|
|
#if SCM_HAVE_T_INT64
|
|
else if (type == SCM_UVEC_U64)
|
|
(((scm_t_uint64*)base)[c_idx]) = scm_to_uint64 (val);
|
|
else if (type == SCM_UVEC_S64)
|
|
(((scm_t_int64*)base)[c_idx]) = scm_to_int64 (val);
|
|
#endif
|
|
else if (type == SCM_UVEC_F32)
|
|
(((float*)base)[c_idx]) = scm_to_double (val);
|
|
else if (type == SCM_UVEC_F64)
|
|
(((double*)base)[c_idx]) = scm_to_double (val);
|
|
}
|
|
|
|
static SCM_C_INLINE SCM
|
|
make_uvec (int type, SCM len, SCM fill)
|
|
{
|
|
size_t c_len = scm_to_unsigned_integer (len, 0, SIZE_MAX / uvec_sizes[type]);
|
|
SCM uvec = alloc_uvec (type, c_len);
|
|
if (!SCM_UNBNDP (fill))
|
|
{
|
|
size_t idx;
|
|
void *base = SCM_UVEC_BASE (uvec);
|
|
for (idx = 0; idx < c_len; idx++)
|
|
uvec_fast_set_x (type, base, idx, fill);
|
|
}
|
|
return uvec;
|
|
}
|
|
|
|
static SCM_C_INLINE SCM
|
|
uvec_length (int type, SCM uvec)
|
|
{
|
|
uvec_assert (type, uvec);
|
|
return scm_from_size_t (SCM_UVEC_LENGTH (uvec));
|
|
}
|
|
|
|
static SCM_C_INLINE SCM
|
|
uvec_ref (int type, SCM uvec, SCM idx)
|
|
{
|
|
size_t c_idx;
|
|
SCM res;
|
|
|
|
uvec_assert (type, uvec);
|
|
c_idx = scm_to_unsigned_integer (idx, 0, SCM_UVEC_LENGTH (uvec)-1);
|
|
res = uvec_fast_ref (type, SCM_UVEC_BASE(uvec), c_idx);
|
|
scm_remember_upto_here_1 (uvec);
|
|
return res;
|
|
}
|
|
|
|
static SCM_C_INLINE SCM
|
|
uvec_set_x (int type, SCM uvec, SCM idx, SCM val)
|
|
{
|
|
size_t c_idx;
|
|
|
|
uvec_assert (type, uvec);
|
|
c_idx = scm_to_unsigned_integer (idx, 0, SCM_UVEC_LENGTH (uvec)-1);
|
|
uvec_fast_set_x (type, SCM_UVEC_BASE(uvec), c_idx, val);
|
|
scm_remember_upto_here_1 (uvec);
|
|
return SCM_UNSPECIFIED;
|
|
}
|
|
|
|
static SCM_C_INLINE SCM
|
|
uvec_to_list (int type, SCM uvec)
|
|
{
|
|
size_t c_idx;
|
|
void *base;
|
|
SCM res = SCM_EOL;
|
|
|
|
uvec_assert (type, uvec);
|
|
c_idx = SCM_UVEC_LENGTH (uvec);
|
|
base = SCM_UVEC_BASE (uvec);
|
|
while (c_idx-- > 0)
|
|
res = scm_cons (uvec_fast_ref (type, base, c_idx), res);
|
|
scm_remember_upto_here_1 (uvec);
|
|
return res;
|
|
}
|
|
|
|
static SCM_C_INLINE SCM
|
|
list_to_uvec (int type, SCM list)
|
|
{
|
|
SCM uvec;
|
|
void *base;
|
|
long idx;
|
|
long len = scm_ilength (list);
|
|
if (len < 0)
|
|
scm_wrong_type_arg_msg (NULL, 0, list, "proper list");
|
|
|
|
uvec = alloc_uvec (type, len);
|
|
base = SCM_UVEC_BASE (uvec);
|
|
idx = 0;
|
|
while (scm_is_pair (list) && idx < len)
|
|
{
|
|
uvec_fast_set_x (type, base, idx, SCM_CAR (list));
|
|
list = SCM_CDR (list);
|
|
idx++;
|
|
}
|
|
return uvec;
|
|
}
|
|
|
|
SCM
|
|
scm_i_read_homogenous_vector (SCM port, char pfx)
|
|
{
|
|
/* We have read '#f', '#u', or '#s'. Next must be a decimal integer
|
|
followed immediately by a list.
|
|
*/
|
|
|
|
int c;
|
|
char tok[80];
|
|
int n_digs;
|
|
SCM list;
|
|
|
|
n_digs = 0;
|
|
while ((c = scm_getc (port)) != EOF && '0' <= c && c <= '9' && n_digs < 80)
|
|
tok[n_digs++] = c;
|
|
|
|
if (c != EOF)
|
|
scm_ungetc (c, port);
|
|
|
|
if (n_digs == 0 && pfx == 'f')
|
|
return SCM_BOOL_F;
|
|
|
|
if (c != '(')
|
|
scm_i_input_error (NULL, port,
|
|
"#~a~a must be followed immediately by a '('",
|
|
scm_list_2 (SCM_MAKE_CHAR (pfx),
|
|
scm_from_locale_stringn (tok, n_digs)));
|
|
|
|
list = scm_read (port);
|
|
|
|
if (n_digs == 1 && strncmp (tok, "8", n_digs) == 0)
|
|
{
|
|
if (pfx == 'u')
|
|
return scm_list_to_u8vector (list);
|
|
else if (pfx == 's')
|
|
return scm_list_to_s8vector (list);
|
|
}
|
|
else if (n_digs == 2 && strncmp (tok, "16", n_digs) == 0)
|
|
{
|
|
if (pfx == 'u')
|
|
return scm_list_to_u16vector (list);
|
|
else if (pfx == 's')
|
|
return scm_list_to_s16vector (list);
|
|
}
|
|
else if (n_digs == 2 && strncmp (tok, "32", n_digs) == 0)
|
|
{
|
|
if (pfx == 'u')
|
|
return scm_list_to_u32vector (list);
|
|
else if (pfx == 's')
|
|
return scm_list_to_s32vector (list);
|
|
else if (pfx == 'f')
|
|
return scm_list_to_f32vector (list);
|
|
}
|
|
else if (n_digs == 2 && strncmp (tok, "64", n_digs) == 0)
|
|
{
|
|
if (pfx == 'u')
|
|
return scm_list_to_u64vector (list);
|
|
else if (pfx == 's')
|
|
return scm_list_to_s64vector (list);
|
|
else if (pfx == 'f')
|
|
return scm_list_to_f64vector (list);
|
|
}
|
|
|
|
scm_i_input_error (NULL, port,
|
|
"unrecognized homogenous vector prefix #~a~a",
|
|
scm_list_2 (SCM_MAKE_CHAR (pfx),
|
|
scm_from_locale_stringn (tok, n_digs)));
|
|
return SCM_BOOL_F;
|
|
}
|
|
|
|
SCM
|
|
scm_i_uniform_vector_prototype (SCM uvec)
|
|
{
|
|
switch (SCM_UVEC_TYPE (uvec))
|
|
{
|
|
case SCM_UVEC_U8:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_S8:
|
|
return SCM_MAKE_CHAR ('\0');
|
|
case SCM_UVEC_U16:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_S16:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_U32:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_S32:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_U64:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_S64:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_F32:
|
|
return SCM_BOOL_F;
|
|
case SCM_UVEC_F64:
|
|
return SCM_BOOL_F;
|
|
default:
|
|
return SCM_BOOL_F;
|
|
}
|
|
}
|
|
|
|
int
|
|
scm_is_uniform_vector (SCM obj)
|
|
{
|
|
return SCM_SMOB_PREDICATE (scm_tc16_uvec, obj);
|
|
}
|
|
|
|
size_t
|
|
scm_c_uniform_vector_length (SCM v)
|
|
{
|
|
if (scm_is_uniform_vector (v))
|
|
return SCM_UVEC_LENGTH (v);
|
|
else
|
|
scm_wrong_type_arg_msg (NULL, 0, v, "uniform vector");
|
|
}
|
|
|
|
size_t
|
|
scm_c_uniform_vector_size (SCM v)
|
|
{
|
|
if (scm_is_uniform_vector (v))
|
|
return SCM_UVEC_LENGTH (v) * uvec_sizes[SCM_UVEC_TYPE (v)];
|
|
else
|
|
scm_wrong_type_arg_msg (NULL, 0, v, "uniform vector");
|
|
}
|
|
|
|
SCM_DEFINE (scm_uniform_vector_p, "uniform-vector?", 1, 0, 0,
|
|
(SCM obj),
|
|
"Return @code{#t} if @var{obj} is a uniform vector.")
|
|
#define FUNC_NAME s_scm_uniform_vector_p
|
|
{
|
|
return scm_from_bool (scm_is_uniform_vector (obj));
|
|
}
|
|
#undef FUNC_NAME
|
|
|
|
SCM_DEFINE (scm_uniform_vector_ref, "uniform-vector-ref", 2, 0, 0,
|
|
(SCM v, SCM idx),
|
|
"Return the element at index @var{idx} of the\n"
|
|
"homogenous numeric vector @var{v}.")
|
|
#define FUNC_NAME s_scm_uniform_vector_ref
|
|
{
|
|
/* Support old argument convention.
|
|
*/
|
|
if (scm_is_pair (idx))
|
|
{
|
|
if (!scm_is_null (SCM_CDR (idx)))
|
|
scm_wrong_num_args (NULL);
|
|
idx = SCM_CAR (idx);
|
|
}
|
|
|
|
if (scm_is_uniform_vector (v))
|
|
return uvec_ref (SCM_UVEC_TYPE (v), v, idx);
|
|
else
|
|
scm_wrong_type_arg_msg (NULL, 0, v, "uniform vector");
|
|
}
|
|
#undef FUNC_NAME
|
|
|
|
SCM_DEFINE (scm_uniform_vector_set_x, "uniform-vector-set!", 3, 0, 0,
|
|
(SCM v, SCM idx, SCM val),
|
|
"Set the element at index @var{idx} of the\n"
|
|
"homogenous numeric vector @var{v} to @var{val}.")
|
|
#define FUNC_NAME s_scm_uniform_vector_set_x
|
|
{
|
|
/* Support old argument convention.
|
|
*/
|
|
if (scm_is_pair (idx))
|
|
{
|
|
if (!scm_is_null (SCM_CDR (idx)))
|
|
scm_wrong_num_args (NULL);
|
|
idx = SCM_CAR (idx);
|
|
}
|
|
|
|
if (scm_is_uniform_vector (v))
|
|
return uvec_set_x (SCM_UVEC_TYPE (v), v, idx, val);
|
|
else
|
|
scm_wrong_type_arg_msg (NULL, 0, v, "uniform vector");
|
|
}
|
|
#undef FUNC_NAME
|
|
|
|
SCM_DEFINE (scm_uniform_vector_to_list, "uniform-vector->list", 1, 0, 0,
|
|
(SCM uvec),
|
|
"Convert the homogeneous numeric vector @var{uvec} to a list.")
|
|
#define FUNC_NAME s_uniform_vector_to_list
|
|
{
|
|
if (scm_is_uniform_vector (uvec))
|
|
return uvec_to_list (SCM_UVEC_TYPE (uvec), uvec);
|
|
else
|
|
scm_wrong_type_arg_msg (NULL, 0, uvec, "uniform vector");
|
|
}
|
|
#undef FUNC_NAME
|
|
|
|
void *
|
|
scm_uniform_vector_elements (SCM uvec)
|
|
{
|
|
if (scm_is_uniform_vector (uvec))
|
|
return SCM_UVEC_BASE (uvec);
|
|
else
|
|
scm_wrong_type_arg_msg (NULL, 0, uvec, "uniform vector");
|
|
}
|
|
|
|
void
|
|
scm_uniform_vector_release (SCM uvec)
|
|
{
|
|
/* Nothing to do right now, but this function might come in handy
|
|
when uniform vectors need to be locked when giving away a pointer
|
|
to their elements.
|
|
*/
|
|
}
|
|
|
|
size_t
|
|
scm_uniform_vector_element_size (SCM uvec)
|
|
{
|
|
if (scm_is_uniform_vector (uvec))
|
|
return uvec_sizes[SCM_UVEC_TYPE (uvec)];
|
|
else
|
|
scm_wrong_type_arg_msg (NULL, 0, uvec, "uniform vector");
|
|
}
|
|
|
|
/* return the size of an element in a uniform array or 0 if type not
|
|
found. */
|
|
size_t
|
|
scm_uniform_element_size (SCM obj)
|
|
{
|
|
size_t result;
|
|
|
|
if (scm_is_uniform_vector (obj))
|
|
return scm_uniform_vector_element_size (obj);
|
|
|
|
switch (SCM_TYP7 (obj))
|
|
{
|
|
case scm_tc7_bvect:
|
|
case scm_tc7_uvect:
|
|
case scm_tc7_ivect:
|
|
result = sizeof (long);
|
|
break;
|
|
|
|
case scm_tc7_svect:
|
|
result = sizeof (short);
|
|
break;
|
|
|
|
#if SCM_SIZEOF_LONG_LONG != 0
|
|
case scm_tc7_llvect:
|
|
result = sizeof (long long);
|
|
break;
|
|
#endif
|
|
|
|
case scm_tc7_fvect:
|
|
result = sizeof (float);
|
|
break;
|
|
|
|
case scm_tc7_dvect:
|
|
result = sizeof (double);
|
|
break;
|
|
|
|
case scm_tc7_cvect:
|
|
result = 2 * sizeof (double);
|
|
break;
|
|
|
|
default:
|
|
result = 0;
|
|
}
|
|
return result;
|
|
}
|
|
|
|
SCM_DEFINE (scm_uniform_vector_length, "uniform-vector-length", 1, 0, 0,
|
|
(SCM v),
|
|
"Return the number of elements in @var{uve}.")
|
|
#define FUNC_NAME s_scm_uniform_vector_length
|
|
{
|
|
if (scm_is_uniform_vector (v))
|
|
return scm_from_size_t (SCM_UVEC_LENGTH (v));
|
|
|
|
SCM_ASRTGO (SCM_NIMP (v), badarg1);
|
|
switch SCM_TYP7 (v)
|
|
{
|
|
default:
|
|
badarg1:SCM_WRONG_TYPE_ARG (1, v);
|
|
case scm_tc7_vector:
|
|
case scm_tc7_wvect:
|
|
return scm_from_size_t (SCM_VECTOR_LENGTH (v));
|
|
case scm_tc7_string:
|
|
return scm_from_size_t (scm_i_string_length (v));
|
|
case scm_tc7_bvect:
|
|
return scm_from_size_t (SCM_BITVECTOR_LENGTH (v));
|
|
case scm_tc7_uvect:
|
|
case scm_tc7_ivect:
|
|
case scm_tc7_fvect:
|
|
case scm_tc7_dvect:
|
|
case scm_tc7_cvect:
|
|
case scm_tc7_svect:
|
|
#if SCM_SIZEOF_LONG_LONG != 0
|
|
case scm_tc7_llvect:
|
|
#endif
|
|
return scm_from_size_t (SCM_UVECTOR_LENGTH (v));
|
|
}
|
|
}
|
|
#undef FUNC_NAME
|
|
|
|
/* ================================================================ */
|
|
/* Exported procedures. */
|
|
/* ================================================================ */
|
|
|
|
#define TYPE SCM_UVEC_U8
|
|
#define TAG u8
|
|
#define CTYPE scm_t_uint8
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_S8
|
|
#define TAG s8
|
|
#define CTYPE scm_t_int8
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_U16
|
|
#define TAG u16
|
|
#define CTYPE scm_t_uint16
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_S16
|
|
#define TAG s16
|
|
#define CTYPE scm_t_int16
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_U32
|
|
#define TAG u32
|
|
#define CTYPE scm_t_uint32
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_S32
|
|
#define TAG s32
|
|
#define CTYPE scm_t_int32
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_U64
|
|
#define TAG u64
|
|
#define CTYPE scm_t_uint64
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_S64
|
|
#define TAG s64
|
|
#define CTYPE scm_t_int64
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_F32
|
|
#define TAG f32
|
|
#define CTYPE float
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
#define TYPE SCM_UVEC_F64
|
|
#define TAG f64
|
|
#define CTYPE double
|
|
#include "libguile/srfi-4.i.c"
|
|
|
|
|
|
/* Create the smob type for homogeneous numeric vectors and install
|
|
the primitives. */
|
|
void
|
|
scm_init_srfi_4 (void)
|
|
{
|
|
scm_tc16_uvec = scm_make_smob_type ("uvec", 0);
|
|
scm_set_smob_equalp (scm_tc16_uvec, uvec_equalp);
|
|
scm_set_smob_free (scm_tc16_uvec, uvec_free);
|
|
scm_set_smob_print (scm_tc16_uvec, uvec_print);
|
|
#include "libguile/srfi-4.x"
|
|
}
|
|
|
|
/* End of srfi-4.c. */
|