/*
 * Copyright 2012-15 Advanced Micro Devices, Inc.
 *
 * Permission is hereby granted, free of charge, to any person obtaining a
 * copy of this software and associated documentation files (the "Software"),
 * to deal in the Software without restriction, including without limitation
 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
 * and/or sell copies of the Software, and to permit persons to whom the
 * Software is furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in
 * all copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
 * THE COPYRIGHT HOLDER(S) OR AUTHOR(S) BE LIABLE FOR ANY CLAIM, DAMAGES OR
 * OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
 * OTHER DEALINGS IN THE SOFTWARE.
 *
 * Authors: AMD
 *
 */

#include "dm_services.h"
#include "include/fixed32_32.h"

static uint64_t u64_div(uint64_t n, uint64_t d)
{
	uint32_t i = 0;
	uint64_t r;
	uint64_t q = div64_u64_rem(n, d, &r);

	for (i = 0; i < 32; ++i) {
		uint64_t sbit = q & (1ULL<<63);

		r <<= 1;
		r |= sbit ? 1 : 0;
		q <<= 1;
		if (r >= d) {
			r -= d;
			q |= 1;
		}
	}

	if (2*r >= d)
		q += 1;
	return q;
}

struct fixed32_32 dal_fixed32_32_from_fraction(uint32_t n, uint32_t d)
{
	struct fixed32_32 fx;

	fx.value = u64_div((uint64_t)n << 32, (uint64_t)d << 32);
	return fx;
}

struct fixed32_32 dal_fixed32_32_add(
	struct fixed32_32 lhs,
	struct fixed32_32 rhs)
{
	struct fixed32_32 fx = {lhs.value + rhs.value};

	ASSERT(fx.value >= rhs.value);
	return fx;
}

struct fixed32_32 dal_fixed32_32_add_int(struct fixed32_32 lhs, uint32_t rhs)
{
	struct fixed32_32 fx = {lhs.value + ((uint64_t)rhs << 32)};

	ASSERT(fx.value >= (uint64_t)rhs << 32);
	return fx;

}
struct fixed32_32 dal_fixed32_32_sub(
	struct fixed32_32 lhs,
	struct fixed32_32 rhs)
{
	struct fixed32_32 fx;

	ASSERT(lhs.value >= rhs.value);
	fx.value = lhs.value - rhs.value;
	return fx;
}

struct fixed32_32 dal_fixed32_32_sub_int(struct fixed32_32 lhs, uint32_t rhs)
{
	struct fixed32_32 fx;

	ASSERT(lhs.value >= ((uint64_t)rhs<<32));
	fx.value = lhs.value - ((uint64_t)rhs<<32);
	return fx;
}

struct fixed32_32 dal_fixed32_32_mul(
	struct fixed32_32 lhs,
	struct fixed32_32 rhs)
{
	struct fixed32_32 fx;
	uint64_t lhs_int = lhs.value>>32;
	uint64_t lhs_frac = (uint32_t)lhs.value;
	uint64_t rhs_int = rhs.value>>32;
	uint64_t rhs_frac = (uint32_t)rhs.value;
	uint64_t ahbh = lhs_int * rhs_int;
	uint64_t ahbl = lhs_int * rhs_frac;
	uint64_t albh = lhs_frac * rhs_int;
	uint64_t albl = lhs_frac * rhs_frac;

	ASSERT((ahbh>>32) == 0);

	fx.value = (ahbh<<32) + ahbl + albh + (albl>>32);
	return fx;

}

struct fixed32_32 dal_fixed32_32_mul_int(struct fixed32_32 lhs, uint32_t rhs)
{
	struct fixed32_32 fx;
	uint64_t lhsi = (lhs.value>>32) * (uint64_t)rhs;
	uint64_t lhsf;

	ASSERT((lhsi>>32) == 0);
	lhsf = ((uint32_t)lhs.value) * (uint64_t)rhs;
	ASSERT((lhsi<<32) + lhsf >= lhsf);
	fx.value = (lhsi<<32) + lhsf;
	return fx;
}

struct fixed32_32 dal_fixed32_32_div(
	struct fixed32_32 lhs,
	struct fixed32_32 rhs)
{
	struct fixed32_32 fx;

	fx.value = u64_div(lhs.value, rhs.value);
	return fx;
}

struct fixed32_32 dal_fixed32_32_div_int(struct fixed32_32 lhs, uint32_t rhs)
{
	struct fixed32_32 fx;

	fx.value = u64_div(lhs.value, (uint64_t)rhs << 32);
	return fx;
}

uint32_t dal_fixed32_32_ceil(struct fixed32_32 v)
{
	ASSERT((uint32_t)v.value ? (v.value >> 32) + 1 >= 1 : true);
	return (v.value>>32) + ((uint32_t)v.value ? 1 : 0);
}

uint32_t dal_fixed32_32_round(struct fixed32_32 v)
{
	ASSERT(v.value + (1ULL<<31) >= (1ULL<<31));
	return (v.value + (1ULL<<31))>>32;
}

