home *** CD-ROM | disk | FTP | other *** search
/ Amiga Developer CD 2.1 / Amiga Developer CD v2.1.iso / NDK / NDK_1.3 / Read-Me1.3 / Workbench1.3 / Ieee.Examples / Ffuz.c next >
Encoding:
C/C++ Source or Header  |  1988-02-11  |  13.5 KB  |  593 lines
  1. /***********************************************************************
  2. *
  3. *            More Fuzzy Numbers on the Amiga
  4. *
  5. *       This Fuz shows a way to use direct calls to the IEEEDP math
  6. *    functions. This way people without the Manx mi.lib & mi32.lib
  7. *    files can compile this mess (just wait 'till you see the assembler
  8. *    version...) and get some fuzzy idea of how to use the MC68881 on
  9. *    Microbotics' Multi-Function Module.
  10. *       Fuz works off of three input variable - rot or rotation factor,
  11. *    vx and vy which are the X & Y vectors. The values that give the
  12. *    best results are floating point numbers of course. Use a value of
  13. *    .1 or less for rot, vy generally likes bigger #'s like 599.9893,
  14. *    and vx can be anything including 0 up to the 9 character limit on
  15. *    the input buffer. - when it's started it'll plot the default values
  16. *    I put in, and takes 20 seconds (roughly) to draw. I chose these
  17. *    values because they're a good example of the type of numerical
  18. *    input you should try. It's a toy. Play with it.
  19. *
  20. *        To compile FFuz.c with Manx, you need to link in two extra
  21. *    libraries, midb.lib and midt.lib that I hacked together for use with
  22. *     Manx. They stand for mathieeedoubbas.lib and mathieeedoubtrans.lib 
  23. *    respectively, which I've gotten damned tired of typing, so if you see
  24. *     midb or midt somewhere in my code, that's what they mean. I'm glad
  25. *    there's a 30 character filename size limit, 'cause gawd knows what
  26. *    weird library or function names that Dale Luck will come up with,
  27. *    you know? I mean, have you ever looked at the function names in the
  28. *    Layers library? I rest my case.
  29. *        Anyway, here's the Manx makefile for FFuz.c
  30.  
  31. CFLAGS    =    +fi +L -S
  32. OBJ    =    FFuz.o
  33.  
  34. ffuz:        $(OBJ)
  35.         ln -o FFuz $(OBJ) -lmidb -lmidt -lmi32 -lc32
  36.  
  37. *            author: David Milligan, for MicroBotics 
  38. *
  39. * FFuz source code is presented courtesy of MicroBotics, Inc., Richardson, TX.
  40. * FFuz is copyright 1987, MicroBotics, Inc. and David Milligan.
  41. * All rights reserved.
  42. ****************************************************************************/
  43.  
  44. #include <exec/types.h>
  45. #include <intuition/intuition.h>
  46. #include "devices/timer.h"
  47. #include <math.h>
  48. #include <stdio.h>
  49.  
  50. #define INTUITION    0x00000001
  51. #define GRAPHICS    0x00000002
  52. #define SCREEN        0x00000004
  53. #define WINDOW        0x00000008
  54. #define MIDB        0x00000010
  55. #define MIDT        0x00000020
  56. #define    NW_IDCMP    MOUSEBUTTONS|GADGETDOWN|GADGETUP|CLOSEWINDOW
  57. #define    NW_FLAGS    BORDERLESS|BACKDROP|ACTIVATE|WINDOWCLOSE
  58.  
  59. #ifndef TRUE
  60. #define TRUE 1
  61. #define FALSE 0
  62. #endif
  63.  
  64. #define CX    320
  65. #define CY    200
  66. #define XTWEAK    170
  67. #define YTWEAK    140
  68.  
  69. /*
  70.  * Gadget definitions
  71.  */
  72. #define ROT    1
  73. #define VX    2
  74. #define VY    3
  75. #define WHOA    4
  76. #define GO    5
  77. #define NUMGADS 5
  78. #define BOOLAFLAGS    RELVERIFY|GADGIMMEDIATE
  79. #define STRAFLAGS    RELVERIFY|GADGIMMEDIATE|STRINGCENTER
  80. #define GADL_EDGE    5
  81. #define GADT_EDGE    30
  82.  
  83. struct  GfxBase       *GfxBase;
  84. struct  IntuitionBase *IntuitionBase;
  85. struct  RastPort      *rp;
  86. struct  ViewPort      *vp;
  87. struct    Window    *w;
  88. struct    Screen    *screen;
  89. struct    IntuiMessage    *message;
  90.  
  91. UBYTE rotbuf[10] =    ".02";    /* This set of values draws a nice display */
  92. UBYTE vxbuf[10] =    "199";
  93. UBYTE vybuf[10] =    "192";
  94. char timebuf[25];
  95.  
  96. /*#define USE_COMPILER*/
  97. #ifdef USE_COMPILER
  98. #define    Add(x,y)    ((x)+(y))
  99. #define    Mul(x,y)    ((x)*(y))
  100. #define    Div(x,y)    ((x)/(y))
  101. #define    Sub(x,y)    ((x)-(y))
  102. #define Flt(x)    (x)
  103. #define Fix(x) (x)
  104. #define Sin sin
  105. #define Cos cos
  106. #define Abs fabs
  107. double    sin(),cos(),fabs();
  108. #else
  109. #define Add  IEEEDPAdd        /* We re-define the Double-Precision math */
  110. #define Mul IEEEDPMul        /* functions to a tolerable length.      */
  111. #define Div IEEEDPDiv
  112. double Add();        /* We want to make sure that we get a double value */
  113. double Mul();        /* returned from all the IEEEDP functions.....*/
  114. double Div();
  115. #define Flt IEEEDPFlt
  116. #define Fix IEEEDPFix
  117. double Flt();
  118.         /* Everything EXCEPT IEEDPFix() - this function is */
  119.         /* what turns the IEEEDP values into usable integer */
  120.         /* values. This hung me up for an embarassing amount */
  121.         /* of time whilst beating on the IEEEDP stuff. It must */
  122.         /* be declared as   int IEEEDPFix()   */
  123. int Fix();
  124. #define Sin IEEEDPSin
  125. #define Cos IEEEDPCos
  126. #define Abs IEEEDPAbs
  127. double Sin();
  128. double Cos();
  129. double Abs();
  130. #endif
  131.  
  132. #ifndef GREENHILLS
  133. long MathIeeeDoubBasBase;
  134. long MathIeeeDoubTransBase;
  135. #endif
  136.  
  137. double atof();
  138.  
  139. /*
  140.  *    color values
  141.  *        pen #      0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15
  142.  */
  143. USHORT red[] =        { 0, 13,11,12,13,13,14,15, 8, 4, 0, 0, 0, 0, 0, 4 };
  144. USHORT green[] =    { 0,  2, 3, 4, 5, 7, 9,13,11, 9, 6, 4, 2, 0, 5,10 };
  145. USHORT blue[] =        { 0,  2, 2, 2, 2, 1, 1, 0, 0, 0, 2, 7,10,14,15,15 };
  146.  
  147. struct TextAttr Topfont =
  148.     {
  149.     "topaz.font", TOPAZ_SIXTY,
  150.     FSF_EXTENDED | FSF_ITALIC | FSF_BOLD,
  151.     FPF_ROMFONT | FPF_PROPORTIONAL
  152.     };
  153. struct TextAttr Toyfont =
  154.     {
  155.     "topaz.font", TOPAZ_EIGHTY,
  156.     FSF_EXTENDED,
  157.     FPF_ROMFONT | FPF_PROPORTIONAL
  158.     };
  159. /*
  160.  *    border line vectors
  161.  */
  162. USHORT gadvecs2[] =    { 0,0,58,0,58,10,0,10,0,0 };
  163. USHORT gadvecs3[] =    { 0,0,84,0,84,12,0,12,0,0 };
  164. USHORT timervec[] =    {0,0,176,0,176,11,0,11,0,0};
  165.  
  166. struct IntuiText timer_text = { 7,0,JAM2,20,380, &Toyfont, &timebuf, NULL };
  167. struct Border timer_brdr = { 5,377, 1,0,JAM2, 5, timervec, NULL };
  168.  
  169. struct Border gad_brdr[NUMGADS] = 
  170. {
  171.     { -2,-3, 15,0,JAM1, 5, gadvecs3, NULL },
  172.     { -2,-3, 13,0,JAM1, 5, gadvecs3, NULL },
  173.     { -2,-3, 11,0,JAM1, 5, gadvecs3, NULL },
  174.     { -2,-1,  5,0,JAM1, 5, gadvecs2, NULL },
  175.     { -2,-1, 10,0,JAM1, 5, gadvecs2, NULL }
  176. };
  177.  
  178. struct IntuiText gad_txt[NUMGADS] =
  179. {
  180.     { 7,0,JAM2,    28,-12, &Toyfont, "rot", NULL },
  181.     { 7,0,JAM2,    34,-12, &Toyfont, "vx", NULL },
  182.     { 7,0,JAM2,    34,-12, &Toyfont, "vy", NULL },
  183.     { 7,4,JAM2,    1,   0, &Topfont, "Whoa!", NULL },
  184.     { 7,9,JAM2,    1,   0, &Topfont, "Start", NULL }
  185. };
  186.  
  187. struct StringInfo gad_info[NUMGADS-2] =
  188. { 
  189.     { rotbuf,  NULL, 0, 10, 0, 0,0,0,2,1, 0, NULL, NULL },
  190.     { vxbuf,   NULL, 0, 10, 0, 0,0,0,2,1, 0, NULL, NULL },
  191.     { vybuf,   NULL, 0, 10, 0, 0,0,0,2,1, 0, NULL, NULL }
  192. };
  193.  
  194. struct Gadget toy_gad[NUMGADS] = 
  195. {
  196.     {        /*    rotation string gadget    */
  197.     NULL,
  198.     GADL_EDGE,GADT_EDGE,
  199.     80,10,
  200.     GADGHCOMP,STRAFLAGS,STRGADGET,
  201.     (APTR)&gad_brdr[0],NULL, &gad_txt[0],
  202.     0, (APTR)&gad_info[0],
  203.     ROT, NULL
  204.     },
  205.     {        /*    x-vector input gadget    */
  206.     &toy_gad[0],
  207.     GADL_EDGE,GADT_EDGE+25,
  208.     80,10,
  209.     GADGHCOMP,STRAFLAGS,STRGADGET,
  210.     (APTR)&gad_brdr[1],NULL, &gad_txt[1],
  211.     0, (APTR)&gad_info[1],
  212.     VX, NULL
  213.     },
  214.     {        /*    y-vector input gadget    */
  215.     &toy_gad[1],
  216.     GADL_EDGE,GADT_EDGE+50,
  217.     80,10,
  218.     GADGHCOMP,STRAFLAGS,STRGADGET,
  219.     (APTR)&gad_brdr[2],NULL, &gad_txt[2],
  220.     0, (APTR)&gad_info[2],
  221.     VY, NULL
  222.     },
  223.     {        /*    stop button    */
  224.     &toy_gad[2],
  225.     GADL_EDGE+11,GADT_EDGE+65,
  226.     58,9,
  227.     GADGHCOMP, BOOLAFLAGS, BOOLGADGET,
  228.     (APTR)&gad_brdr[3],NULL, &gad_txt[3],
  229.     0, NULL,
  230.     WHOA, NULL
  231.     },
  232.     {        /*    start button     */
  233.     &toy_gad[3],
  234.     GADL_EDGE+11,GADT_EDGE+80,
  235.     58,9, 
  236.     GADGHCOMP, BOOLAFLAGS, BOOLGADGET, 
  237.     (APTR)&gad_brdr[4],NULL, &gad_txt[4],
  238.     0, NULL, 
  239.     GO, NULL
  240.     }
  241. };
  242.  
  243. struct NewScreen ns =
  244.     {
  245.     0, 0, 640, 400,
  246.     4, 0, 0,
  247.     HIRES|LACE, CUSTOMSCREEN,
  248.     &Toyfont,NULL, NULL
  249.     };
  250.  
  251. struct NewWindow nw =
  252.     {
  253.     0, 12, 640, 388, 7, 0,
  254.     NW_IDCMP,
  255.     NW_FLAGS,
  256.     &toy_gad[4], NULL,
  257.     "               MicroBotics' Fuzzy Number Plotter, by Dmil",
  258.     NULL, NULL, 0, 0, 0, 0, CUSTOMSCREEN
  259.     };
  260.  
  261. struct timeval t1;
  262. struct timeval t2;
  263.  
  264. extern struct MsgPort *CreatePort();
  265. extern struct IORequest *CreateExtIO();
  266.  
  267.  
  268. unsigned int mask = 0;
  269. BOOL idling;
  270.     /* A variable for everything I could think of, and a couple more */
  271.     /* just in case it wasn't confusing enough.    */
  272. int from_wb,col,cnt,nx,ny,n;
  273. double rot,vx,vy,inc,deginc,degvx,degvy,absvx,absvy,xtweak,ytweak;
  274. double PI,DEG,sintemp,sintempvx,costemp,costempvy,vnx,vny,n1,n2;
  275.  
  276. /************************************************
  277.  *         Main Program             *
  278.  ************************************************/
  279. main(argc)
  280. int argc;
  281. {
  282.     ULONG  class;
  283.     USHORT code;
  284.     register struct Gadget *gad;
  285.     register int gadgid;
  286.     
  287.     col = 1;
  288.     inc = nx = ny = cnt = 0;
  289.  
  290.     if(argc)        /* started from CLI or Workbench? */
  291.       from_wb = 0;
  292.     else
  293.       from_wb = -1;
  294.  
  295.     std_open();    /* open all the librarys, screens, etc. */
  296.  
  297.     xtweak = Flt(XTWEAK);
  298.     ytweak = Flt(YTWEAK);
  299.  
  300.             /* This could've be #defined , but I did that */
  301.             /* in the last program. Tempting fate... */
  302.     PI = 3.14159265358979323846;
  303.     n = 80;
  304.     n1 = Flt(n);
  305.     DEG = Div(PI,n1);
  306.  
  307.     rot = atof(rotbuf); /* use default values for initial plot */
  308.     vx = atof(vxbuf);
  309.     vy = atof(vybuf);
  310.     SetBPen(rp, 0);
  311.     SetAPen(rp,col);
  312.     Move(rp,CX,CY+YTWEAK);
  313.     idling = FALSE;
  314.  
  315.     GetSysTime(&t1);
  316.  
  317. /************************************************
  318.  *         Main Fuzzy Area            *
  319.  ************************************************/
  320.  
  321. FOREVER 
  322.   {
  323.     if((message = (struct IntuiMessage *)GetMsg(w->UserPort)) != NULL)
  324.     {
  325.     class = message->Class;
  326.     code  = message->Code;
  327.     gad = message->IAddress;
  328.     ReplyMsg(message);
  329.        switch (class)
  330.         {
  331.              case CLOSEWINDOW:
  332.         stage_left();        /* Anyone remember Snagglepuss? */
  333.         break;
  334.          case GADGETDOWN:
  335.          case GADGETUP:
  336.         gadgid = gad->GadgetID;
  337.         if(class == GADGETUP)
  338.            switch (gadgid)
  339.          {
  340.             case WHOA:
  341.                ny = YTWEAK;
  342.                GetSysTime(&t2);
  343.                showetime();
  344.                break;
  345.             case GO:
  346.                rot = atof(rotbuf);
  347.                vx = atof(vxbuf);
  348.                vy = atof(vybuf);
  349.                SetRast(rp,0);
  350.                RefreshWindowFrame(w);
  351.                col = 1;
  352.                SetAPen(rp,col);
  353.                Move(rp,CX,CY+YTWEAK);
  354.                inc =0;
  355.                cnt = 0;
  356.                ny = 0;
  357.                GetSysTime(&t1);
  358.                idling = FALSE;
  359.                break;
  360.  
  361.             default:
  362.               break;
  363.             }
  364.        default:
  365.         break;
  366.         }
  367.  
  368.     }
  369.     do_fuzzer();
  370.   }
  371.  stage_left();    /* main exit */
  372. }    /* end of main */
  373.  
  374.  
  375.     /*
  376.      * Here we are, all set to bang out fuzzy numbers & draw them in
  377.      * loud, garish colors. In interlace too, so that all you out 
  378.      * there without long persistance monitors (I've got one) get the
  379.      * added eyestrain of flickering teensy dots.
  380.      *  I've broken the math functions into individual calls stored
  381.      * in lots of variables to avoid extremely long lines of code.
  382.      * You may not like this approach. Bummer. I do.
  383.      */
  384. do_fuzzer()
  385. {
  386.        if(ny < YTWEAK)
  387.     {
  388.     inc = Add(inc,rot);
  389.  
  390.     deginc = Mul(inc,DEG);
  391.     sintemp = Sin(deginc);
  392.     costemp = Cos(deginc);
  393.     degvy = Mul(deginc,vy);
  394.     costempvy = Cos(degvy);
  395.     absvy = Abs(costempvy);
  396.     vny = Mul(costemp,absvy);
  397.     if(vx == 0)
  398.       {
  399.        nx = Fix(Mul(sintemp,xtweak));
  400.       } else {    
  401.      degvx = Mul(deginc,vx);
  402.      sintempvx = Sin(degvx);
  403.      absvx = Abs(sintempvx);
  404.      vnx = Mul(sintemp,absvx);
  405.     nx = Fix(Mul(vnx,xtweak));
  406.      }
  407.     ny = Fix(Mul(vny,ytweak));
  408.     cnt ++;
  409.     if(cnt >141)
  410.      {
  411.       cnt = 0;
  412.       col++;
  413.       if(col == 15) col = 1;
  414.       SetAPen(rp,col);
  415.      }
  416.     Draw(rp,nx+CX,ny+CY);
  417.          if(ny >= YTWEAK && !idling)
  418.         {
  419.         GetSysTime(&t2);
  420.          showetime();
  421.         }
  422.     }
  423. }
  424.  
  425. std_open()
  426. {
  427.     register int i;
  428.  
  429.     if(!(IntuitionBase = (struct IntuitionBase *)
  430.             OpenLibrary("intuition.library", 0)))
  431.         errexit(" Can't open Intuition");
  432.  
  433.     mask |= INTUITION;
  434.  
  435.     if(!(GfxBase = (struct GfxBase *)OpenLibrary("graphics.library", 0)))
  436.         errexit(" Can't open Gfx lib");
  437.  
  438.     mask |= GRAPHICS;
  439.  
  440. #ifndef GREENHILLS
  441.  
  442.     if(!(MathIeeeDoubBasBase = OpenLibrary("mathieeedoubbas.library",0)))
  443.         errexit(" no basic ieee library\n");
  444. #endif
  445.  
  446.     mask |= MIDB;
  447.  
  448. #ifndef GREENHILLS
  449.     if(!(MathIeeeDoubTransBase = OpenLibrary("mathieeedoubtrans.library",0)))
  450.         errexit(" no transcendental ieee library\n");
  451. #endif
  452.  
  453.     mask |= MIDT;
  454.  
  455.     if(!(screen = (struct Screen *)OpenScreen(&ns))) 
  456.         errexit(" Can't open the screen");
  457.  
  458.     mask |= SCREEN;
  459.  
  460.     nw.Screen = screen;
  461.  
  462.     ShowTitle(screen,FALSE);    /* turn off the screen title bar */
  463.     vp = &screen->ViewPort;
  464.         
  465. /*
  466.  * install some loud garish color
  467.  */
  468.     for(i = 0; i < 16; i++)
  469.         SetRGB4(vp,i,red[i],green[i],blue[i]);
  470.  
  471.     if(!(w = (struct Window *)OpenWindow(&nw)))
  472.         errexit(" Can't open the window");
  473.  
  474.     mask |= WINDOW;
  475.  
  476.     rp = w->RPort;
  477. }
  478.  
  479. errexit(str)
  480. register UBYTE *str;
  481. {
  482.     if(!from_wb)        /* started from CLI, so it's safe to */
  483.       printf("%s\n",str);    /* printf the error message    */
  484.  
  485.     stage_left();
  486. }
  487.  
  488. stage_left()
  489. {
  490. /*
  491.  *    use the masked flag bits to close what we opened
  492.  */
  493.  
  494.     if (mask & WINDOW)
  495.         CloseWindow(w);
  496.     if (mask & SCREEN)
  497.         CloseScreen(screen);
  498. #ifndef GREENHILLS
  499.     if (mask & MIDT)
  500.         CloseLibrary(MathIeeeDoubTransBase);
  501.     if (mask & MIDB)
  502.         CloseLibrary(MathIeeeDoubBasBase);
  503. #endif
  504.     if (mask & GRAPHICS)
  505.         CloseLibrary(GfxBase);
  506.     if (mask & INTUITION)
  507.         CloseLibrary(IntuitionBase);
  508.     exit();
  509. }
  510.  
  511. showetime()
  512. {
  513.     if(!idling)
  514.        sprintf(timebuf,"%ld.%ld seconds"
  515.         ,t2.tv_secs-t1.tv_secs,t1.tv_micro-t2.tv_micro);
  516.        DrawBorder(rp,&timer_brdr,0L,0L);
  517.        PrintIText(rp,&timer_text,0L,0L);
  518.        idling = TRUE;
  519. }
  520.  
  521. struct timerequest *PrepareTimer(precision)
  522. SHORT precision;
  523. {
  524.  
  525.     int error;
  526.     SHORT whichunit;
  527.  
  528.     struct MsgPort *timerport;
  529.     struct timerequest *timermsg;
  530.     
  531.         timerport = CreatePort(0,0);
  532.         if (timerport == NULL) 
  533.         return(NULL);    /* Error during CreatePort */
  534.  
  535.     timermsg = (struct timerequest *)CreateExtIO(
  536.                     timerport,sizeof(struct timerequest));
  537.     if (timermsg == NULL)
  538.         {
  539.         DeletePort(timerport);
  540.         return(NULL);    /* Error during CreateExtIO */
  541.         }
  542.     
  543.     if(precision)    /* if true, use precision timer  ( under 1 second ) */
  544.         whichunit = UNIT_MICROHZ;
  545.     else
  546.         whichunit = UNIT_VBLANK;
  547.  
  548.     error = OpenDevice(TIMERNAME, whichunit, timermsg, 0);
  549.     if (error != 0)
  550.         {
  551.         DeleteExtIO(timermsg,sizeof(struct timerequest));
  552.         DeletePort(timerport);
  553.         return(NULL);    /* Error during OpenDevice */
  554.         }
  555.     return(timermsg);
  556. }
  557.  
  558. int
  559. GetSysTime(tv)
  560. struct timeval *tv;
  561. {
  562.     struct timerequest *tr;
  563.     tr = PrepareTimer(TRUE);    /* MUST use Precise timer for this */
  564.     if(tr == 0) return(-1);        /* non zero return says error */
  565.  
  566.     tr->tr_node.io_Command = TR_GETSYSTIME;
  567.     DoIO( tr );
  568.  
  569.     tv->tv_secs = tr->tr_time.tv_secs;
  570.     tv->tv_micro = tr->tr_time.tv_micro;
  571.  
  572.     DeleteTimer(tr);
  573.     return(0);
  574.  
  575. }
  576.  
  577. int
  578. DeleteTimer(tr)
  579. struct timerequest *tr;
  580. {
  581.     struct MsgPort *tp;
  582.  
  583.     tp = tr->tr_node.io_Message.mn_ReplyPort;
  584.     if(tr != 0)
  585.     {
  586.         CloseDevice(tr);
  587.         DeleteExtIO(tr,sizeof(struct timerequest));
  588.     }
  589.     if(tp != 0)
  590.         DeletePort(tp);
  591.     return(0);        
  592. }
  593.