flint-sys 0.9.0

Bindings to the FLINT C library
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34

/*
    Copyright (C) 2021 Daniel Schultz

    This file is part of FLINT.

    FLINT is free software: you can redistribute it and/or modify it under
    the terms of the GNU Lesser General Public License (LGPL) as published
    by the Free Software Foundation; either version 3 of the License, or
    (at your option) any later version.  See <https://www.gnu.org/licenses/>.
*/

#include "thread_pool.h"

/*
    if f is continuous with
        f(0) = 0
        f'(x) = alpha for 0 < x < a
        f'(x) = beta for a < x < a + b
    return solution for x to f(x) = (a*alpha + b*beta)*yn/yd
    for 0 <= yn/yd <= 1
*/
ulong _thread_pool_find_work_2(
    ulong a, ulong alpha,
    ulong b, ulong beta,
    ulong yn, ulong yd)
{
    /* very low priority TODO: this can overflow only in very extreme cases */
    ulong y = yn*(a*alpha + b*beta)/yd;

    if (y <= a*alpha)
        return y/alpha;

    return a + (y - a*alpha)/beta;
}