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-2000, 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 the "LEGAL", "LICENSE", "CREDITS" and "README" files for all the
  23  *  related information and legal notices. It is expected that any future
  24  *  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 <math.h>
  36 #include "misc.h"
  37 #include "var.h"
  38 #include "sect.h"
  39 #include "nat.h"
  40 #include "item.h"
  41 #include "news.h"
  42 #include "file.h"
  43 #include "xy.h"
  44 #include "optlist.h"
  45 #include "update.h"
  46 #include "subs.h"
  47
  48 float   levels[MAXNOC][4];
  49
  50 /*
  51  * hap and edu avg mean that the weight on current happiness is
  52  *  (cur_hap * hap_avg + hap_prod * etu) / (hap_avg + etu);
  53  * same for education.
  54  * right now, happiness has 1 day (48 etu) average, prod of 10 from
  55  * initial level of 0 yields (1) 1.42, (6) 6.03, (12) 8.42, (18) 9.37
  56  *
  57  * education has 4 day (192 etu) average, prod of 10 from initial
  58  * level of 0 yields (1) 0.4, (6) 2.2, (12) 3.9, (18) 5.2.
  59  */
  60
  61 extern float    hap_avg;
  62 extern float    edu_avg;
  63 extern float    ally_factor;
  64
  65 /*
  66  * for values below the "easy level" values, production is
  67  * as normal.  For values above "easy", production gets harder
  68  * based on an equation in "limit_level()" in update/nat.c.
  69  * Basically, the smaller the the values for "level_log", the
  70  * smaller return on investment above level_easy[] values.
  71  */
  72 /*
  73  * Damn! I hate this, but ...
  74  * The values here for tech are *not* the real ones.
  75  * They are changed later in the limit_level routine.
  76  */
  77                         /*tech   res   edu   hap */
  78 float   level_easy[4] = { 0.75, 0.75, 5.00, 5.00 };
  79 float   level_log[4] =  { 1.75, 2.00, 4.00, 6.00 };
  80
  81 /*
  82  * technique to limit the sharpers who turn entire countries
  83  * into tech plants overnight...
  84  */
  85
  86 double
  87 logx(double d, double base)
  88 {
  89         if (base == 1.0)
  90                 return d;
  91         return log10(d) / log10(base);
  92 }
  93
  94 double
  95 limit_level(double level, int type, int flag)
  96 {
  97         double  above_easy;
  98         double  above;
  99         double  logbase;
 100         double  easy;
 101
 102 /*
 103  * Begin ugly hack.
 104  */
 105         extern float    easy_tech, tech_log_base;
 106
 107         level_easy[0] = easy_tech;
 108         level_log[0] = tech_log_base;
 109 /*
 110  * End ugly hack.
 111  */
 112
 113         if (level > level_easy[type]) {
 114                 logbase = level_log[type];
 115                 easy = level_easy[type];
 116                 above_easy = level - easy;
 117                 if (flag)
 118                         above = above_easy / logx(logbase + above_easy,logbase);
 119                 else
 120                         above = logx(above_easy + 1.0, logbase);
 121                 if (above > 250) above = 250;
 122                 return ((above) < 0) ? easy : (easy + above);
 123         } else
 124                 return level;
 125 }
 126
 127 void
 128 prod_nat(int etu)
 129 {
 130         extern  long money[MAXNOC];
 131         extern  long pops[MAXNOC];
 132         extern  double hap_cons, edu_cons;
 133         extern  long sea_money[MAXNOC];
 134         extern  long lnd_money[MAXNOC];
 135         extern  long air_money[MAXNOC];
 136         struct  natstr *np;
 137         float   hap;
 138         float   edu;
 139         float   hap_edu;
 140         long    pop;
 141         double  rlev;
 142         double  tlev;
 143         double  tech[MAXNOC];
 144         double  res[MAXNOC];
 145         double  newvalue;
 146         natid   n;
 147         int     cn,cont;
 148
 149         for (n=0; NULL != (np = getnatp(n)); n++) {
 150                 if ((np->nat_stat & STAT_NORM) == 0)
 151                         continue;
 152                 /*
 153                  * hap_edu: the more education people have, the
 154                  * more happiness they want.
 155                  */
 156                 hap_edu = np->nat_level[NAT_ELEV];
 157                 hap_edu = 1.5 - ((hap_edu + 10.0) / (hap_edu + 20.0));
 158                 pop = pops[n] + 1;
 159                 /*
 160                  * get per-population happiness and education
 161                  * see what the total per-civilian production is
 162                  * for this time period.
 163                  */
 164                 hap = levels[n][NAT_HLEV] * hap_edu * hap_cons /
 165                         ((float)pop * etu);
 166                 edu = levels[n][NAT_ELEV] * edu_cons /
 167                         ((float)pop * etu);
 168                 wu((natid)0, n, "%3.0f happiness, %3.0f education produced\n",
 169                         levels[n][NAT_HLEV], levels[n][NAT_ELEV]);
 170                 hap = limit_level(hap, NAT_HLEV, 1);
 171                 edu = limit_level(edu, NAT_ELEV, 1);
 172                 /*
 173                  * change the "moving average"...old happiness and
 174                  * education levels are weighted heavier than current
 175                  * production.
 176                  */
 177                 newvalue = (np->nat_level[NAT_HLEV] * hap_avg + hap * etu) /
 178                         (hap_avg + etu);
 179                 np->nat_level[NAT_HLEV] = newvalue;
 180                 newvalue = (np->nat_level[NAT_ELEV] * edu_avg + edu * etu) /
 181                         (edu_avg + etu);
 182                 np->nat_level[NAT_ELEV] = newvalue;
 183                 /*
 184                  * limit tech/research production
 185                  */
 186                 levels[n][NAT_TLEV] =
 187                         limit_level(levels[n][NAT_TLEV] / 1,
 188                                 NAT_TLEV, 0) * 1;
 189                 levels[n][NAT_RLEV] =
 190                         limit_level(levels[n][NAT_RLEV] / 1,
 191                                 NAT_RLEV, 0) * 1;
 192                 wu((natid)0, n,
 193                         "total pop is %d, yielding %4.2f hap, %4.2f edu\n",
 194                         pop - 1, hap, edu);
 195         }
 196         if (ally_factor > 0.0)
 197                 share_incr(res, tech);
 198         else {
 199                 bzero((s_char *)res, sizeof(res));
 200                 bzero((s_char *)tech, sizeof(tech));
 201         }
 202         for (n=0; NULL != (np = getnatp(n)); n++) {
 203                 if ((np->nat_stat & STAT_NORM) == 0)
 204                         continue;
 205                 tlev = levels[n][NAT_TLEV];
 206                 rlev = levels[n][NAT_RLEV];
 207                 if (tech[n] != 0.0 || res[n] != 0.0) {
 208                         wu((natid)0, n,
 209                                 "%5.4f technology (%5.4f + %5.4f), ",
 210                                 tlev + tech[n], tlev, tech[n]);
 211                         wu((natid)0, n,
 212                                 "%5.4f research (%5.4f + %5.4f) produced\n",
 213                                 rlev + res[n], rlev, res[n]);
 214                 } else
 215                         wu((natid)0, n,
 216                                 "%5.4f tech, %5.4f research produced\n",
 217                                      tlev, rlev);
 218                 rlev += res[n];
 219                 tlev += tech[n];
 220                 if (rlev != 0.0)
 221                         np->nat_level[NAT_RLEV] += rlev;
 222                 if (tlev != 0.0)
 223                         np->nat_level[NAT_TLEV] += tlev;
 224                 if ((sea_money[n] != 0) || (air_money[n] != 0) ||
 225                         (lnd_money[n] != 0))
 226                         wu((natid)0, n,
 227                         "Army delta $%d, Navy delta $%d, Air force delta $%d\n",
 228                                 lnd_money[n], sea_money[n], air_money[n]);
 229                 wu((natid)0, n, "money delta was $%d for this update\n",
 230                         np->nat_money - money[n]);
 231                 if (opt_LOSE_CONTACT) {
 232                     for (cn=0; cn <= MAXNOC; cn++) {
 233                         cont = getcontact(np, cn);
 234                         if (cont > 0) {
 235                                 logerror ("country %d at level %d with country %d.\n", n, cont, cn);
 236                                 setcont(n, cn, cont-1);
 237                         }
 238                     }
 239                 }
 240         }
 241 }
 242
 243 /*
 244  * find out everyones increment
 245  */
 246 void
 247 share_incr(register double *res, register double *tech)
 248 {
 249         register struct natstr *np;
 250         register struct natstr *other;
 251         register natid i;
 252         register natid j;
 253         int rnc;
 254         int tnc;
 255
 256         for (i=0; NULL != (np = getnatp(i)); i++) {
 257                 res[i] = tech[i] = 0.0;
 258                 if ((np->nat_stat & STAT_INUSE) == 0)
 259                         continue;
 260                 if (np->nat_stat & STAT_GOD)
 261                   continue;
 262                 if (np->nat_stat == VIS)
 263                   continue;
 264                 rnc = tnc = 0;
 265                 for (j=0; NULL != (other = getnatp(j)); j++) {
 266                         if (j == i)
 267                                 continue;
 268                         if (other->nat_stat & STAT_GOD)
 269                           continue;
 270                         if (other->nat_stat == VIS)
 271                           continue;
 272                         if ((other->nat_stat & STAT_INUSE) == 0)
 273                           continue;
 274                         if (opt_HIDDEN) {
 275                             if (!getcontact(np, j))
 276                                 continue;
 277                         }
 278                         if (!opt_ALL_BLEED) {
 279                           if (getrel(np, j) != ALLIED)
 280                             continue;
 281                           if (getrel(other, i) != ALLIED)
 282                             continue;
 283                           res[i] += levels[j][NAT_RLEV];
 284                           tech[i] += levels[j][NAT_TLEV];
 285                           rnc++;
 286                           tnc++;
 287                         } else {
 288                           if (levels[j][NAT_TLEV] > 0.001) {
 289                             tech[i] += levels[j][NAT_TLEV];
 290                             tnc++;
 291                           }
 292                           if (levels[j][NAT_RLEV] > 0.001) {
 293                             res[i] += levels[j][NAT_RLEV];
 294                             rnc++;
 295                           }
 296                         }
 297                 }
 298                 if (rnc == 0 && tnc == 0)
 299                   continue;
 300                 if (rnc > 0) {
 301                   res[i] /= rnc * ally_factor;
 302                 }
 303                 if (tnc > 0) {
 304                   tech[i] /= tnc * ally_factor;
 305                 }
 306 /*              logerror("Country #%d gets %g res from %d allies, %g tech from %d allies", i, res[i], rnc, tech[i], tnc);*/
 307         }
 308 }
 309