void mul_alloc_real_matrix_product_smart(const real_array_t * a, const real_array_t * b, real_array_t* dest)
{
if ((a->ndims == 1) && (b->ndims == 2)) {
simple_alloc_1d_real_array(dest,b->dim_size[1]);
mul_real_vector_matrix(a,b,dest);
} else if ((a->ndims == 2) && (b->ndims == 1)) {
simple_alloc_1d_real_array(dest,a->dim_size[0]);
mul_real_matrix_vector(a,b,dest);
} else if ((a->ndims == 2) && (b->ndims == 2)) {
simple_alloc_2d_real_array(dest,a->dim_size[0],b->dim_size[1]);
mul_real_matrix_product(a,b,dest);
} else {
printf("Invalid size of matrix\n");
}
}
void range_alloc_real_array(modelica_real start, modelica_real stop, modelica_real inc, real_array_t* dest)
{
int n;
n = floor((stop-start)/inc)+1;
simple_alloc_1d_real_array(dest,n);
range_real_array(start,stop,inc,dest);
}
real_array_t mul_alloc_real_matrix_product_smart(const real_array_t a, const real_array_t b)
{
real_array_t dest;
if((a.ndims == 1) && (b.ndims == 2)) {
simple_alloc_1d_real_array(&dest,b.dim_size[1]);
mul_real_vector_matrix(&a,&b,&dest);
} else if((a.ndims == 2) && (b.ndims == 1)) {
simple_alloc_1d_real_array(&dest,a.dim_size[0]);
mul_real_matrix_vector(&a,&b,&dest);
} else if((a.ndims == 2) && (b.ndims == 2)) {
simple_alloc_2d_real_array(&dest,a.dim_size[0],b.dim_size[1]);
mul_real_matrix_product(&a,&b,&dest);
} else {
omc_assert_macro(0 == "Invalid size of matrix");
}
return dest;
}
void array_alloc_scalar_real_array(real_array_t* dest, int n, modelica_real first,...)
{
int i;
va_list ap;
simple_alloc_1d_real_array(dest,n);
va_start(ap,first);
put_real_element(first,0,dest);
for (i = 1; i < n; ++i) {
put_real_element(va_arg(ap,modelica_real),i,dest);
}
va_end(ap);
}
/* Creates a real array from a range with a start, stop and step value.
* Ex: 1.0:2.0:6.0 => {1.0,3.0,5.0} */
void create_real_array_from_range(real_array_t *dest, modelica_real start, modelica_real step, modelica_real stop)
{
size_t elements;
size_t i;
modelica_real (*comp_func)(modelica_real, modelica_real);
omc_assert_macro(step != 0);
comp_func = (step > 0) ? &real_le : &real_ge;
elements = comp_func(start, stop) ? (((stop - start) / step) + 1) : 0;
/* fprintf(stderr, "start %g step %g stop %g elements %d\n", start, step, stop, elements); */
simple_alloc_1d_real_array(dest, elements);
for(i = 0; i < elements; start += step, ++i) {
real_set(dest, i, start);
}
}
