1 /* Implementation of the EOSHIFT intrinsic
2 Copyright (C) 2002-2017 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"
30 #if defined (HAVE_GFC_INTEGER_16)
33 eoshift3 (gfc_array_char
* const restrict ret
,
34 const gfc_array_char
* const restrict array
,
35 const gfc_array_i16
* const restrict h
,
36 const gfc_array_char
* const restrict bound
,
37 const GFC_INTEGER_16
* const restrict pwhich
,
38 const char * filler
, index_type filler_len
)
40 /* r.* indicates the return array. */
41 index_type rstride
[GFC_MAX_DIMENSIONS
];
46 /* s.* indicates the source array. */
47 index_type sstride
[GFC_MAX_DIMENSIONS
];
52 /* h.* indicates the shift array. */
53 index_type hstride
[GFC_MAX_DIMENSIONS
];
55 const GFC_INTEGER_16
*hptr
;
56 /* b.* indicates the bound array. */
57 index_type bstride
[GFC_MAX_DIMENSIONS
];
61 index_type count
[GFC_MAX_DIMENSIONS
];
62 index_type extent
[GFC_MAX_DIMENSIONS
];
72 /* The compiler cannot figure out that these are set, initialize
73 them to avoid warnings. */
78 arraysize
= size0 ((array_t
*) array
);
79 size
= GFC_DESCRIPTOR_SIZE(array
);
86 if (ret
->base_addr
== NULL
)
90 ret
->base_addr
= xmallocarray (arraysize
, size
);
92 ret
->dtype
= array
->dtype
;
93 for (i
= 0; i
< GFC_DESCRIPTOR_RANK (array
); i
++)
97 ub
= GFC_DESCRIPTOR_EXTENT(array
,i
) - 1;
102 str
= GFC_DESCRIPTOR_EXTENT(ret
,i
-1)
103 * GFC_DESCRIPTOR_STRIDE(ret
,i
-1);
105 GFC_DIMENSION_SET(ret
->dim
[i
], 0, ub
, str
);
108 /* xmallocarray allocates a single byte for zero size. */
109 ret
->base_addr
= xmallocarray (arraysize
, size
);
112 else if (unlikely (compile_options
.bounds_check
))
114 bounds_equal_extents ((array_t
*) ret
, (array_t
*) array
,
115 "return value", "EOSHIFT");
118 if (unlikely (compile_options
.bounds_check
))
120 bounds_reduced_extents ((array_t
*) h
, (array_t
*) array
, which
,
121 "SHIFT argument", "EOSHIFT");
130 for (dim
= 0; dim
< GFC_DESCRIPTOR_RANK (array
); dim
++)
134 roffset
= GFC_DESCRIPTOR_STRIDE_BYTES(ret
,dim
);
137 soffset
= GFC_DESCRIPTOR_STRIDE_BYTES(array
,dim
);
140 len
= GFC_DESCRIPTOR_EXTENT(array
,dim
);
145 extent
[n
] = GFC_DESCRIPTOR_EXTENT(array
,dim
);
146 rstride
[n
] = GFC_DESCRIPTOR_STRIDE_BYTES(ret
,dim
);
147 sstride
[n
] = GFC_DESCRIPTOR_STRIDE_BYTES(array
,dim
);
149 hstride
[n
] = GFC_DESCRIPTOR_STRIDE(h
,n
);
151 bstride
[n
] = GFC_DESCRIPTOR_STRIDE_BYTES(bound
,n
);
163 if (bound
&& bstride
[0] == 0)
166 dim
= GFC_DESCRIPTOR_RANK (array
);
167 rstride0
= rstride
[0];
168 sstride0
= sstride
[0];
169 hstride0
= hstride
[0];
170 bstride0
= bstride
[0];
171 rptr
= ret
->base_addr
;
172 sptr
= array
->base_addr
;
175 bptr
= bound
->base_addr
;
181 /* Do the shift for this dimension. */
183 if (( sh
>= 0 ? sh
: -sh
) > len
)
189 delta
= (sh
>= 0) ? sh
: -sh
;
193 src
= &sptr
[delta
* soffset
];
199 dest
= &rptr
[delta
* roffset
];
201 for (n
= 0; n
< len
- delta
; n
++)
203 memcpy (dest
, src
, size
);
214 memcpy (dest
, bptr
, size
);
223 memset (dest
, filler
[0], size
);
225 for (i
= 0; i
< size
; i
+= filler_len
)
226 memcpy (&dest
[i
], filler
, filler_len
);
231 /* Advance to the next section. */
238 while (count
[n
] == extent
[n
])
240 /* When we get to the end of a dimension, reset it and increment
241 the next dimension. */
243 /* We could precalculate these products, but this is a less
244 frequently used path so probably not worth it. */
245 rptr
-= rstride
[n
] * extent
[n
];
246 sptr
-= sstride
[n
] * extent
[n
];
247 hptr
-= hstride
[n
] * extent
[n
];
248 bptr
-= bstride
[n
] * extent
[n
];
252 /* Break out of the loop. */
268 extern void eoshift3_16 (gfc_array_char
* const restrict
,
269 const gfc_array_char
* const restrict
,
270 const gfc_array_i16
* const restrict
,
271 const gfc_array_char
* const restrict
,
272 const GFC_INTEGER_16
*);
273 export_proto(eoshift3_16
);
276 eoshift3_16 (gfc_array_char
* const restrict ret
,
277 const gfc_array_char
* const restrict array
,
278 const gfc_array_i16
* const restrict h
,
279 const gfc_array_char
* const restrict bound
,
280 const GFC_INTEGER_16
* const restrict pwhich
)
282 eoshift3 (ret
, array
, h
, bound
, pwhich
, "\0", 1);
286 extern void eoshift3_16_char (gfc_array_char
* const restrict
,
288 const gfc_array_char
* const restrict
,
289 const gfc_array_i16
* const restrict
,
290 const gfc_array_char
* const restrict
,
291 const GFC_INTEGER_16
* const restrict
,
292 GFC_INTEGER_4
, GFC_INTEGER_4
);
293 export_proto(eoshift3_16_char
);
296 eoshift3_16_char (gfc_array_char
* const restrict ret
,
297 GFC_INTEGER_4 ret_length
__attribute__((unused
)),
298 const gfc_array_char
* const restrict array
,
299 const gfc_array_i16
* const restrict h
,
300 const gfc_array_char
* const restrict bound
,
301 const GFC_INTEGER_16
* const restrict pwhich
,
302 GFC_INTEGER_4 array_length
__attribute__((unused
)),
303 GFC_INTEGER_4 bound_length
__attribute__((unused
)))
305 eoshift3 (ret
, array
, h
, bound
, pwhich
, " ", 1);
309 extern void eoshift3_16_char4 (gfc_array_char
* const restrict
,
311 const gfc_array_char
* const restrict
,
312 const gfc_array_i16
* const restrict
,
313 const gfc_array_char
* const restrict
,
314 const GFC_INTEGER_16
* const restrict
,
315 GFC_INTEGER_4
, GFC_INTEGER_4
);
316 export_proto(eoshift3_16_char4
);
319 eoshift3_16_char4 (gfc_array_char
* const restrict ret
,
320 GFC_INTEGER_4 ret_length
__attribute__((unused
)),
321 const gfc_array_char
* const restrict array
,
322 const gfc_array_i16
* const restrict h
,
323 const gfc_array_char
* const restrict bound
,
324 const GFC_INTEGER_16
* const restrict pwhich
,
325 GFC_INTEGER_4 array_length
__attribute__((unused
)),
326 GFC_INTEGER_4 bound_length
__attribute__((unused
)))
328 static const gfc_char4_t space
= (unsigned char) ' ';
329 eoshift3 (ret
, array
, h
, bound
, pwhich
,
330 (const char *) &space
, sizeof (gfc_char4_t
));