test_logprod.c

#include "hope.h"
#include "liknorm.h"
#include "stdio.h"

int main()
{

    struct LikNormMachine *machine = liknorm_create_machine(10);
    double y = 5.0;
    double tau = 0.15;
    double eta = 2.1;
    liknorm_set_poisson(machine, y);
    liknorm_set_prior(machine, tau, eta);

    CLOSE(liknorm_logprod(machine, 700), -1.0142320547350045e+304);
    CLOSE(liknorm_logprod(machine, 50), -5.184705528587072e+21);
    CLOSE(liknorm_logprod(machine, 12.0), -162701.74640927234);
    CLOSE(liknorm_logprod(machine, 0.59), -18.996086183827522);
    CLOSE(liknorm_logprod(machine, -3.0), -43.37977733679752);
    CLOSE(liknorm_logprod(machine, -10.9), -107.65575872666366);
    CLOSE(liknorm_logprod(machine, -20), -193.35499027049082);
    CLOSE(liknorm_logprod(machine, -50), -563.8549902684297);
    CLOSE(liknorm_logprod(machine, -500), -22321.35499026843);

    double k = 1;
    double n = 10;
    liknorm_set_binomial(machine, k, n);
    CLOSE(liknorm_logprod(machine, 0.0), -21.196385238253026);
    CLOSE(liknorm_logprod(machine, 1.0), -24.372530307835795);
    CLOSE(liknorm_logprod(machine, 10.0), -90.76536742165445);
    CLOSE(liknorm_logprod(machine, 50.0), -546.7649134326535);
    CLOSE(liknorm_logprod(machine, 100.0), -1454.2649134326534);
    CLOSE(liknorm_logprod(machine, 500.0), -22214.264913432653);
    CLOSE(liknorm_logprod(machine, 5000.0), -1909514.2649134325);

    CLOSE(liknorm_logprod(machine, -10.0), -52.76536742164575);
    CLOSE(liknorm_logprod(machine, -50.0), -356.7649134326536);
    CLOSE(liknorm_logprod(machine, -100.0), -1074.2649134326537);
    CLOSE(liknorm_logprod(machine, -500.0), -20314.264913432653);
    CLOSE(liknorm_logprod(machine, -5000.0), -1890514.2649134325);

    k = 3;
    double r = 5;
    liknorm_set_prior(machine, 1, 0);
    liknorm_set_nbinomial(machine, k, r);
    CLOSE(liknorm_logprod(machine, -0.10536051565782628), -9.19814790278881);

    liknorm_destroy_machine(machine);

    return hope_status();
}