00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036 #include <cray/portdefs.h>
00037 #include <cray/dopevec.h>
00038 #include <liberrno.h>
00039 #include <float.h>
00040 #include <inttypes.h>
00041 #include <stdlib.h>
00042 #include <string.h>
00043
00044
00045 #define MAX_NARY_DIMS MAXDIM
00046 #define BIGDEFAULTSZ 64
00047
00048 #define ERROR _lerror
00049 #define FALSE 0
00050 #define TRUE 1
00051
00052
00053
00054
00055 #define GET_RANK_FROM_DESC(array) (array->n_dim)
00056 #define GET_EXTENT_FROM_DESC(array,dim) (array->dimension[dim].extent)
00057 #define GET_ASSOCIATED_FROM_DESC(array) (array->assoc)
00058 #define GET_ALEN_FROM_DESC(array) (array->base_addr.a.el_len)
00059 #define GET_DV_ASCII_FROM_DESC(array) (array->type_lens.type == DVTYPE_ASCII)
00060 #define GET_DV_LOGICAL_FROM_DESC(array) (array->type_lens.type == DVTYPE_LOGICAL)
00061
00062 #define GET_ADDRESS_FROM_DESC(array) \
00063 ((array->type_lens.type == DVTYPE_ASCII) ? \
00064 ((char *) _fcdtocp (array->base_addr.charptr)) : \
00065 ((char *) array->base_addr.a.ptr))
00066
00067 #define GET_BYTEALIGNED_FROM_DESC(array) \
00068 ((array->type_lens.type == DVTYPE_ASCII) || \
00069 (array->type_lens.type == DVTYPE_DERIVEDBYTE))
00070
00071 #define GET_ELEMENT_SZ_FROM_DESC(array) \
00072 ((array->type_lens.type == DVTYPE_ASCII)? \
00073 (_fcdlen(array->base_addr.charptr)) : \
00074 ((array->base_addr.a.el_len)>>3))
00075
00076 static size_t GET_STRIDE_FROM_DESC(DopeVectorType * array, int32_t dim)
00077 {
00078 size_t rval;
00079 rval = array->dimension[dim].stride_mult;
00080 switch (array->type_lens.type) {
00081 case DVTYPE_ASCII:
00082 case DVTYPE_DERIVEDBYTE:
00083 return (rval);
00084
00085 case DVTYPE_DERIVEDWORD:
00086 return (rval << 2 );
00087
00088 default:
00089 switch (array->type_lens.int_len) {
00090 case 8:
00091 return (rval);
00092
00093 case 16:
00094 return (rval << 1);
00095
00096 default:
00097 return ( rval <<2 );
00098 }
00099 }
00100 }
00101
00102 #define MK_STRIDE(byte_aligned,typ_sz) \
00103 ((byte_aligned || typ_sz == 1) ? \
00104 (str_sz = typ_sz) : \
00105 ((typ_sz == 2) ? \
00106 (typ_sz/(BYTES_PER_WORD/2)) :\
00107 (typ_sz/BYTES_PER_WORD)))
00108
00109
00110 #define SET_RANK_IN_DESC(array,r) (array->n_dim = r)
00111 #define SET_ADDRESS_IN_DESC(array,d) (array->base_addr.a.ptr = d)
00112 #define SET_CHARPTR_IN_DESC(array,cr,sz) (array->base_addr.charptr = _cptofcd (cr, sz))
00113 #define SET_LBOUND_IN_DESC(array,dim,b) (array->dimension[dim].low_bound = b)
00114 #define SET_EXTENT_IN_DESC(array,dim,e) (array->dimension[dim].extent = e)
00115 #define SET_STRMULT_IN_DESC(array,dim,e)(array->dimension[dim].stride_mult = e)
00116 #define SET_ASSOCIATED_IN_DESC(array) (array->assoc = TRUE)
00117 #define SET_CONTIG_IN_DESC(array) (array->a_contig = TRUE)
00118 #define SET_ORIG_BS_IN_DESC(array,a) (array->orig_base = a)
00119 #define SET_ORIG_SZ_IN_DESC(array,s) (array->orig_size = s)
00120 #define SET_ALEN_IN_DESC(array,l) (array->base_addr.a.el_len = l)
00121
00122
00123
00124
00125 #if defined(_LITTLE_ENDIAN)
00126 #define OFFSET_TO_TF_BYTE(by) 0
00127 #else
00128 #define OFFSET_TO_TF_BYTE(by) (by-1)
00129 #endif
00130
00131
00132
00133
00134 #define i8 int64_t
00135 #define i4 int32_t
00136 #define i2 int16_t
00137 #define i1 int8_t
00138
00139 #define ui8 uint64_t
00140 #define ui4 uint32_t
00141 #define ui2 uint16_t
00142 #define ui1 uint8_t
00143
00144 #define l8 uint64_t
00145 #define l4 uint32_t
00146 #define l2 uint16_t
00147 #define l1 uint8_t
00148
00149 #define r16 long double
00150 #define r8 double
00151 #define r4 float
00152
00153 typedef struct { float r, i; } complex;
00154 typedef struct { double r, i; } dcomplex;
00155 typedef struct { long double r, i ; } qcomplex;
00156
00157 #define c32 qcomplex
00158 #define c16 dcomplex
00159 #define c8 complex
00160
00161 #define ALIGNED_r16(address) (((size_t)address & 0x7) == 0)
00162 #define ALIGNED_r8(address) (((size_t)address & 0x7) == 0)
00163 #define ALIGNED_i8(address) (((size_t)address & 0x7) == 0)
00164 #define ALIGNED_i4(address) (((size_t)address & 0x3) == 0)
00165 #define ALIGNED_r4(address) (((size_t)address & 0x3) == 0)
00166 #define ALIGNED_i2(address) (((size_t)address & 0x1) == 0)
00167 #define ALIGNED_i1(address) (TRUE)
00168
00169
00170