Function fpdiv

src/big-decimal/utils.ts:70–72

(a: bigint, b: bigint, bits: number)

/** Fixed-point divide on the base-2 grid: (a << bits) / b */
export function fpdiv(a: bigint, b: bigint, bits: number): bigint {
  return (a << BigInt(bits)) / b;
}

/**
 * Fixed-point square root via Newton/Heron iteration on the base-2 grid.