/* Optimized version of the standard memset() function.
   This file is part of the GNU C Library.
   Copyright (C) 2000, 2001, 2002, 2003 Free Software Foundation, Inc.
   Contributed by Dan Pop for Itanium <Dan.Pop@cern.ch>.
   Rewritten for McKinley by Sverre Jarp, HP Labs/CERN <Sverre.Jarp@cern.ch>

   The GNU C Library is free software; you can redistribute it and/or
   modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   version 2.1 of the License, or (at your option) any later version.

   The GNU C Library is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public
   License along with the GNU C Library; if not, see
   <http://www.gnu.org/licenses/>.  */

/* Return: dest

   Inputs:
        in0:    dest
        in1:    value
        in2:    count

   The algorithm is fairly straightforward: set byte by byte until we
   we get to a 16B-aligned address, then loop on 128 B chunks using an
   early store as prefetching, then loop on 32B chucks, then clear remaining
   words, finally clear remaining bytes.
   Since a stf.spill f0 can store 16B in one go, we use this instruction
   to get peak speed when value = 0.  */

#include <sysdep.h>
#undef ret

#define dest		in0
#define value		in1
#define	cnt		in2

#define tmp		r31
#define save_lc		r30
#define ptr0		r29
#define ptr1		r28
#define ptr2		r27
#define ptr3		r26
#define ptr9		r24
#define	loopcnt		r23
#define linecnt		r22
#define bytecnt		r21

#define fvalue		f6

/* This routine uses only scratch predicate registers (p6 - p15) */
#define p_scr		p6	/* default register for same-cycle branches */
#define p_nz		p7
#define p_zr		p8
#define p_unalgn	p9
#define p_y		p11
#define p_n		p12
#define p_yy		p13
#define p_nn		p14

#define movi0		mov

#define MIN1		15
#define MIN1P1HALF	8
#define LINE_SIZE	128
#define LSIZE_SH        7			/* shift amount */
#define PREF_AHEAD	8

#define USE_FLP
#if defined(USE_INT)
#define store		st8
#define myval           value
#elif defined(USE_FLP)
#define store		stf8
#define myval		fvalue
#endif

.align	64
ENTRY(memset)
{ .mmi
	.prologue
	alloc	tmp = ar.pfs, 3, 0, 0, 0
	lfetch.nt1 [dest]
	.save   ar.lc, save_lc
	movi0	save_lc = ar.lc
} { .mmi
	.body
	mov	ret0 = dest		/* return value */
	cmp.ne	p_nz, p_zr = value, r0	/* use stf.spill if value is zero */
	cmp.eq	p_scr, p0 = cnt, r0
;; }
{ .mmi
	and	ptr2 = -(MIN1+1), dest	/* aligned address */
	and	tmp = MIN1, dest	/* prepare to check for alignment */
	tbit.nz p_y, p_n = dest, 0	/* Do we have an odd address? (M_B_U) */
} { .mib
	mov	ptr1 = dest
	mux1	value = value, @brcst	/* create 8 identical bytes in word */
(p_scr)	br.ret.dpnt.many rp		/* return immediately if count = 0 */
;; }
{ .mib
	cmp.ne	p_unalgn, p0 = tmp, r0
} { .mib				/* NB: # of bytes to move is 1 higher */
	sub	bytecnt = (MIN1+1), tmp	/*     than loopcnt */
	cmp.gt	p_scr, p0 = 16, cnt		/* is it a minimalistic task? */
(p_scr)	br.cond.dptk.many .move_bytes_unaligned	/* go move just a few (M_B_U) */
;; }
{ .mmi
(p_unalgn) add	ptr1 = (MIN1+1), ptr2		/* after alignment */
(p_unalgn) add	ptr2 = MIN1P1HALF, ptr2		/* after alignment */
(p_unalgn) tbit.nz.unc p_y, p_n = bytecnt, 3	/* should we do a st8 ? */
;; }
{ .mib
(p_y)	add	cnt = -8, cnt
(p_unalgn) tbit.nz.unc p_yy, p_nn = bytecnt, 2	/* should we do a st4 ? */
} { .mib
(p_y)	st8	[ptr2] = value, -4
(p_n)	add	ptr2 = 4, ptr2
;; }
{ .mib
(p_yy)	add	cnt = -4, cnt
(p_unalgn) tbit.nz.unc p_y, p_n = bytecnt, 1	/* should we do a st2 ? */
} { .mib
(p_yy)	st4	[ptr2] = value, -2
(p_nn)	add	ptr2 = 2, ptr2
;; }
{ .mmi
	mov	tmp = LINE_SIZE+1		/* for compare */
(p_y)	add	cnt = -2, cnt
(p_unalgn) tbit.nz.unc p_yy, p_nn = bytecnt, 0	/* should we do a st1 ? */
} { .mmi
	setf.sig fvalue=value			/* transfer value to FLP side */
(p_y)	st2	[ptr2] = value, -1
(p_n)	add	ptr2 = 1, ptr2
;; }

{ .mmi
(p_yy)	st1	[ptr2] = value
	cmp.gt	p_scr, p0 = tmp, cnt		/* is it a minimalistic task? */
} { .mbb
(p_yy)	add	cnt = -1, cnt
(p_scr)	br.cond.dpnt.many .fraction_of_line	/* go move just a few */
;; }

{ .mib
	nop.m 0
	shr.u	linecnt = cnt, LSIZE_SH
(p_zr)	br.cond.dptk.many .l1b			/* Jump to use stf.spill */
;; }

#ifndef GAS_ALIGN_BREAKS_UNWIND_INFO
	.align 32 /* --------  L1A: store ahead into cache lines; fill later */
#endif
{ .mmi
	and	tmp = -(LINE_SIZE), cnt		/* compute end of range */
	mov	ptr9 = ptr1			/* used for prefetching */
	and	cnt = (LINE_SIZE-1), cnt	/* remainder */
} { .mmi
	mov	loopcnt = PREF_AHEAD-1		/* default prefetch loop */
	cmp.gt	p_scr, p0 = PREF_AHEAD, linecnt	/* check against actual value */
;; }
{ .mmi
(p_scr)	add	loopcnt = -1, linecnt		/* start of stores */
	add	ptr2 = 8, ptr1			/* (beyond prefetch stores) */
	add	ptr1 = tmp, ptr1		/* first address beyond total */
;; }						/* range */
{ .mmi
	add	tmp = -1, linecnt		/* next loop count */
	movi0	ar.lc = loopcnt
;; }
.pref_l1a:
{ .mib
	store [ptr9] = myval, 128	/* Do stores one cache line apart */
	nop.i	0
	br.cloop.dptk.few .pref_l1a
;; }
{ .mmi
	add	ptr0 = 16, ptr2		/* Two stores in parallel */
	movi0	ar.lc = tmp
;; }
.l1ax:
 { .mmi
	store [ptr2] = myval, 8
	store [ptr0] = myval, 8
 ;; }
 { .mmi
	store [ptr2] = myval, 24
	store [ptr0] = myval, 24
 ;; }
 { .mmi
	store [ptr2] = myval, 8
	store [ptr0] = myval, 8
 ;; }
 { .mmi
	store [ptr2] = myval, 24
	store [ptr0] = myval, 24
 ;; }
 { .mmi
	store [ptr2] = myval, 8
	store [ptr0] = myval, 8
 ;; }
 { .mmi
	store [ptr2] = myval, 24
	store [ptr0] = myval, 24
 ;; }
 { .mmi
	store [ptr2] = myval, 8
	store [ptr0] = myval, 32
	cmp.lt	p_scr, p0 = ptr9, ptr1	/* do we need more prefetching? */
 ;; }
{ .mmb
	store [ptr2] = myval, 24
(p_scr)	store [ptr9] = myval, 128
	br.cloop.dptk.few .l1ax
;; }
{ .mbb
	cmp.le  p_scr, p0 = 8, cnt		/* just a few bytes left ? */
(p_scr) br.cond.dpnt.many  .fraction_of_line	/* Branch no. 2 */
	br.cond.dpnt.many  .move_bytes_from_alignment	/* Branch no. 3 */
;; }

#ifdef GAS_ALIGN_BREAKS_UNWIND_INFO
	{ nop 0 }
#else
	.align 32
#endif
.l1b:	/* ------------------  L1B: store ahead into cache lines; fill later */
{ .mmi
	and	tmp = -(LINE_SIZE), cnt		/* compute end of range */
	mov	ptr9 = ptr1			/* used for prefetching */
	and	cnt = (LINE_SIZE-1), cnt	/* remainder */
} { .mmi
	mov	loopcnt = PREF_AHEAD-1		/* default prefetch loop */
	cmp.gt	p_scr, p0 = PREF_AHEAD, linecnt	/* check against actual value */
;; }
{ .mmi
(p_scr)	add	loopcnt = -1, linecnt
	add	ptr2 = 16, ptr1	/* start of stores (beyond prefetch stores) */
	add	ptr1 = tmp, ptr1	/* first address beyond total range */
;; }
{ .mmi
	add	tmp = -1, linecnt	/* next loop count */
	movi0	ar.lc = loopcnt
;; }
.pref_l1b:
{ .mib
	stf.spill [ptr9] = f0, 128	/* Do stores one cache line apart */
	nop.i   0
	br.cloop.dptk.few .pref_l1b
;; }
{ .mmi
	add	ptr0 = 16, ptr2		/* Two stores in parallel */
	movi0	ar.lc = tmp
;; }
.l1bx:
 { .mmi
	stf.spill [ptr2] = f0, 32
	stf.spill [ptr0] = f0, 32
 ;; }
 { .mmi
	stf.spill [ptr2] = f0, 32
	stf.spill [ptr0] = f0, 32
 ;; }
 { .mmi
	stf.spill [ptr2] = f0, 32
	stf.spill [ptr0] = f0, 64
	cmp.lt	p_scr, p0 = ptr9, ptr1	/* do we need more prefetching? */
 ;; }
{ .mmb
	stf.spill [ptr2] = f0, 32
(p_scr)	stf.spill [ptr9] = f0, 128
	br.cloop.dptk.few .l1bx
;; }
{ .mib
	cmp.gt  p_scr, p0 = 8, cnt	/* just a few bytes left ? */
(p_scr)	br.cond.dpnt.many  .move_bytes_from_alignment
;; }

.fraction_of_line:
{ .mib
	add	ptr2 = 16, ptr1
	shr.u	loopcnt = cnt, 5	/* loopcnt = cnt / 32 */
;; }
{ .mib
	cmp.eq	p_scr, p0 = loopcnt, r0
	add	loopcnt = -1, loopcnt
(p_scr)	br.cond.dpnt.many store_words
;; }
{ .mib
	and	cnt = 0x1f, cnt		/* compute the remaining cnt */
	movi0   ar.lc = loopcnt
;; }
#ifndef GAS_ALIGN_BREAKS_UNWIND_INFO
	.align 32
#endif
.l2:	/* ----------------------------  L2A:  store 32B in 2 cycles */
{ .mmb
	store	[ptr1] = myval, 8
	store	[ptr2] = myval, 8
;; } { .mmb
	store	[ptr1] = myval, 24
	store	[ptr2] = myval, 24
	br.cloop.dptk.many .l2
;; }
store_words:
{ .mib
	cmp.gt	p_scr, p0 = 8, cnt		/* just a few bytes left ? */
(p_scr)	br.cond.dpnt.many .move_bytes_from_alignment	/* Branch */
;; }

{ .mmi
	store	[ptr1] = myval, 8		/* store */
	cmp.le	p_y, p_n = 16, cnt		/* */
	add	cnt = -8, cnt			/* subtract */
;; }
{ .mmi
(p_y)	store	[ptr1] = myval, 8		/* store */
(p_y)	cmp.le.unc p_yy, p_nn = 16, cnt		/* */
(p_y)	add	cnt = -8, cnt			/* subtract */
;; }
{ .mmi						/* store */
(p_yy)	store	[ptr1] = myval, 8		/* */
(p_yy)	add	cnt = -8, cnt			/* subtract */
;; }

.move_bytes_from_alignment:
{ .mib
	cmp.eq	p_scr, p0 = cnt, r0
	tbit.nz.unc p_y, p0 = cnt, 2	/* should we terminate with a st4 ? */
(p_scr)	br.cond.dpnt.few .restore_and_exit
;; }
{ .mib
(p_y)	st4	[ptr1] = value, 4
	tbit.nz.unc p_yy, p0 = cnt, 1	/* should we terminate with a st2 ? */
;; }
{ .mib
(p_yy)	st2	[ptr1] = value, 2
	tbit.nz.unc p_y, p0 = cnt, 0
;; }

{ .mib
(p_y)	st1	[ptr1] = value
;; }
.restore_and_exit:
{ .mib
	nop.m	0
	movi0	ar.lc = save_lc
	br.ret.sptk.many rp
;; }

.move_bytes_unaligned:
{ .mmi
       .pred.rel "mutex",p_y, p_n
       .pred.rel "mutex",p_yy, p_nn
(p_n)	cmp.le  p_yy, p_nn = 4, cnt
(p_y)	cmp.le  p_yy, p_nn = 5, cnt
(p_n)	add	ptr2 = 2, ptr1
} { .mmi
(p_y)	add	ptr2 = 3, ptr1
(p_y)	st1	[ptr1] = value, 1	/* fill 1 (odd-aligned) byte */
(p_y)	add	cnt = -1, cnt		/* [15, 14 (or less) left] */
;; }
{ .mmi
(p_yy)	cmp.le.unc p_y, p0 = 8, cnt
	add	ptr3 = ptr1, cnt	/* prepare last store */
	movi0	ar.lc = save_lc
} { .mmi
(p_yy)	st2	[ptr1] = value, 4	/* fill 2 (aligned) bytes */
(p_yy)	st2	[ptr2] = value, 4	/* fill 2 (aligned) bytes */
(p_yy)	add	cnt = -4, cnt		/* [11, 10 (o less) left] */
;; }
{ .mmi
(p_y)	cmp.le.unc p_yy, p0 = 8, cnt
	add	ptr3 = -1, ptr3		/* last store */
	tbit.nz p_scr, p0 = cnt, 1	/* will there be a st2 at the end ? */
} { .mmi
(p_y)	st2	[ptr1] = value, 4	/* fill 2 (aligned) bytes */
(p_y)	st2	[ptr2] = value, 4	/* fill 2 (aligned) bytes */
(p_y)	add	cnt = -4, cnt		/* [7, 6 (or less) left] */
;; }
{ .mmi
(p_yy)	st2	[ptr1] = value, 4	/* fill 2 (aligned) bytes */
(p_yy)	st2	[ptr2] = value, 4	/* fill 2 (aligned) bytes */
					/* [3, 2 (or less) left] */
	tbit.nz p_y, p0 = cnt, 0	/* will there be a st1 at the end ? */
} { .mmi
(p_yy)	add	cnt = -4, cnt
;; }
{ .mmb
(p_scr)	st2	[ptr1] = value		/* fill 2 (aligned) bytes */
(p_y)	st1	[ptr3] = value		/* fill last byte (using ptr3) */
	br.ret.sptk.many rp
;; }
END(memset)
libc_hidden_def (memset)