1 `/* Implementation of the EOSHIFT intrinsic
2 Copyright (C) 2002-2018 Free Software Foundation, Inc.
3 Contributed by Paul Brook <paul@nowt.org>
5 This file is part of the GNU Fortran runtime library (libgfortran).
7 Libgfortran is free software; you can redistribute it and/or
8 modify it under the terms of the GNU General Public
9 License as published by the Free Software Foundation; either
10 version 3 of the License, or (at your option) any later version.
12 Libgfortran is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 GNU General Public License for more details.
17 Under Section 7 of GPL version 3, you are granted additional
18 permissions described in the GCC Runtime Library Exception, version
19 3.1, as published by the Free Software Foundation.
21 You should have received a copy of the GNU General Public License and
22 a copy of the GCC Runtime Library Exception along with this program;
23 see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
24 <http://www.gnu.org/licenses/>. */
26 #include "libgfortran.h"
31 `#if defined (HAVE_'atype_name`)
34 eoshift1 (gfc_array_char * const restrict ret,
35 const gfc_array_char * const restrict array,
36 const 'atype` * const restrict h,
37 const char * const restrict pbound,
38 const 'atype_name` * const restrict pwhich,
39 const char * filler, index_type filler_len)
41 /* r.* indicates the return array. */
42 index_type rstride[GFC_MAX_DIMENSIONS];
47 /* s.* indicates the source array. */
48 index_type sstride[GFC_MAX_DIMENSIONS];
53 /* h.* indicates the shift array. */
54 index_type hstride[GFC_MAX_DIMENSIONS];
56 const 'atype_name` *hptr;
58 index_type count[GFC_MAX_DIMENSIONS];
59 index_type extent[GFC_MAX_DIMENSIONS];
69 /* The compiler cannot figure out that these are set, initialize
70 them to avoid warnings. */
75 size = GFC_DESCRIPTOR_SIZE(array);
85 arraysize = size0 ((array_t *) array);
86 if (ret->base_addr == NULL)
91 ret->dtype = array->dtype;
92 for (i = 0; i < GFC_DESCRIPTOR_RANK (array); i++)
96 ub = GFC_DESCRIPTOR_EXTENT(array,i) - 1;
101 str = GFC_DESCRIPTOR_EXTENT(ret,i-1)
102 * GFC_DESCRIPTOR_STRIDE(ret,i-1);
104 GFC_DIMENSION_SET(ret->dim[i], 0, ub, str);
107 /* xmallocarray allocates a single byte for zero size. */
108 ret->base_addr = xmallocarray (arraysize, size);
111 else if (unlikely (compile_options.bounds_check))
113 bounds_equal_extents ((array_t *) ret, (array_t *) array,
114 "return value", "EOSHIFT");
117 if (unlikely (compile_options.bounds_check))
119 bounds_reduced_extents ((array_t *) h, (array_t *) array, which,
120 "SHIFT argument", "EOSHIFT");
127 for (dim = 0; dim < GFC_DESCRIPTOR_RANK (array); dim++)
131 roffset = GFC_DESCRIPTOR_STRIDE_BYTES(ret,dim);
134 soffset = GFC_DESCRIPTOR_STRIDE_BYTES(array,dim);
137 len = GFC_DESCRIPTOR_EXTENT(array,dim);
142 extent[n] = GFC_DESCRIPTOR_EXTENT(array,dim);
143 rstride[n] = GFC_DESCRIPTOR_STRIDE_BYTES(ret,dim);
144 sstride[n] = GFC_DESCRIPTOR_STRIDE_BYTES(array,dim);
146 hstride[n] = GFC_DESCRIPTOR_STRIDE(h,n);
157 dim = GFC_DESCRIPTOR_RANK (array);
158 rstride0 = rstride[0];
159 sstride0 = sstride[0];
160 hstride0 = hstride[0];
161 rptr = ret->base_addr;
162 sptr = array->base_addr;
167 /* Do the shift for this dimension. */
169 if (( sh >= 0 ? sh : -sh ) > len)
175 delta = (sh >= 0) ? sh: -sh;
179 src = &sptr[delta * soffset];
185 dest = &rptr[delta * roffset];
188 /* If the elements are contiguous, perform a single block move. */
189 if (soffset == size && roffset == size)
191 size_t chunk = size * (len - delta);
192 memcpy (dest, src, chunk);
197 for (n = 0; n < len - delta; n++)
199 memcpy (dest, src, size);
211 memcpy (dest, pbound, size);
220 memset (dest, filler[0], size);
222 for (i = 0; i < size; i += filler_len)
223 memcpy (&dest[i], filler, filler_len);
228 /* Advance to the next section. */
234 while (count[n] == extent[n])
236 /* When we get to the end of a dimension, reset it and increment
237 the next dimension. */
239 /* We could precalculate these products, but this is a less
240 frequently used path so probably not worth it. */
241 rptr -= rstride[n] * extent[n];
242 sptr -= sstride[n] * extent[n];
243 hptr -= hstride[n] * extent[n];
247 /* Break out of the loop. */
262 void eoshift1_'atype_kind` (gfc_array_char * const restrict,
263 const gfc_array_char * const restrict,
264 const 'atype` * const restrict, const char * const restrict,
265 const 'atype_name` * const restrict);
266 export_proto(eoshift1_'atype_kind`);
269 eoshift1_'atype_kind` (gfc_array_char * const restrict ret,
270 const gfc_array_char * const restrict array,
271 const 'atype` * const restrict h,
272 const char * const restrict pbound,
273 const 'atype_name` * const restrict pwhich)
275 eoshift1 (ret, array, h, pbound, pwhich, "\0", 1);
279 void eoshift1_'atype_kind`_char (gfc_array_char * const restrict,
281 const gfc_array_char * const restrict,
282 const 'atype` * const restrict,
283 const char * const restrict,
284 const 'atype_name` * const restrict,
285 GFC_INTEGER_4, GFC_INTEGER_4);
286 export_proto(eoshift1_'atype_kind`_char);
289 eoshift1_'atype_kind`_char (gfc_array_char * const restrict ret,
290 GFC_INTEGER_4 ret_length __attribute__((unused)),
291 const gfc_array_char * const restrict array,
292 const 'atype` * const restrict h,
293 const char * const restrict pbound,
294 const 'atype_name` * const restrict pwhich,
295 GFC_INTEGER_4 array_length __attribute__((unused)),
296 GFC_INTEGER_4 bound_length __attribute__((unused)))
298 eoshift1 (ret, array, h, pbound, pwhich, " ", 1);
302 void eoshift1_'atype_kind`_char4 (gfc_array_char * const restrict,
304 const gfc_array_char * const restrict,
305 const 'atype` * const restrict,
306 const char * const restrict,
307 const 'atype_name` * const restrict,
308 GFC_INTEGER_4, GFC_INTEGER_4);
309 export_proto(eoshift1_'atype_kind`_char4);
312 eoshift1_'atype_kind`_char4 (gfc_array_char * const restrict ret,
313 GFC_INTEGER_4 ret_length __attribute__((unused)),
314 const gfc_array_char * const restrict array,
315 const 'atype` * const restrict h,
316 const char * const restrict pbound,
317 const 'atype_name` * const restrict pwhich,
318 GFC_INTEGER_4 array_length __attribute__((unused)),
319 GFC_INTEGER_4 bound_length __attribute__((unused)))
321 static const gfc_char4_t space = (unsigned char) ''` ''`;
322 eoshift1 (ret, array, h, pbound, pwhich,
323 (const char *) &space, sizeof (gfc_char4_t));