git.pond.sub.org Git - empserver/blob - src/lib/update/nat.c

   1 /*
   2  *  Empire - A multi-player, client/server Internet based war game.
   3  *  Copyright (C) 1986-2007, Dave Pare, Jeff Bailey, Thomas Ruschak,
   4  *                           Ken Stevens, Steve McClure
   5  *
   6  *  This program is free software; you can redistribute it and/or modify
   7  *  it under the terms of the GNU General Public License as published by
   8  *  the Free Software Foundation; either version 2 of the License, or
   9  *  (at your option) any later version.
  10  *
  11  *  This program is distributed in the hope that it will be useful,
  12  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14  *  GNU General Public License for more details.
  15  *
  16  *  You should have received a copy of the GNU General Public License
  17  *  along with this program; if not, write to the Free Software
  18  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  19  *
  20  *  ---
  21  *
  22  *  See files README, COPYING and CREDITS in the root of the source
  23  *  tree for related information and legal notices.  It is expected
  24  *  that future projects/authors will amend these files as needed.
  25  *
  26  *  ---
  27  *
  28  *  nat.c: Accumulate tech, edu, research and happiness.
  29  *
  30  *  Known contributors to this file:
  31  *     Dave Pare, 1989
  32  *     Steve McClure, 1997
  33  */
  34
  35 #include <config.h>
  36
  37 #include <math.h>
  38 #include "budg.h"
  39 #include "item.h"
  40 #include "update.h"
  41
  42 /*
  43  * hap and edu avg mean that the weight on current happiness is
  44  *  (cur_hap * hap_avg + hap_prod * etu) / (hap_avg + etu);
  45  * same for education.
  46  * right now, happiness has 1 day (48 etu) average, prod of 10 from
  47  * initial level of 0 yields (1) 1.42, (6) 6.03, (12) 8.42, (18) 9.37
  48  *
  49  * education has 4 day (192 etu) average, prod of 10 from initial
  50  * level of 0 yields (1) 0.4, (6) 2.2, (12) 3.9, (18) 5.2.
  51  */
  52
  53 static void share_incr(double *, double *);
  54
  55 /*
  56  * for values below the "easy level" values, production is
  57  * as normal.  For values above "easy", production gets harder
  58  * based on an equation in "limit_level()" in update/nat.c.
  59  * Basically, the smaller the the values for "level_log", the
  60  * smaller return on investment above level_easy[] values.
  61  */
  62 /*
  63  * Damn! I hate this, but ...
  64  * The values here for tech are *not* the real ones.
  65  * They are changed later in the limit_level routine.
  66  */
  67                         /*tech   res   edu   hap */
  68 static float level_easy[4] = { 0.75, 0.75, 5.00, 5.00 };
  69 static float level_log[4] = { 1.75, 2.00, 4.00, 6.00 };
  70
  71 float levels[MAXNOC][4];
  72
  73 /*
  74  * technique to limit the sharpers who turn entire countries
  75  * into tech plants overnight...
  76  */
  77
  78 double
  79 logx(double d, double base)
  80 {
  81     if (base == 1.0)
  82         return d;
  83     return log10(d) / log10(base);
  84 }
  85
  86 static double
  87 limit_level(double level, int type, int flag)
  88 {
  89     double above_easy;
  90     double above;
  91     double logbase;
  92     double easy;
  93
  94 /*
  95  * Begin ugly hack.
  96  */
  97     level_easy[0] = easy_tech;
  98     level_log[0] = tech_log_base;
  99 /*
 100  * End ugly hack.
 101  */
 102
 103     if (level > level_easy[type]) {
 104         logbase = level_log[type];
 105         easy = level_easy[type];
 106         above_easy = level - easy;
 107         if (flag)
 108             above = above_easy / logx(logbase + above_easy, logbase);
 109         else
 110             above = logx(above_easy + 1.0, logbase);
 111         if (above > 250)
 112             above = 250;
 113         return above < 0 ? easy : easy + above;
 114     } else
 115         return level;
 116 }
 117
 118 void
 119 prod_nat(int etu)
 120 {
 121     struct natstr *np;
 122     float hap;
 123     float edu;
 124     float hap_edu;
 125     long pop;
 126     double rlev;
 127     double tlev;
 128     double tech[MAXNOC];
 129     double res[MAXNOC];
 130     double newvalue;
 131     natid n;
 132     int cn;
 133     struct natstr *cnp;
 134
 135     for (n = 0; NULL != (np = getnatp(n)); n++) {
 136         grant_btus(np, etu_per_update - np->nat_access);
 137         np->nat_access = 0;
 138         if (np->nat_stat < STAT_ACTIVE)
 139             continue;
 140         /*
 141          * hap_edu: the more education people have, the
 142          * more happiness they want.
 143          */
 144         hap_edu = np->nat_level[NAT_ELEV];
 145         hap_edu = 1.5 - ((hap_edu + 10.0) / (hap_edu + 20.0));
 146         pop = pops[n] + 1;
 147         /*
 148          * get per-population happiness and education
 149          * see what the total per-civilian production is
 150          * for this time period.
 151          */
 152         hap = levels[n][NAT_HLEV] * hap_edu * hap_cons /
 153             ((float)pop * etu);
 154         edu = levels[n][NAT_ELEV] * edu_cons / ((float)pop * etu);
 155         wu((natid)0, n, "%3.0f happiness, %3.0f education produced\n",
 156            levels[n][NAT_HLEV], levels[n][NAT_ELEV]);
 157         hap = limit_level(hap, NAT_HLEV, 1);
 158         edu = limit_level(edu, NAT_ELEV, 1);
 159         /*
 160          * change the "moving average"...old happiness and
 161          * education levels are weighted heavier than current
 162          * production.
 163          */
 164         newvalue = (np->nat_level[NAT_HLEV] * hap_avg + hap * etu) /
 165             (hap_avg + etu);
 166         np->nat_level[NAT_HLEV] = newvalue;
 167         newvalue = (np->nat_level[NAT_ELEV] * edu_avg + edu * etu) /
 168             (edu_avg + etu);
 169         np->nat_level[NAT_ELEV] = newvalue;
 170         /*
 171          * limit tech/research production
 172          */
 173         levels[n][NAT_TLEV] =
 174             limit_level(levels[n][NAT_TLEV] / 1, NAT_TLEV, 0) * 1;
 175         levels[n][NAT_RLEV] =
 176             limit_level(levels[n][NAT_RLEV] / 1, NAT_RLEV, 0) * 1;
 177         wu((natid)0, n,
 178            "total pop was %ld, yielding %4.2f hap, %4.2f edu\n",
 179            pop - 1, hap, edu);
 180     }
 181     if (ally_factor > 0.0)
 182         share_incr(res, tech);
 183     else {
 184         memset(res, 0, sizeof(res));
 185         memset(tech, 0, sizeof(tech));
 186     }
 187     for (n = 0; NULL != (np = getnatp(n)); n++) {
 188         if (np->nat_stat < STAT_ACTIVE)
 189             continue;
 190         tlev = levels[n][NAT_TLEV];
 191         rlev = levels[n][NAT_RLEV];
 192         if (tech[n] != 0.0 || res[n] != 0.0) {
 193             wu((natid)0, n,
 194                "%5.4f technology (%5.4f + %5.4f), ",
 195                tlev + tech[n], tlev, tech[n]);
 196             wu((natid)0, n,
 197                "%5.4f research (%5.4f + %5.4f) produced\n",
 198                rlev + res[n], rlev, res[n]);
 199         } else
 200             wu((natid)0, n,
 201                "%5.4f tech, %5.4f research produced\n", tlev, rlev);
 202         rlev += res[n];
 203         tlev += tech[n];
 204         if (rlev != 0.0)
 205             np->nat_level[NAT_RLEV] += rlev;
 206         if (tlev != 0.0)
 207             np->nat_level[NAT_TLEV] += tlev;
 208         if ((sea_money[n] != 0) || (air_money[n] != 0) ||
 209             (lnd_money[n] != 0))
 210             wu((natid)0, n,
 211                "Army delta $%ld, Navy delta $%ld, Air force delta $%ld\n",
 212                lnd_money[n], sea_money[n], air_money[n]);
 213         wu((natid)0, n, "money delta was $%ld for this update\n",
 214            np->nat_money - money[n]);
 215         if (opt_LOSE_CONTACT) {
 216             for (cn = 1; cn < MAXNOC; cn++) {
 217                 if ((cnp = getnatp(cn)) != NULL)
 218                     agecontact(cnp);
 219             }
 220         }
 221     }
 222 }
 223
 224 /*
 225  * find out everyones increment
 226  */
 227 static void
 228 share_incr(double *res, double *tech)
 229 {
 230     struct natstr *np;
 231     struct natstr *other;
 232     natid i;
 233     natid j;
 234     int rnc;
 235     int tnc;
 236
 237     for (i = 0; NULL != (np = getnatp(i)); i++) {
 238         res[i] = tech[i] = 0.0;
 239         if (np->nat_stat < STAT_SANCT || np->nat_stat == STAT_GOD)
 240             continue;
 241         rnc = tnc = 0;
 242         for (j = 0; NULL != (other = getnatp(j)); j++) {
 243             if (j == i)
 244                 continue;
 245             if (other->nat_stat != STAT_ACTIVE)
 246                 continue;
 247             if (opt_HIDDEN) {
 248                 if (!getcontact(np, j))
 249                     continue;
 250             }
 251             if (!opt_ALL_BLEED) {
 252                 if (getrel(np, j) != ALLIED)
 253                     continue;
 254                 if (getrel(other, i) != ALLIED)
 255                     continue;
 256                 res[i] += levels[j][NAT_RLEV];
 257                 tech[i] += levels[j][NAT_TLEV];
 258                 rnc++;
 259                 tnc++;
 260             } else {
 261                 if (levels[j][NAT_TLEV] > 0.001) {
 262                     tech[i] += levels[j][NAT_TLEV];
 263                     tnc++;
 264                 }
 265                 if (levels[j][NAT_RLEV] > 0.001) {
 266                     res[i] += levels[j][NAT_RLEV];
 267                     rnc++;
 268                 }
 269             }
 270         }
 271         if (rnc > 0) {
 272             res[i] /= rnc * ally_factor;
 273         }
 274         if (tnc > 0) {
 275             tech[i] /= tnc * ally_factor;
 276         }
 277 /*              logerror("Country #%d gets %g res from %d allies, %g tech from %d allies", i, res[i], rnc, tech[i], tnc);*/
 278     }
 279 }