/****************************************************************************/
/*	Copyright (c) 2009, Bendor Research Pty. Ltd. All rights reserved.		*/
/*	This software is released under the GNU General Public License, version	*/
/*	3, with the exemptions that 1) the header files installed by this 	 	*/
/*	package can be included in closed-source code and 2) the libraries		*/
/*	built by this package can be statically linked into a closed-source 	*/
/*	binary. See the file LICENSE for the details.							*/
/*	This software is provided "AS IS", without any warranty or guarantee.	*/
/*	By using this software you acknowledge that neither Bendor Research nor	*/
/*	any of its personnel or affiliates can be held liable for any damages,	*/
/*	including, but not limited to, loss of business, loss of money, loss of */
/*	life, loss of reputation, arising out of this software being faulty or	*/
/*	not performing as advertised, even if Bendor Research was aware of such	*/
/*	faults in the software. You use this software at your own risk, any and */
/*	all liabilities that arise out of your using this software are strictly	*/
/*	yours.																	*/
/****************************************************************************/

#include <arm_func.h>

/*
*	Division by constant. At entry R0 contains a shift amount on the lowest
*	bits and a signum on its top bit. R1 is the multiplier and R2 is the
*	value to divide. We do the following:
*
*	int	sdiv32( int sh1, int sig, unsigned mul, int x )
*	{
*	int	t;
*
*		t = ( (long long) x * mul ) >> 32;
*		t = ( t >> sh1 ) + ( x < 0 ? 1 : 0 );
*		t = sig ? -t : +t;
*		return( t );
*	}
*
*	The math details of the algorithm and the way sh1, sig and mul are
*	calculated can be found in the article:
*
*	Torbjörn Granlund, Peter L. Montgomery:
*	Division by Invariant Integers using Multiplication
*/

	afunc	sdiv32
	umull	r12,r3,r2,r1
	mov		r3,r3,asr r0
	sub		r3,r3,r2,asr #31
	eor		r3,r3,r0,asr #31
	sub		r0,r3,r0,asr #31
	bx		lr

