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
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.
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.
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
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.
28 * nat.c: Accumulate tech, edu, research and happiness.
30 * Known contributors to this file:
49 * hap and edu avg mean that the weight on current happiness is
50 * (cur_hap * hap_avg + hap_prod * etu) / (hap_avg + etu);
52 * right now, happiness has 1 day (48 etu) average, prod of 10 from
53 * initial level of 0 yields (1) 1.42, (6) 6.03, (12) 8.42, (18) 9.37
55 * education has 4 day (192 etu) average, prod of 10 from initial
56 * level of 0 yields (1) 0.4, (6) 2.2, (12) 3.9, (18) 5.2.
61 extern float ally_factor;
63 static void share_incr(register double *, register double *);
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.
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.
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 };
81 float levels[MAXNOC][4];
84 * technique to limit the sharpers who turn entire countries
85 * into tech plants overnight...
89 logx(double d, double base)
93 return log10(d) / log10(base);
97 limit_level(double level, int type, int flag)
107 level_easy[0] = easy_tech;
108 level_log[0] = tech_log_base;
113 if (level > level_easy[type]) {
114 logbase = level_log[type];
115 easy = level_easy[type];
116 above_easy = level - easy;
118 above = above_easy / logx(logbase + above_easy, logbase);
120 above = logx(above_easy + 1.0, logbase);
123 return ((above) < 0) ? easy : (easy + above);
131 extern long money[MAXNOC];
132 extern long pops[MAXNOC];
133 extern long sea_money[MAXNOC];
134 extern long lnd_money[MAXNOC];
135 extern long air_money[MAXNOC];
149 for (n = 0; NULL != (np = getnatp(n)); n++) {
150 if ((np->nat_stat & STAT_NORM) == 0)
153 * hap_edu: the more education people have, the
154 * more happiness they want.
156 hap_edu = np->nat_level[NAT_ELEV];
157 hap_edu = 1.5 - ((hap_edu + 10.0) / (hap_edu + 20.0));
160 * get per-population happiness and education
161 * see what the total per-civilian production is
162 * for this time period.
164 hap = levels[n][NAT_HLEV] * hap_edu * hap_cons /
166 edu = levels[n][NAT_ELEV] * edu_cons / ((float)pop * etu);
167 wu((natid)0, n, "%3.0f happiness, %3.0f education produced\n",
168 levels[n][NAT_HLEV], levels[n][NAT_ELEV]);
169 hap = limit_level(hap, NAT_HLEV, 1);
170 edu = limit_level(edu, NAT_ELEV, 1);
172 * change the "moving average"...old happiness and
173 * education levels are weighted heavier than current
176 newvalue = (np->nat_level[NAT_HLEV] * hap_avg + hap * etu) /
178 np->nat_level[NAT_HLEV] = newvalue;
179 newvalue = (np->nat_level[NAT_ELEV] * edu_avg + edu * etu) /
181 np->nat_level[NAT_ELEV] = newvalue;
183 * limit tech/research production
185 levels[n][NAT_TLEV] =
186 limit_level(levels[n][NAT_TLEV] / 1, NAT_TLEV, 0) * 1;
187 levels[n][NAT_RLEV] =
188 limit_level(levels[n][NAT_RLEV] / 1, NAT_RLEV, 0) * 1;
190 "total pop is %d, yielding %4.2f hap, %4.2f edu\n",
193 if (ally_factor > 0.0)
194 share_incr(res, tech);
196 memset(res, 0, sizeof(res));
197 memset(tech, 0, sizeof(tech));
199 for (n = 0; NULL != (np = getnatp(n)); n++) {
200 if ((np->nat_stat & STAT_NORM) == 0)
202 tlev = levels[n][NAT_TLEV];
203 rlev = levels[n][NAT_RLEV];
204 if (tech[n] != 0.0 || res[n] != 0.0) {
206 "%5.4f technology (%5.4f + %5.4f), ",
207 tlev + tech[n], tlev, tech[n]);
209 "%5.4f research (%5.4f + %5.4f) produced\n",
210 rlev + res[n], rlev, res[n]);
213 "%5.4f tech, %5.4f research produced\n", tlev, rlev);
217 np->nat_level[NAT_RLEV] += rlev;
219 np->nat_level[NAT_TLEV] += tlev;
220 if ((sea_money[n] != 0) || (air_money[n] != 0) ||
223 "Army delta $%d, Navy delta $%d, Air force delta $%d\n",
224 lnd_money[n], sea_money[n], air_money[n]);
225 wu((natid)0, n, "money delta was $%d for this update\n",
226 np->nat_money - money[n]);
227 if (opt_LOSE_CONTACT) {
228 for (cn = 0; cn <= MAXNOC; cn++) {
229 cont = getcontact(np, cn);
231 logerror("country %d at level %d with country %d.\n",
233 setcont(n, cn, cont - 1);
241 * find out everyones increment
244 share_incr(register double *res, register double *tech)
246 register struct natstr *np;
247 register struct natstr *other;
253 for (i = 0; NULL != (np = getnatp(i)); i++) {
254 res[i] = tech[i] = 0.0;
255 if ((np->nat_stat & STAT_INUSE) == 0)
257 if (np->nat_stat & STAT_GOD)
259 if (np->nat_stat == VIS)
262 for (j = 0; NULL != (other = getnatp(j)); j++) {
265 if (other->nat_stat & STAT_GOD)
267 if (other->nat_stat == VIS)
269 if ((other->nat_stat & STAT_INUSE) == 0)
272 if (!getcontact(np, j))
275 if (!opt_ALL_BLEED) {
276 if (getrel(np, j) != ALLIED)
278 if (getrel(other, i) != ALLIED)
280 res[i] += levels[j][NAT_RLEV];
281 tech[i] += levels[j][NAT_TLEV];
285 if (levels[j][NAT_TLEV] > 0.001) {
286 tech[i] += levels[j][NAT_TLEV];
289 if (levels[j][NAT_RLEV] > 0.001) {
290 res[i] += levels[j][NAT_RLEV];
295 if (rnc == 0 && tnc == 0)
298 res[i] /= rnc * ally_factor;
301 tech[i] /= tnc * ally_factor;
303 /* logerror("Country #%d gets %g res from %d allies, %g tech from %d allies", i, res[i], rnc, tech[i], tnc);*/