void index_boolean_array(const boolean_array_t* source,
const index_spec_t* source_spec,
boolean_array_t* dest)
{
_index_t* idx_vec1;
_index_t* idx_vec2;
_index_t* idx_size;
int j;
int i;
assert(base_array_ok(source));
assert(base_array_ok(dest));
assert(index_spec_ok(source_spec));
assert(index_spec_fit_base_array(source_spec,source));
for(i = 0, j = 0; i < source->ndims; ++i) {
if((source_spec->index_type[i] == 'W')
||
(source_spec->index_type[i] == 'A')) {
++j;
}
}
assert(j == dest->ndims);
idx_vec1 = size_alloc(source->ndims); /*indices in the source array*/
idx_vec2 = size_alloc(dest->ndims); /* indices in the destination array*/
idx_size = size_alloc(source_spec->ndims);
for(i = 0; i < source->ndims; ++i) {
idx_vec1[i] = 0;
}
for(i = 0; i < source_spec->ndims; ++i) {
if(source_spec->index[i]) {
idx_size[i] = imax(source_spec->dim_size[i],1);
} else {
idx_size[i] = source->dim_size[i];
}
}
do {
for(i = 0, j = 0; i < source->ndims; ++i) {
if((source_spec->index_type[i] == 'W')
||
(source_spec->index_type[i] == 'A')) {
idx_vec2[j] = idx_vec1[i];
j++;
}
}
boolean_set(dest, calc_base_index(dest->ndims, idx_vec2, dest),
boolean_get(*source,
calc_base_index_spec(source->ndims, idx_vec1,
source, source_spec)));
} while(0 == next_index(source->ndims, idx_vec1, idx_size));
}
void indexed_assign_string_array(const string_array_t * source,
string_array_t* dest,
const index_spec_t* dest_spec)
{
_index_t* idx_vec1;
_index_t* idx_vec2;
_index_t* idx_size;
int i,j;
assert(base_array_ok(source));
assert(base_array_ok(dest));
assert(index_spec_ok(dest_spec));
assert(index_spec_fit_base_array(dest_spec, dest));
for(i = 0,j = 0; i < dest_spec->ndims; ++i) {
if(dest_spec->dim_size[i] != 0) {
++j;
}
}
assert(j == source->ndims);
idx_vec1 = size_alloc(dest->ndims);
idx_vec2 = size_alloc(source->ndims);
idx_size = size_alloc(dest_spec->ndims);
for(i = 0; i < dest_spec->ndims; ++i) {
idx_vec1[i] = 0;
if(dest_spec->index[i] != NULL) {
idx_size[i] = imax(dest_spec->dim_size[i],1);
} else {
idx_size[i] = dest->dim_size[i];
}
}
do {
for(i = 0, j = 0; i < dest_spec->ndims; ++i) {
if(dest_spec->dim_size[i] != 0) {
idx_vec2[j] = idx_vec1[i];
++j;
}
}
string_set(dest, calc_base_index_spec(dest->ndims, idx_vec1,
dest, dest_spec),
string_get(*source, calc_base_index(source->ndims,
idx_vec2, source)));
} while(0 == next_index(dest_spec->ndims, idx_vec1, idx_size));
}
void indexed_assign_boolean_array(const boolean_array_t source, boolean_array_t* dest,
const index_spec_t* dest_spec)
{
_index_t *idx_vec1, *idx_size;
int j;
indexed_assign_base_array_size_alloc(&source, dest, dest_spec, &idx_vec1, &idx_size);
j = 0;
do {
boolean_set(dest,
calc_base_index_spec(dest->ndims, idx_vec1, dest, dest_spec),
boolean_get(source, j));
j++;
} while(0 == next_index(dest_spec->ndims, idx_vec1, idx_size));
omc_assert_macro(j == base_array_nr_of_elements(source));
}
void indexed_assign_real_array(const real_array_t source, real_array_t* dest,
const index_spec_t* dest_spec)
{
_index_t* idx_vec1;
_index_t* idx_size;
int i,j;
omc_assert_macro(base_array_ok(&source));
omc_assert_macro(base_array_ok(dest));
omc_assert_macro(index_spec_ok(dest_spec));
omc_assert_macro(index_spec_fit_base_array(dest_spec, dest));
for(i = 0,j = 0; i < dest_spec->ndims; ++i) {
if(dest_spec->dim_size[i] != 0) {
++j;
}
}
omc_assert_macro(j == source.ndims);
idx_vec1 = size_alloc(dest->ndims);
idx_size = size_alloc(dest_spec->ndims);
for(i = 0; i < dest_spec->ndims; ++i) {
idx_vec1[i] = 0;
if(dest_spec->index[i] != NULL) { /* is 'S' or 'A' */
idx_size[i] = imax(dest_spec->dim_size[i],1);
} else { /* is 'W' */
idx_size[i] = dest->dim_size[i];
}
}
j = 0;
do {
real_set(dest,
calc_base_index_spec(dest->ndims, idx_vec1, dest, dest_spec),
real_get(source, j));
j++;
} while(0 == next_index(dest_spec->ndims, idx_vec1, idx_size));
omc_assert_macro(j == base_array_nr_of_elements(source));
}
void index_real_array(const real_array_t * source,
const index_spec_t* source_spec,
real_array_t* dest)
{
_index_t* idx_vec1;
_index_t* idx_size;
int j;
int i;
state mem_state;
assert(base_array_ok(source));
assert(base_array_ok(dest));
assert(index_spec_ok(source_spec));
assert(index_spec_fit_base_array(source_spec,source));
/*for (i = 0, j = 0; i < source->ndims; ++i)
{
printf("source_spec->index_type[%d] = %c\n", i, source_spec->index_type[i]);
if ((source_spec->index_type[i] == 'W')
||
(source_spec->index_type[i] == 'A'))
++j;
}
assert(j == dest->ndims);*/
for (i = 0,j = 0; i < source_spec->ndims; ++i) {
if (source_spec->dim_size[i] != 0) {
++j;
}
}
assert(j == dest->ndims);
mem_state = get_memory_state();
idx_vec1 = size_alloc(0,source->ndims); /*indices in the source array*/
/* idx_vec2 = size_alloc(0,dest->ndims); / * indices in the destination array* / */
idx_size = size_alloc(0,source_spec->ndims);
for (i = 0; i < source->ndims; ++i) {
idx_vec1[i] = 0;
}
for (i = 0; i < source_spec->ndims; ++i) {
if (source_spec->index[i] != NULL) { /* is 'S' or 'A' */
idx_size[i] = imax(source_spec->dim_size[i],1); /* the imax() is not needed, because there is (idx[d] >= size[d]) in the next_index(), but ... */
} else { /* is 'W' */
idx_size[i] = source->dim_size[i];
}
}
j = 0;
do {
/*
for (i = 0, j = 0; i < source->ndims; ++i) {
if ((source_spec->index_type[i] == 'W')
||
(source_spec->index_type[i] == 'A')) {
idx_vec2[j] = idx_vec1[i];
j++;
}
}
*/
real_set(dest, j, /* calc_base_index(dest->ndims, idx_vec2, dest), */
real_get(source,
calc_base_index_spec(source->ndims, idx_vec1,
source, source_spec)));
j++;
} while (0 == next_index(source->ndims, idx_vec1, idx_size));
assert(j == base_array_nr_of_elements(dest));
restore_memory_state(mem_state);
}
void indexed_assign_real_array(const real_array_t * source,
real_array_t* dest,
const index_spec_t* dest_spec)
{
_index_t* idx_vec1;
_index_t* idx_size;
int i,j;
state mem_state;
assert(base_array_ok(source));
assert(base_array_ok(dest));
assert(index_spec_ok(dest_spec));
assert(index_spec_fit_base_array(dest_spec, dest));
for (i = 0,j = 0; i < dest_spec->ndims; ++i) {
if (dest_spec->dim_size[i] != 0) {
++j;
}
}
assert(j == source->ndims);
mem_state = get_memory_state();
idx_vec1 = size_alloc(0,dest->ndims);
/* idx_vec2 = size_alloc(0,source->ndims); */
idx_size = size_alloc(0,dest_spec->ndims);
for (i = 0; i < dest_spec->ndims; ++i) {
idx_vec1[i] = 0;
if (dest_spec->index[i] != NULL) { /* is 'S' or 'A' */
idx_size[i] = imax(dest_spec->dim_size[i],1);
} else { /* is 'W' */
idx_size[i] = dest->dim_size[i];
}
}
j = 0;
do {
/*for (i = 0, j = 0; i < dest_spec->ndims; ++i) {
if (dest_spec->dim_size[i] != 0) {
idx_vec2[j] = idx_vec1[i];
++j;
}
}*/
/*for (i = 0, j = 0; i < dest_spec->ndims; ++i) {
if ((dest_spec->index_type[i] == 'W')
||
(dest_spec->index_type[i] == 'A')) {
idx_vec2[j] = idx_vec1[i];
++j;
}
}*/
real_set(dest,
calc_base_index_spec(dest->ndims, idx_vec1, dest, dest_spec),
real_get(source, j));
j++;
/* calc_base_index(source->ndims, idx_vec1, source)));*/
/* calc_base_index(source->ndims, idx_vec2, source)));*/
} while (0 == next_index(dest_spec->ndims, idx_vec1, idx_size));
assert(j == base_array_nr_of_elements(source));
restore_memory_state(mem_state);
}
static inline modelica_real *calc_real_index_spec(int ndims, const _index_t *idx_vec,
const real_array_t *arr,
const index_spec_t *spec)
{
return real_ptrget(arr, calc_base_index_spec(ndims, idx_vec, arr, spec));
}
static inline modelica_boolean *calc_boolean_index_spec(int ndims, const _index_t* idx_vec,
const boolean_array_t *arr,
const index_spec_t *spec)
{
return boolean_ptrget(arr, calc_base_index_spec(ndims, idx_vec, arr, spec));
}
static inline modelica_string *calc_string_index_spec(int ndims, const _index_t *idx_vec,
const string_array_t *arr,
const index_spec_t *spec)
{
return string_ptrget(arr, calc_base_index_spec(ndims, idx_vec, arr, spec));
}
