/*      $NetBSD: phaser.c,v 1.15 2009/08/12 08:54:54 dholland Exp $     */

/*
 * Copyright (c) 1980, 1993
 *      The Regents of the University of California.  All rights reserved.
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions
 * are met:
 * 1. Redistributions of source code must retain the above copyright
 *    notice, this list of conditions and the following disclaimer.
 * 2. Redistributions in binary form must reproduce the above copyright
 *    notice, this list of conditions and the following disclaimer in the
 *    documentation and/or other materials provided with the distribution.
 * 3. Neither the name of the University nor the names of its contributors
 *    may be used to endorse or promote products derived from this software
 *    without specific prior written permission.
 *
 * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
 * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
 * SUCH DAMAGE.
 */

#include <sys/cdefs.h>
#ifndef lint
#if 0
static char sccsid[] = "@(#)phaser.c    8.1 (Berkeley) 5/31/93";
#else
__RCSID("$NetBSD: phaser.c,v 1.15 2009/08/12 08:54:54 dholland Exp $");
#endif
#endif /* not lint */

#include <stdio.h>
#include <math.h>
#include "trek.h"
#include "getpar.h"

/* factors for phaser hits; see description below */

#define ALPHA           3.0             /* spread */
#define BETA            3.0             /* franf() */
#define GAMMA           0.30            /* cos(angle) */
#define EPSILON         150.0           /* dist ** 2 */
#define OMEGA           10.596          /* overall scaling factor */

/* OMEGA ~= 100 * (ALPHA + 1) * (BETA + 1) / (EPSILON + 1) */

/*
**  Phaser Control
**
**      There are up to NBANKS phaser banks which may be fired
**      simultaneously.  There are two modes, "manual" and
**      "automatic".  In manual mode, you specify exactly which
**      direction you want each bank to be aimed, the number
**      of units to fire, and the spread angle.  In automatic
**      mode, you give only the total number of units to fire.
**
**      The spread is specified as a number between zero and
**      one, with zero being minimum spread and one being maximum
**      spread.  You  will normally want zero spread, unless your
**      short range scanners are out, in which case you probably
**      don't know exactly where the Klingons are.  In that case,
**      you really don't have any choice except to specify a
**      fairly large spread.
**
**      Phasers spread slightly, even if you specify zero spread.
**
**      Uses trace flag 30
*/

static struct cvntab Matab[] = {
        { "m",          "anual",        (cmdfun) 1,     0 },
        { "a",          "utomatic",     (cmdfun) 0,     0 },
        { NULL,         NULL,           NULL,           0 }
};

struct banks {
        int     units;
        double  angle;
        double  spread;
};



/*ARGSUSED*/
void
phaser(int v __unused)
{
        int             i;
        int             j;
        struct kling    *k;
        double          dx, dy;
        double          anglefactor, distfactor;
        struct banks    *b;
        int             manual, flag, extra = 0;
        int             hit;
        double          tot;
        int             n;
        int             hitreqd[NBANKS];
        struct banks    bank[NBANKS];
        const struct cvntab     *ptr;

        if (Ship.cond == DOCKED) {
                printf("Phasers cannot fire through starbase shields\n");
                return;
        }
        if (damaged(PHASER)) {
                out(PHASER);
                return;
        }
        if (Ship.shldup) {
                printf("Sulu: Captain, we cannot fire through shields.\n");
                return;
        }
        if (Ship.cloaked) {
                printf("Sulu: Captain, surely you must realize that we cannot "
                       "fire\n");
                printf("  phasers with the cloaking device up.\n");
                return;
        }

        /* decide if we want manual or automatic mode */
        manual = 0;
        if (testnl()) {
                if (damaged(COMPUTER)) {
                        printf("%s", Device[COMPUTER].name);
                        manual++;
                } else if (damaged(SRSCAN)) {
                        printf("%s", Device[SRSCAN].name);
                        manual++;
                }
                if (manual)
                        printf(" damaged, manual mode selected\n");
        }

        if (!manual) {
                ptr = getcodpar("Manual or automatic", Matab);
                manual = (long) ptr->value;
        }
        if (!manual && damaged(COMPUTER)) {
                printf("Computer damaged, manual selected\n");
                skiptonl(0);
                manual++;
        }

        /* initialize the bank[] array */
        flag = 1;
        for (i = 0; i < NBANKS; i++)
                bank[i].units = 0;
        if (manual) {
                /* collect manual mode statistics */
                while (flag) {
                        printf("%d units available\n", Ship.energy);
                        extra = 0;
                        flag = 0;
                        for (i = 0; i < NBANKS; i++) {
                                b = &bank[i];
                                printf("\nBank %d:\n", i);
                                hit = getintpar("units");
                                if (hit < 0)
                                        return;
                                if (hit == 0)
                                        break;
                                extra += hit;
                                if (extra > Ship.energy) {
                                        printf("available energy exceeded.  ");
                                        skiptonl(0);
                                        flag++;
                                        break;
                                }
                                b->units = hit;
                                hit = getintpar("course");
                                if (hit < 0 || hit > 360)
                                        return;
                                b->angle = hit * 0.0174532925;
                                b->spread = getfltpar("spread");
                                if (b->spread < 0 || b->spread > 1)
                                        return;
                        }
                        Ship.energy -= extra;
                }
                extra = 0;
        } else {
                /* automatic distribution of power */
                if (Etc.nkling <= 0) {
                        printf("Sulu: But there are no Klingons in this "
                               "quadrant\n");
                        return;
                }
                printf("Phasers locked on target.  ");
                while (flag) {
                        printf("%d units available\n", Ship.energy);
                        hit = getintpar("Units to fire");
                        if (hit <= 0)
                                return;
                        if (hit > Ship.energy) {
                                printf("available energy exceeded.  ");
                                skiptonl(0);
                                continue;
                        }
                        flag = 0;
                        Ship.energy -= hit;
                        extra = hit;
                        n = Etc.nkling;
                        if (n > NBANKS)
                                n = NBANKS;
                        tot = n * (n + 1) / 2;
                        for (i = 0; i < n; i++) {
                                k = &Etc.klingon[i];
                                b = &bank[i];
                                distfactor = k->dist;
                                anglefactor = ALPHA * BETA * OMEGA /
                                        (distfactor * distfactor + EPSILON);
                                anglefactor *= GAMMA;
                                distfactor = k->power;
                                distfactor /= anglefactor;
                                hitreqd[i] = distfactor + 0.5;
                                dx = Ship.sectx - k->x;
                                dy = k->y - Ship.secty;
                                b->angle = atan2(dy, dx);
                                b->spread = 0.0;
                                b->units = ((n - i) / tot) * extra;
#ifdef xTRACE
                                if (Trace) {
                                        printf("b%d hr%d u%d df%.2f af%.2f\n",
                                                i, hitreqd[i], b->units,
                                                distfactor, anglefactor);
                                }
#endif
                                extra -= b->units;
                                hit = b->units - hitreqd[i];
                                if (hit > 0) {
                                        extra += hit;
                                        b->units -= hit;
                                }
                        }

                        /* give out any extra energy we might have around */
                        if (extra > 0) {
                                for (i = 0; i < n; i++) {
                                        b = &bank[i];
                                        hit = hitreqd[i] - b->units;
                                        if (hit <= 0)
                                                continue;
                                        if (hit >= extra) {
                                                b->units += extra;
                                                extra = 0;
                                                break;
                                        }
                                        b->units = hitreqd[i];
                                        extra -= hit;
                                }
                                if (extra > 0)
                                        printf("%d units overkill\n", extra);
                        }
                }
        }

#ifdef xTRACE
        if (Trace) {
                for (i = 0; i < NBANKS; i++) {
                        b = &bank[i];
                        printf("b%d u%d", i, b->units);
                        if (b->units > 0)
                                printf(" a%.2f s%.2f\n", b->angle, b->spread);
                        else
                                printf("\n");
                }
        }
#endif

        /* actually fire the shots */
        Move.free = 0;
        for (i = 0; i < NBANKS; i++) {
                b = &bank[i];
                if (b->units <= 0) {
                        continue;
                }
                printf("\nPhaser bank %d fires:\n", i);
                n = Etc.nkling;
                k = Etc.klingon;
                for (j = 0; j < n; j++) {
                        if (b->units <= 0)
                                break;
                        /*
                        ** The formula for hit is as follows:
                        **
                        **  zap = OMEGA * [(sigma + ALPHA) * (rho + BETA)]
                        **      / (dist ** 2 + EPSILON)]
                        **      * [cos(delta * sigma) + GAMMA]
                        **      * hit
                        **
                        ** where sigma is the spread factor,
                        ** rho is a random number (0 -> 1),
                        ** GAMMA is a crud factor for angle (essentially
                        **      cruds up the spread factor),
                        ** delta is the difference in radians between the
                        **      angle you are shooting at and the actual
                        **      angle of the klingon,
                        ** ALPHA scales down the significance of sigma,
                        ** BETA scales down the significance of rho,
                        ** OMEGA is the magic number which makes everything
                        **      up to "* hit" between zero and one,
                        ** dist is the distance to the klingon
                        ** hit is the number of units in the bank, and
                        ** zap is the amount of the actual hit.
                        **
                        ** Everything up through dist squared should maximize
                        ** at 1.0, so that the distance factor is never
                        ** greater than one.  Conveniently, cos() is
                        ** never greater than one, but the same restric-
                        ** tion applies.
                        */
                        distfactor = BETA + franf();
                        distfactor *= ALPHA + b->spread;
                        distfactor *= OMEGA;
                        anglefactor = k->dist;
                        distfactor /= anglefactor * anglefactor + EPSILON;
                        distfactor *= b->units;
                        dx = Ship.sectx - k->x;
                        dy = k->y - Ship.secty;
                        anglefactor = atan2(dy, dx) - b->angle;
                        anglefactor = cos((anglefactor * b->spread) + GAMMA);
                        if (anglefactor < 0.0) {
                                k++;
                                continue;
                        }
                        hit = anglefactor * distfactor + 0.5;
                        k->power -= hit;
                        printf("%d unit hit on Klingon", hit);
                        if (!damaged(SRSCAN))
                                printf(" at %d,%d", k->x, k->y);
                        printf("\n");
                        b->units -= hit;
                        if (k->power <= 0) {
                                killk(k->x, k->y);
                                continue;
                        }
                        k++;
                }
        }

        /* compute overkill */
        for (i = 0; i < NBANKS; i++)
                extra += bank[i].units;
        if (extra > 0)
                printf("\n%d units expended on empty space\n", extra);
}