JJSierraM/squiggle.c
1
0
Fork 0
forked from personal/squiggle.c
Code Issues Pull Requests Packages Projects Releases Wiki Activity

formatting pass nuclear probabilities

Browse Source
This commit is contained in:
NunoSempere 2023-07-23 23:54:15 +02:00
parent 8b32d2d9b6
commit b76a630d58
1 changed files with 23 additions and 28 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
examples/08_nuclear_war/example.c
View File

@ -3,12 +3,11 @@
#include <stdio.h> #include <stdio.h>
#include <stdlib.h> #include <stdlib.h>
double probability_of_dying_nuno(uint64_t* seed)
{
double probability_of_dying_nuno(uint64_t* seed){
double first_year_russian_nuclear_weapons = 1953; double first_year_russian_nuclear_weapons = 1953;
double current_year = 2023; double current_year = 2023;
double laplace_probability_nuclear_exchange_next_year = sample_beta(current_year-first_year_russian_nuclear_weapons, 0, seed); double laplace_probability_nuclear_exchange_next_year = sample_beta(current_year - first_year_russian_nuclear_weapons, 0, seed);
double london_hit_conditional_on_russia_nuclear_weapon_usage = sample_beta(7.67, 69.65, seed); double london_hit_conditional_on_russia_nuclear_weapon_usage = sample_beta(7.67, 69.65, seed);
// I.e., a beta distribution with a range of 0.05 to 0.16 into: https://nunosempere.com/blog/2023/03/15/fit-beta/ // I.e., a beta distribution with a range of 0.05 to 0.16 into: https://nunosempere.com/blog/2023/03/15/fit-beta/
@ -18,22 +17,19 @@ double probability_of_dying_nuno(uint64_t* seed){
// 0.2 to 0.8, i.e., 20% to 80%, again using the previous tool // 0.2 to 0.8, i.e., 20% to 80%, again using the previous tool
double proportion_which_die_if_bomb_drops_in_london = sample_beta(10.00, 2.45, seed); // 60% to 95% double proportion_which_die_if_bomb_drops_in_london = sample_beta(10.00, 2.45, seed); // 60% to 95%
double probability_of_dying = laplace_probability_nuclear_exchange_next_year * double probability_of_dying = laplace_probability_nuclear_exchange_next_year * london_hit_conditional_on_russia_nuclear_weapon_usage * informed_actor_not_able_to_escape * proportion_which_die_if_bomb_drops_in_london;
london_hit_conditional_on_russia_nuclear_weapon_usage *
informed_actor_not_able_to_escape *
proportion_which_die_if_bomb_drops_in_london;
return probability_of_dying; return probability_of_dying;
} }
double probability_of_dying_eli(uint64_t* seed){ double probability_of_dying_eli(uint64_t* seed)
double russia_nato_nuclear_exchange_in_next_month=sample_beta(1.30, 1182.99,seed); // .0001 to .003 {
double russia_nato_nuclear_exchange_in_next_month = sample_beta(1.30, 1182.99, seed); // .0001 to .003
double london_hit_conditional = sample_beta(3.47, 8.97, seed); // 0.1 to 0.5 double london_hit_conditional = sample_beta(3.47, 8.97, seed); // 0.1 to 0.5
double informed_actors_not_able_to_escape = sample_beta(2.73, 5.67, seed); // .1 to .6 double informed_actors_not_able_to_escape = sample_beta(2.73, 5.67, seed); // .1 to .6
double proportion_which_die_if_bomb_drops_in_london = sample_beta(3.00, 1.46, seed); // 0.3 to 0.95; double proportion_which_die_if_bomb_drops_in_london = sample_beta(3.00, 1.46, seed); // 0.3 to 0.95;
double probability_of_dying = russia_nato_nuclear_exchange_in_next_month*london_hit_conditional*informed_actors_not_able_to_escape*proportion_which_die_if_bomb_drops_in_london; double probability_of_dying = russia_nato_nuclear_exchange_in_next_month * london_hit_conditional * informed_actors_not_able_to_escape * proportion_which_die_if_bomb_drops_in_london;
return probability_of_dying; return probability_of_dying;
} }
int main() int main()
@ -42,6 +38,5 @@ int main()
uint64_t* seed = malloc(sizeof(uint64_t)); uint64_t* seed = malloc(sizeof(uint64_t));
*seed = 1000; // xorshift can't start with 0 *seed = 1000; // xorshift can't start with 0
free(seed); free(seed);
} }