Viewing file:  PR100914.c (5.94 KB) -rw-r--r--
Select action/file-type:
 (+) |  (+) |  (+) | Code (+) | Session (+) |  (+) | SDB (+) |  (+) |  (+) |  (+) |  (+) |  (+) |
/* Test the fix for PR100914 */
#include <assert.h>
#include <complex.h>
#include <stdbool.h>
#include <stdio.h>
#include <math.h>
#include <ISO_Fortran_binding.h>
#define _CFI_type_mask 0xFF
#define _CFI_type_kind_shift 8
#define _CFI_decode_type(NAME) (signed char)((NAME) & CFI_type_mask)
#define _CFI_decode_kind(NAME) (signed char)(((NAME) >> CFI_type_kind_shift) & CFI_type_mask)
#define _CFI_encode_type(TYPE, KIND) (int16_t)\
((((KIND) & CFI_type_mask) << CFI_type_kind_shift)\
| ((TYPE) & CFI_type_mask))
#undef CMPLXF
#define CMPLXF(x, y) ((float complex)((float)(x) + I * (float)(y)))
#undef CMPLX
#define CMPLX(x, y) ((double complex)((double)(x) + (double complex)I * (double)(y)))
#undef CMPLXL
#define CMPLXL(x, y) ((long double complex)((long double)(x) + (long double complex)I * (long double)(y)))
#undef CMPLX
#define CMPLX(x, y) ((_Float128 _Complex )((double)(x) + (double complex)I * (double)(y)))
#define N 11
#define M 7
typedef float _Complex c_float_complex;
typedef double _Complex c_double_complex;
typedef long double _Complex c_long_double_complex;
typedef _Float128 _Complex c_float128_complex;
bool c_vrfy_c_float_complex (const CFI_cdesc_t *restrict);
bool c_vrfy_c_double_complex (const CFI_cdesc_t *restrict);
bool c_vrfy_c_long_double_complex (const CFI_cdesc_t *restrict);
bool c_vrfy_c_float128_complex (const CFI_cdesc_t *restrict);
bool c_vrfy_complex (const CFI_cdesc_t *restrict);
bool c_vrfy_desc (const CFI_cdesc_t*restrict, const CFI_type_t, const signed char, const size_t, const size_t);
void check_tk (const CFI_cdesc_t*restrict, const CFI_type_t, const signed char, const size_t, const size_t);
bool
c_vrfy_c_float_complex (const CFI_cdesc_t *restrict auxp)
{
CFI_index_t i, lb, ub, ex;
size_t sz;
c_float_complex *ip = NULL;
assert (auxp);
assert (auxp->base_addr);
assert (auxp->elem_len>0);
lb = auxp->dim[0].lower_bound;
ex = auxp->dim[0].extent;
assert (ex==11);
sz = (size_t)auxp->elem_len / sizeof (c_float_complex);
assert (sz==1);
ub = ex + lb - 1;
ip = (c_float_complex*)auxp->base_addr;
for (i=0; i<ex; i++, ip+=sz)
if ((cabsf (*ip-(c_float_complex)(CMPLXF((i+1), (2*(i+1)))))>(float)0.0))
return false;
for (i=lb; i<ub+1; i++)
{
ip = (c_float_complex*)CFI_address(auxp, &i);
if ((cabsf (*ip-(c_float_complex)(CMPLXF((i-lb+1), (2*(i-lb+1)))))>(float)0.0))
return false;
}
return true;
}
bool
c_vrfy_c_double_complex (const CFI_cdesc_t *restrict auxp)
{
CFI_index_t i, lb, ub, ex;
size_t sz;
c_double_complex *ip = NULL;
assert (auxp);
assert (auxp->base_addr);
assert (auxp->elem_len>0);
lb = auxp->dim[0].lower_bound;
ex = auxp->dim[0].extent;
assert (ex==11);
sz = (size_t)auxp->elem_len / sizeof (c_double_complex);
assert (sz==1);
ub = ex + lb - 1;
ip = (c_double_complex*)auxp->base_addr;
for (i=0; i<ex; i++, ip+=sz)
if ((cabs (*ip-(c_double_complex)(CMPLX((i+1), (2*(i+1)))))>(double)0.0))
return false;
for (i=lb; i<ub+1; i++)
{
ip = (c_double_complex*)CFI_address(auxp, &i);
if ((cabs (*ip-(c_double_complex)(CMPLX((i-lb+1), (2*(i-lb+1)))))>(double)0.0))
return false;
}
return true;
}
bool
c_vrfy_c_long_double_complex (const CFI_cdesc_t *restrict auxp)
{
CFI_index_t i, lb, ub, ex;
size_t sz;
c_long_double_complex *ip = NULL;
assert (auxp);
assert (auxp->base_addr);
assert (auxp->elem_len>0);
lb = auxp->dim[0].lower_bound;
ex = auxp->dim[0].extent;
assert (ex==11);
sz = (size_t)auxp->elem_len / sizeof (c_long_double_complex);
assert (sz==1);
ub = ex + lb - 1;
ip = (c_long_double_complex*)auxp->base_addr;
for (i=0; i<ex; i++, ip+=sz)
if ((cabsl (*ip-(c_long_double_complex)(CMPLXL((i+1), (2*(i+1)))))>(long double)0.0))
return false;
for (i=lb; i<ub+1; i++)
{
ip = (c_long_double_complex*)CFI_address(auxp, &i);
if ((cabsl (*ip-(c_long_double_complex)(CMPLXL((i-lb+1), (2*(i-lb+1)))))>(long double)0.0))
return false;
}
return true;
}
bool
c_vrfy_c_float128_complex (const CFI_cdesc_t *restrict auxp)
{
CFI_index_t i, lb, ub, ex;
size_t sz;
c_float128_complex *ip = NULL;
assert (auxp);
assert (auxp->base_addr);
assert (auxp->elem_len>0);
lb = auxp->dim[0].lower_bound;
ex = auxp->dim[0].extent;
assert (ex==11);
sz = (size_t)auxp->elem_len / sizeof (c_float128_complex);
assert (sz==1);
ub = ex + lb - 1;
ip = (c_float128_complex*)auxp->base_addr;
for (i=0; i<ex; i++, ip+=sz)
if ((cabs ((double complex)(*ip-(c_float128_complex)(CMPLX((i+1), (2*(i+1))))))>(double)0.0))
return false;
for (i=lb; i<ub+1; i++)
{
ip = (c_float128_complex*)CFI_address(auxp, &i);
if ((cabs ((double complex)(*ip-(c_float128_complex)(CMPLX((i-lb+1), (2*(i-lb+1))))))>(double)0.0))
return false;
}
return true;
}
bool
c_vrfy_complex (const CFI_cdesc_t *restrict auxp)
{
signed char type, kind;
assert (auxp);
type = _CFI_decode_type(auxp->type);
kind = _CFI_decode_kind(auxp->type);
assert (type == CFI_type_Complex);
switch (kind)
{
case 4:
return c_vrfy_c_float_complex (auxp);
break;
case 8:
return c_vrfy_c_double_complex (auxp);
break;
case 10:
return c_vrfy_c_long_double_complex (auxp);
break;
case 16:
return c_vrfy_c_float128_complex (auxp);
break;
default:
assert (false);
}
return true;
}
void
check_tk (const CFI_cdesc_t *restrict auxp, const CFI_type_t type, const signed char kind, const size_t elem_len, const size_t nelem)
{
signed char ityp, iknd;
assert (auxp);
assert (auxp->elem_len==elem_len*nelem);
assert (auxp->rank==1);
assert (auxp->dim[0].sm>0);
assert ((size_t)auxp->dim[0].sm==elem_len*nelem);
/* */
assert (auxp->type==type);
ityp = _CFI_decode_type(auxp->type);
assert (ityp == CFI_type_Complex);
iknd = _CFI_decode_kind(auxp->type);
assert (_CFI_decode_type(type)==ityp);
assert (kind==iknd);
assert (c_vrfy_complex (auxp));
return;
}
// Local Variables:
// mode: C
// End:
|
