#include "handy.h" #include "lut.h" /* lut.c this code (c)Copyright 1992 by Noah Spurrier */ /* These are the new Look Up Tables for Acid Warp. Many improvements are to be found. These tables are totally portable. Older versions worked only on 320x200 screen sizes. The tables could not handle larger screens, and the tables could not be enlarged because some would grow geometrically. These tables should handle screens up to about 2048x2048 (as currently hard-coded), but they can be easily expanded to handle larger screens as video display technoledgy advances. (uh huh) On the other hand, who knows when I might have access to a slide printer. The lut_dist() function is about 10% slower than the previous non-portable version, but it's still lots faster than trying to evaluate sqrt(x^2 + y^2). Also no math libraries need be linked using these funcions. */ #define SIN_TABLE_SIZE 1024 #define FRAME_SIZE 1023 int Sin_Table [SIN_TABLE_SIZE] = { 0, 3, 6, 9, 13, 16, 19, 22, 25, 28, 31, 34, 38, 41, 44, 47, 50, 53, 56, 59, 63, 66, 69, 72, 75, 78, 81, 84, 87, 91, 94, 97, 100, 103, 106, 109, 112, 115, 118, 121, 124, 127, 130, 133, 136, 139, 142, 145, 148, 151, 154, 157, 160, 163, 166, 169, 172, 175, 178, 181, 184, 187, 190, 193, 196, 199, 202, 204, 207, 210, 213, 216, 219, 222, 224, 227, 230, 233, 236, 238, 241, 244, 247, 249, 252, 255, 258, 260, 263, 266, 268, 271, 274, 276, 279, 282, 284, 287, 289, 292, 294, 297, 300, 302, 305, 307, 310, 312, 315, 317, 320, 322, 324, 327, 329, 332, 334, 336, 339, 341, 343, 346, 348, 350, 353, 355, 357, 359, 362, 364, 366, 368, 370, 373, 375, 377, 379, 381, 383, 385, 387, 389, 391, 393, 395, 397, 399, 401, 403, 405, 407, 409, 411, 413, 414, 416, 418, 420, 422, 423, 425, 427, 429, 430, 432, 434, 435, 437, 439, 440, 442, 443, 445, 446, 448, 449, 451, 452, 454, 455, 457, 458, 459, 461, 462, 464, 465, 466, 467, 469, 470, 471, 472, 474, 475, 476, 477, 478, 479, 480, 481, 482, 483, 484, 485, 486, 487, 488, 489, 490, 491, 492, 493, 493, 494, 495, 496, 497, 497, 498, 499, 499, 500, 501, 501, 502, 502, 503, 504, 504, 505, 505, 506, 506, 506, 507, 507, 508, 508, 508, 509, 509, 509, 509, 510, 510, 510, 510, 510, 511, 511, 511, 511, 511, 511, 511, 511, 511, 511, 511, 511, 511, 511, 510, 510, 510, 510, 510, 510, 509, 509, 509, 508, 508, 508, 507, 507, 507, 506, 506, 505, 505, 504, 504, 503, 503, 502, 502, 501, 500, 500, 499, 498, 498, 497, 496, 495, 495, 494, 493, 492, 491, 491, 490, 489, 488, 487, 486, 485, 484, 483, 482, 481, 480, 479, 478, 476, 475, 474, 473, 472, 471, 469, 468, 467, 465, 464, 463, 462, 460, 459, 457, 456, 455, 453, 452, 450, 449, 447, 446, 444, 443, 441, 439, 438, 436, 434, 433, 431, 429, 428, 426, 424, 423, 421, 419, 417, 415, 413, 412, 410, 408, 406, 404, 402, 400, 398, 396, 394, 392, 390, 388, 386, 384, 382, 380, 378, 376, 374, 371, 369, 367, 365, 363, 360, 358, 356, 354, 351, 349, 347, 345, 342, 340, 338, 335, 333, 330, 328, 326, 323, 321, 318, 316, 313, 311, 308, 306, 303, 301, 298, 296, 293, 291, 288, 285, 283, 280, 278, 275, 272, 270, 267, 264, 262, 259, 256, 253, 251, 248, 245, 242, 240, 237, 234, 231, 229, 226, 223, 220, 217, 214, 212, 209, 206, 203, 200, 197, 194, 191, 188, 186, 183, 180, 177, 174, 171, 168, 165, 162, 159, 156, 153, 150, 147, 144, 141, 138, 135, 132, 129, 126, 123, 120, 117, 114, 111, 107, 104, 101, 98, 95, 92, 89, 86, 83, 80, 77, 73, 70, 67, 64, 61, 58, 55, 52, 49, 45, 42, 39, 36, 33, 30, 27, 24, 20, 17, 14, 11, 8, 5, 2, -2, -5, -8, -11, -14, -17, -20, -24, -27, -30, -33, -36, -39, -42, -45, -49, -52, -55, -58, -61, -64, -67, -70, -73, -77, -80, -83, -86, -89, -92, -95, -98, -101, -104, -107, -111, -114, -117, -120, -123, -126, -129, -132, -135, -138, -141, -144, -147, -150, -153, -156, -159, -162, -165, -168, -171, -174, -177, -180, -183, -186, -188, -191, -194, -197, -200, -203, -206, -209, -212, -214, -217, -220, -223, -226, -229, -231, -234, -237, -240, -242, -245, -248, -251, -253, -256, -259, -262, -264, -267, -270, -272, -275, -278, -280, -283, -285, -288, -291, -293, -296, -298, -301, -303, -306, -308, -311, -313, -316, -318, -321, -323, -326, -328, -330, -333, -335, -338, -340, -342, -345, -347, -349, -351, -354, -356, -358, -360, -363, -365, -367, -369, -371, -374, -376, -378, -380, -382, -384, -386, -388, -390, -392, -394, -396, -398, -400, -402, -404, -406, -408, -410, -412, -413, -415, -417, -419, -421, -423, -424, -426, -428, -429, -431, -433, -434, -436, -438, -439, -441, -443, -444, -446, -447, -449, -450, -452, -453, -455, -456, -457, -459, -460, -462, -463, -464, -465, -467, -468, -469, -471, -472, -473, -474, -475, -476, -478, -479, -480, -481, -482, -483, -484, -485, -486, -487, -488, -489, -490, -491, -491, -492, -493, -494, -495, -495, -496, -497, -498, -498, -499, -500, -500, -501, -502, -502, -503, -503, -504, -504, -505, -505, -506, -506, -507, -507, -507, -508, -508, -508, -509, -509, -509, -510, -510, -510, -510, -510, -510, -511, -511, -511, -511, -511, -511, -511, -511, -511, -511, -511, -511, -511, -511, -510, -510, -510, -510, -510, -509, -509, -509, -509, -508, -508, -508, -507, -507, -506, -506, -506, -505, -505, -504, -504, -503, -502, -502, -501, -501, -500, -499, -499, -498, -497, -497, -496, -495, -494, -493, -493, -492, -491, -490, -489, -488, -487, -486, -485, -484, -483, -482, -481, -480, -479, -478, -477, -476, -475, -474, -472, -471, -470, -469, -467, -466, -465, -464, -462, -461, -459, -458, -457, -455, -454, -452, -451, -449, -448, -446, -445, -443, -442, -440, -439, -437, -435, -434, -432, -430, -429, -427, -425, -423, -422, -420, -418, -416, -414, -413, -411, -409, -407, -405, -403, -401, -399, -397, -395, -393, -391, -389, -387, -385, -383, -381, -379, -377, -375, -373, -370, -368, -366, -364, -362, -359, -357, -355, -353, -350, -348, -346, -343, -341, -339, -336, -334, -332, -329, -327, -324, -322, -320, -317, -315, -312, -310, -307, -305, -302, -300, -297, -294, -292, -289, -287, -284, -282, -279, -276, -274, -271, -268, -266, -263, -260, -258, -255, -252, -249, -247, -244, -241, -238, -236, -233, -230, -227, -224, -222, -219, -216, -213, -210, -207, -204, -202, -199, -196, -193, -190, -187, -184, -181, -178, -175, -172, -169, -166, -163, -160, -157, -154, -151, -148, -145, -142, -139, -136, -133, -130, -127, -124, -121, -118, -115, -112, -109, -106, -103, -100, -97, -94, -91, -87, -84, -81, -78, -75, -72, -69, -66, -63, -59, -56, -53, -50, -47, -44, -41, -38, -34, -31, -28, -25, -22, -19, -16, -13, -9, -6, -3, 0, }; long lut_sin (long angle) { if (angle < 0) angle = -angle + ANGLE_UNIT_HALF; angle %= ANGLE_UNIT; return (long)Sin_Table [angle * 4]; } /* Defined above as a macro in "lut.h" long lut_cos (long angle) { return lut_sin (angle + ANGLE_UNIT_QUART); } */ int Frame_Edge_Angle [FRAME_SIZE + 1] = { /* Experimental. Larger angle unit, looks smoother. Angle unit = 1024 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 18, 18, 18, 18, 18, 18, 19, 19, 19, 19, 19, 19, 19, 20, 20, 20, 20, 20, 20, 21, 21, 21, 21, 21, 21, 21, 22, 22, 22, 22, 22, 22, 23, 23, 23, 23, 23, 23, 24, 24, 24, 24, 24, 24, 24, 25, 25, 25, 25, 25, 25, 26, 26, 26, 26, 26, 26, 26, 27, 27, 27, 27, 27, 27, 28, 28, 28, 28, 28, 28, 28, 29, 29, 29, 29, 29, 29, 30, 30, 30, 30, 30, 30, 30, 31, 31, 31, 31, 31, 31, 32, 32, 32, 32, 32, 32, 32, 33, 33, 33, 33, 33, 33, 34, 34, 34, 34, 34, 34, 34, 35, 35, 35, 35, 35, 35, 36, 36, 36, 36, 36, 36, 36, 37, 37, 37, 37, 37, 37, 37, 38, 38, 38, 38, 38, 38, 39, 39, 39, 39, 39, 39, 39, 40, 40, 40, 40, 40, 40, 40, 41, 41, 41, 41, 41, 41, 42, 42, 42, 42, 42, 42, 42, 43, 43, 43, 43, 43, 43, 43, 44, 44, 44, 44, 44, 44, 44, 45, 45, 45, 45, 45, 45, 46, 46, 46, 46, 46, 46, 46, 47, 47, 47, 47, 47, 47, 47, 48, 48, 48, 48, 48, 48, 48, 49, 49, 49, 49, 49, 49, 49, 50, 50, 50, 50, 50, 50, 50, 51, 51, 51, 51, 51, 51, 51, 52, 52, 52, 52, 52, 52, 52, 53, 53, 53, 53, 53, 53, 53, 54, 54, 54, 54, 54, 54, 54, 55, 55, 55, 55, 55, 55, 55, 56, 56, 56, 56, 56, 56, 56, 57, 57, 57, 57, 57, 57, 57, 58, 58, 58, 58, 58, 58, 58, 59, 59, 59, 59, 59, 59, 59, 59, 60, 60, 60, 60, 60, 60, 60, 61, 61, 61, 61, 61, 61, 61, 62, 62, 62, 62, 62, 62, 62, 63, 63, 63, 63, 63, 63, 63, 63, 64, 64, 64, 64, 64, 64, 64, 65, 65, 65, 65, 65, 65, 65, 66, 66, 66, 66, 66, 66, 66, 66, 67, 67, 67, 67, 67, 67, 67, 68, 68, 68, 68, 68, 68, 68, 68, 69, 69, 69, 69, 69, 69, 69, 69, 70, 70, 70, 70, 70, 70, 70, 71, 71, 71, 71, 71, 71, 71, 71, 72, 72, 72, 72, 72, 72, 72, 72, 73, 73, 73, 73, 73, 73, 73, 74, 74, 74, 74, 74, 74, 74, 74, 75, 75, 75, 75, 75, 75, 75, 75, 76, 76, 76, 76, 76, 76, 76, 76, 77, 77, 77, 77, 77, 77, 77, 77, 78, 78, 78, 78, 78, 78, 78, 78, 79, 79, 79, 79, 79, 79, 79, 79, 80, 80, 80, 80, 80, 80, 80, 80, 81, 81, 81, 81, 81, 81, 81, 81, 82, 82, 82, 82, 82, 82, 82, 82, 83, 83, 83, 83, 83, 83, 83, 83, 83, 84, 84, 84, 84, 84, 84, 84, 84, 85, 85, 85, 85, 85, 85, 85, 85, 86, 86, 86, 86, 86, 86, 86, 86, 86, 87, 87, 87, 87, 87, 87, 87, 87, 88, 88, 88, 88, 88, 88, 88, 88, 88, 89, 89, 89, 89, 89, 89, 89, 89, 90, 90, 90, 90, 90, 90, 90, 90, 90, 91, 91, 91, 91, 91, 91, 91, 91, 91, 92, 92, 92, 92, 92, 92, 92, 92, 92, 93, 93, 93, 93, 93, 93, 93, 93, 93, 94, 94, 94, 94, 94, 94, 94, 94, 94, 95, 95, 95, 95, 95, 95, 95, 95, 95, 96, 96, 96, 96, 96, 96, 96, 96, 96, 97, 97, 97, 97, 97, 97, 97, 97, 97, 98, 98, 98, 98, 98, 98, 98, 98, 98, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 100, 100, 100, 100, 100, 100, 100, 100, 100, 101, 101, 101, 101, 101, 101, 101, 101, 101, 101, 102, 102, 102, 102, 102, 102, 102, 102, 102, 103, 103, 103, 103, 103, 103, 103, 103, 103, 103, 104, 104, 104, 104, 104, 104, 104, 104, 104, 104, 105, 105, 105, 105, 105, 105, 105, 105, 105, 105, 106, 106, 106, 106, 106, 106, 106, 106, 106, 106, 107, 107, 107, 107, 107, 107, 107, 107, 107, 107, 108, 108, 108, 108, 108, 108, 108, 108, 108, 108, 109, 109, 109, 109, 109, 109, 109, 109, 109, 109, 110, 110, 110, 110, 110, 110, 110, 110, 110, 110, 111, 111, 111, 111, 111, 111, 111, 111, 111, 111, 111, 112, 112, 112, 112, 112, 112, 112, 112, 112, 112, 113, 113, 113, 113, 113, 113, 113, 113, 113, 113, 113, 114, 114, 114, 114, 114, 114, 114, 114, 114, 114, 114, 115, 115, 115, 115, 115, 115, 115, 115, 115, 115, 115, 116, 116, 116, 116, 116, 116, 116, 116, 116, 116, 116, 117, 117, 117, 117, 117, 117, 117, 117, 117, 117, 117, 118, 118, 118, 118, 118, 118, 118, 118, 118, 118, 118, 119, 119, 119, 119, 119, 119, 119, 119, 119, 119, 119, 119, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 121, 121, 121, 121, 121, 121, 121, 121, 121, 121, 121, 121, 122, 122, 122, 122, 122, 122, 122, 122, 122, 122, 122, 122, 123, 123, 123, 123, 123, 123, 123, 123, 123, 123, 123, 124, 124, 124, 124, 124, 124, 124, 124, 124, 124, 124, 124, 125, 125, 125, 125, 125, 125, 125, 125, 125, 125, 125, 125, 125, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, */ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 23, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 28, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 29, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31 }; /* End Frame_Edge_Angle[] */ int Frame_Edge_Distance [FRAME_SIZE + 1] = { 1024,1024,1024,1024,1024,1024,1024,1024,1024,1024, 1024,1024,1024,1024,1024,1024,1024,1024,1024,1024, 1024,1024,1024,1024,1024,1024,1024,1024,1024,1024, 1024,1024,1024,1024,1024,1024,1024,1024,1024,1024, 1024,1024,1024,1024,1024,1024,1025,1025,1025,1025, 1025,1025,1025,1025,1025,1025,1025,1025,1025,1025, 1025,1025,1025,1025,1025,1026,1026,1026,1026,1026, 1026,1026,1026,1026,1026,1026,1026,1026,1026,1027, 1027,1027,1027,1027,1027,1027,1027,1027,1027,1027, 1027,1028,1028,1028,1028,1028,1028,1028,1028,1028, 1028,1028,1029,1029,1029,1029,1029,1029,1029,1029, 1029,1029,1030,1030,1030,1030,1030,1030,1030,1030, 1031,1031,1031,1031,1031,1031,1031,1031,1031,1032, 1032,1032,1032,1032,1032,1032,1032,1033,1033,1033, 1033,1033,1033,1033,1034,1034,1034,1034,1034,1034, 1034,1035,1035,1035,1035,1035,1035,1035,1036,1036, 1036,1036,1036,1036,1037,1037,1037,1037,1037,1037, 1038,1038,1038,1038,1038,1038,1039,1039,1039,1039, 1039,1039,1040,1040,1040,1040,1040,1040,1041,1041, 1041,1041,1041,1042,1042,1042,1042,1042,1042,1043, 1043,1043,1043,1043,1044,1044,1044,1044,1044,1045, 1045,1045,1045,1045,1046,1046,1046,1046,1046,1047, 1047,1047,1047,1048,1048,1048,1048,1048,1049,1049, 1049,1049,1049,1050,1050,1050,1050,1051,1051,1051, 1051,1051,1052,1052,1052,1052,1053,1053,1053,1053, 1054,1054,1054,1054,1055,1055,1055,1055,1056,1056, 1056,1056,1056,1057,1057,1057,1057,1058,1058,1058, 1058,1059,1059,1059,1060,1060,1060,1060,1061,1061, 1061,1061,1062,1062,1062,1062,1063,1063,1063,1064, 1064,1064,1064,1065,1065,1065,1065,1066,1066,1066, 1067,1067,1067,1067,1068,1068,1068,1069,1069,1069, 1069,1070,1070,1070,1071,1071,1071,1071,1072,1072, 1072,1073,1073,1073,1074,1074,1074,1074,1075,1075, 1075,1076,1076,1076,1077,1077,1077,1078,1078,1078, 1078,1079,1079,1079,1080,1080,1080,1081,1081,1081, 1082,1082,1082,1083,1083,1083,1084,1084,1084,1085, 1085,1085,1086,1086,1086,1087,1087,1087,1088,1088, 1088,1089,1089,1089,1090,1090,1090,1091,1091,1091, 1092,1092,1092,1093,1093,1093,1094,1094,1095,1095, 1095,1096,1096,1096,1097,1097,1097,1098,1098,1098, 1099,1099,1100,1100,1100,1101,1101,1101,1102,1102, 1103,1103,1103,1104,1104,1104,1105,1105,1106,1106, 1106,1107,1107,1107,1108,1108,1109,1109,1109,1110, 1110,1111,1111,1111,1112,1112,1112,1113,1113,1114, 1114,1114,1115,1115,1116,1116,1116,1117,1117,1118, 1118,1118,1119,1119,1120,1120,1120,1121,1121,1122, 1122,1122,1123,1123,1124,1124,1125,1125,1125,1126, 1126,1127,1127,1127,1128,1128,1129,1129,1130,1130, 1130,1131,1131,1132,1132,1133,1133,1133,1134,1134, 1135,1135,1136,1136,1136,1137,1137,1138,1138,1139, 1139,1139,1140,1140,1141,1141,1142,1142,1143,1143, 1143,1144,1144,1145,1145,1146,1146,1147,1147,1148, 1148,1148,1149,1149,1150,1150,1151,1151,1152,1152, 1153,1153,1153,1154,1154,1155,1155,1156,1156,1157, 1157,1158,1158,1159,1159,1160,1160,1160,1161,1161, 1162,1162,1163,1163,1164,1164,1165,1165,1166,1166, 1167,1167,1168,1168,1169,1169,1170,1170,1170,1171, 1171,1172,1172,1173,1173,1174,1174,1175,1175,1176, 1176,1177,1177,1178,1178,1179,1179,1180,1180,1181, 1181,1182,1182,1183,1183,1184,1184,1185,1185,1186, 1186,1187,1187,1188,1188,1189,1189,1190,1190,1191, 1191,1192,1192,1193,1193,1194,1195,1195,1196,1196, 1197,1197,1198,1198,1199,1199,1200,1200,1201,1201, 1202,1202,1203,1203,1204,1204,1205,1205,1206,1207, 1207,1208,1208,1209,1209,1210,1210,1211,1211,1212, 1212,1213,1213,1214,1215,1215,1216,1216,1217,1217, 1218,1218,1219,1219,1220,1220,1221,1222,1222,1223, 1223,1224,1224,1225,1225,1226,1227,1227,1228,1228, 1229,1229,1230,1230,1231,1231,1232,1233,1233,1234, 1234,1235,1235,1236,1237,1237,1238,1238,1239,1239, 1240,1240,1241,1242,1242,1243,1243,1244,1244,1245, 1246,1246,1247,1247,1248,1248,1249,1250,1250,1251, 1251,1252,1252,1253,1254,1254,1255,1255,1256,1256, 1257,1258,1258,1259,1259,1260,1261,1261,1262,1262, 1263,1263,1264,1265,1265,1266,1266,1267,1268,1268, 1269,1269,1270,1271,1271,1272,1272,1273,1274,1274, 1275,1275,1276,1277,1277,1278,1278,1279,1280,1280, 1281,1281,1282,1283,1283,1284,1284,1285,1286,1286, 1287,1287,1288,1289,1289,1290,1290,1291,1292,1292, 1293,1293,1294,1295,1295,1296,1296,1297,1298,1298, 1299,1300,1300,1301,1301,1302,1303,1303,1304,1305, 1305,1306,1306,1307,1308,1308,1309,1309,1310,1311, 1311,1312,1313,1313,1314,1314,1315,1316,1316,1317, 1318,1318,1319,1320,1320,1321,1321,1322,1323,1323, 1324,1325,1325,1326,1326,1327,1328,1328,1329,1330, 1330,1331,1332,1332,1333,1334,1334,1335,1335,1336, 1337,1337,1338,1339,1339,1340,1341,1341,1342,1343, 1343,1344,1344,1345,1346,1346,1347,1348,1348,1349, 1350,1350,1351,1352,1352,1353,1354,1354,1355,1356, 1356,1357,1358,1358,1359,1360,1360,1361,1361,1362, 1363,1363,1364,1365,1365,1366,1367,1367,1368,1369, 1369,1370,1371,1371,1372,1373,1373,1374,1375,1375, 1376,1377,1377,1378,1379,1379,1380,1381,1381,1382, 1383,1383,1384,1385,1385,1386,1387,1388,1388,1389, 1390,1390,1391,1392,1392,1393,1394,1394,1395,1396, 1396,1397,1398,1398,1399,1400,1400,1401,1402,1402, 1403,1404,1404,1405,1406,1407,1407,1408,1409,1409, 1410,1411,1411,1412,1413,1413,1414,1415,1416,1416, 1417,1418,1418,1419,1420,1420,1421,1422,1422,1423, 1424,1425,1425,1426,1427,1427,1428,1429,1429,1430, 1431,1431,1432,1433,1434,1434,1435,1436,1436,1437, 1438,1438,1439,1440,1441,1441,1442,1443,1443,1444, 1445,1446,1446,1447 }; /* End Frame_Edge_Distance[] */ long lut_angle (long x, long y) { int quadrant; long angle; long swap; /* Get and preserve quadrant info before we convert everything into quadrant one for processing. */ quadrant = 1; if ((x*y) > 0) { if (x < 0) { quadrant = 3; x = -x; y = -y; } } else { if (x < 0) { quadrant = 2; swap = -x; x = y; y = swap; } if (y < 0) { quadrant = 4; swap = -y; y = x; x = swap; } } if (x < y) { if (y) angle = ANGLE_UNIT_QUART - Frame_Edge_Angle[FRAME_SIZE * x / y] - 1; else angle = ANGLE_UNIT_QUART; } else { if (x) angle = Frame_Edge_Angle[FRAME_SIZE * y / x]; else angle = 0; } /* Now put the angle back in the proper quadrant */ switch (quadrant) { case 1: return angle; case 2: return angle + ANGLE_UNIT_QUART; case 3: return angle + ANGLE_UNIT_HALF ; case 4: return angle + ANGLE_UNIT_THREE_QUARTERS; } return -1; } /* End lut_angle() */ /* Finds the distance between two points; the (0,0) and (x,y) using a look up table. This is the new and improved version. Notice how is gleams. Find intersection with the edge of the virtual "frame" or virtual screen. Use look-up-table to determine the the distance of that point from the origin. Use similar triangles to scale that distance down to the actual distance of x,y from the origin. I.e. the ratio of the frame intersection distance to the size of the frame EQUALS the ratio of the x,y distance to y. Multiply by y to get x,y distance. _________. intersection | /| FrameIntrsctDist Dist(x,y) | x,y./ | ------------------ == ----------- | /| | FRAME_SIZE y | / | | |/____|___| The usual formula from high school is (sqrt(x^2 + y^2)). This is slower. Requires doubles. And requires linking math library. frame_intersect = FRAME_SIZE * x / y; find intersection frame_dist = Edge_Distance [frame_intersect]; find intersect distance in table return frame_dist * y / FRAME_SIZE; scale using ratios to get x,y distance */ long lut_dist (long x, long y) { x = ABS(x); y = ABS(y); /* Keep it in the first quadrant. */ if (y == 0) return x; /* Simple cases. Also avoid division by zero. */ if (x == 0) return y; /* Is the Intersection with top or with the side of the Frame? */ if (x < y) return Frame_Edge_Distance[FRAME_SIZE * x / y] * y / FRAME_SIZE; else return Frame_Edge_Distance[FRAME_SIZE * y / x] * x / FRAME_SIZE; } /* End lut_dist() */ #define SQUARE_TABLE_SIZE 200 long SQUARE_TABLE [SQUARE_TABLE_SIZE] = { 0, 1, 4, 9, 16, 25, 36, 49, 64, 81, 100, 121, 144, 169, 196, 225, 256, 289, 324, 361, 400, 441, 484, 529, 576, 625, 676, 729, 784, 841, 900, 961, 1024, 1089, 1156, 1225, 1296, 1369, 1444, 1521, 1600, 1681, 1764, 1849, 1936, 2025, 2116, 2209, 2304, 2401, 2500, 2601, 2704, 2809, 2916, 3025, 3136, 3249, 3364, 3481, 3600, 3721, 3844, 3969, 4096, 4225, 4356, 4489, 4624, 4761, 4900, 5041, 5184, 5329, 5476, 5625, 5776, 5929, 6084, 6241, 6400, 6561, 6724, 6889, 7056, 7225, 7396, 7569, 7744, 7921, 8100, 8281, 8464, 8649, 8836, 9025, 9216, 9409, 9604, 9801, 10000, 10201, 10404, 10609, 10816, 11025, 11236, 11449, 11664, 11881, 12100, 12321, 12544, 12769, 12996, 13225, 13456, 13689, 13924, 14161, 14400, 14641, 14884, 15129, 15376, 15625, 15876, 16129, 16384, 16641, 16900, 17161, 17424, 17689, 17956, 18225, 18496, 18769, 19044, 19321, 19600, 19881, 20164, 20449, 20736, 21025, 21316, 21609, 21904, 22201, 22500, 22801, 23104, 23409, 23716, 24025, 24336, 24649, 24964, 25281, 25600, 25921, 26244, 26569, 26896, 27225, 27556, 27889, 28224, 28561, 28900, 29241, 29584, 29929, 30276, 30625, 30976, 31329, 31684, 32041, 32400, 32761, 33124, 33489, 33856, 34225, 34596, 34969, 35344, 35721, 36100, 36481, 36864, 37249, 37636, 38025, 38416, 38809, 39204, 39601, }; /* Some developement test routines follow. They are precompiled out. */ #if 0 void main (void) { int x, y; printf ("\n * "); x = -100; for (y = -160; y < 161; ++y) printf ("%d ", lut_angle(y, x)); printf ("\n"); x = -160; for (y = -100; y < 101; ++y) printf ("%d ", lut_angle(x, y)); } #endif #if 0 void main (void) { long x; long p; long c; long z; START(); /* Benchmark macro from handy.h */ for (x = 0; x < 800; ++x) for (c = 0; c < 100; ++c) for (p = 0; p < 100; ++p) z = lut_dist_x (p, c); printf ("\nTime for old lut_dist: %ld\n", MARK()); START(); for (x = 0; x < 800; ++x) for (c = 0; c < 100; ++c) for (p = 0; p < 100; ++p) z = lut_dist (p, c); printf ("\nTime for new lut_dist: %ld\n", MARK()); } #endif #if 0 Fosil code may be brought back from extinction. This is how I used to implement the lut_dist() function. This is actually faster, but it's harder to understand, harder to rescale for larger screens, and table size grows geometrically with screen size. Can't get much more than 320x240. #define ATAN_TABLE_HI_SIZE 1000 #define ATAN_TABLE_LO_SIZE 960 #define SQRT_TABLE_HI_SIZE 2000 /*891*/ #define SQRT_TABLE_LO_SIZE 280 long SQRT_TABLE_LO [SQRT_TABLE_LO_SIZE] = { 0, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 17 }; int SQRT_TABLE_HI [SQRT_TABLE_HI_SIZE] = { 0, 6, 9, 11, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 24, 25, 26, 27, 28, 28, 29, 30, 30, 31, 32, 32, 33, 33, 34, 35, 35, 36, 36, 37, 37, 38, 38, 39, 39, 40, 40, 41, 41, 42, 42, 43, 43, 44, 44, 45, 45, 46, 46, 46, 47, 47, 48, 48, 49, 49, 49, 50, 50, 51, 51, 51, 52, 52, 53, 53, 53, 54, 54, 54, 55, 55, 55, 56, 56, 57, 57, 57, 58, 58, 58, 59, 59, 59, 60, 60, 60, 61, 61, 61, 62, 62, 62, 63, 63, 63, 64, 64, 64, 64, 65, 65, 65, 66, 66, 66, 67, 67, 67, 68, 68, 68, 68, 69, 69, 69, 70, 70, 70, 70, 71, 71, 71, 72, 72, 72, 72, 73, 73, 73, 73, 74, 74, 74, 75, 75, 75, 75, 76, 76, 76, 76, 77, 77, 77, 77, 78, 78, 78, 78, 79, 79, 79, 79, 80, 80, 80, 80, 81, 81, 81, 81, 82, 82, 82, 82, 83, 83, 83, 83, 84, 84, 84, 84, 85, 85, 85, 85, 86, 86, 86, 86, 86, 87, 87, 87, 87, 88, 88, 88, 88, 89, 89, 89, 89, 89, 90, 90, 90, 90, 91, 91, 91, 91, 91, 92, 92, 92, 92, 93, 93, 93, 93, 93, 94, 94, 94, 94, 94, 95, 95, 95, 95, 95, 96, 96, 96, 96, 97, 97, 97, 97, 97, 98, 98, 98, 98, 98, 99, 99, 99, 99, 99, 100, 100, 100, 100, 100, 101, 101, 101, 101, 101, 102, 102, 102, 102, 102, 103, 103, 103, 103, 103, 104, 104, 104, 104, 104, 104, 105, 105, 105, 105, 105, 106, 106, 106, 106, 106, 107, 107, 107, 107, 107, 108, 108, 108, 108, 108, 108, 109, 109, 109, 109, 109, 110, 110, 110, 110, 110, 110, 111, 111, 111, 111, 111, 112, 112, 112, 112, 112, 112, 113, 113, 113, 113, 113, 113, 114, 114, 114, 114, 114, 115, 115, 115, 115, 115, 115, 116, 116, 116, 116, 116, 116, 117, 117, 117, 117, 117, 117, 118, 118, 118, 118, 118, 118, 119, 119, 119, 119, 119, 119, 120, 120, 120, 120, 120, 120, 121, 121, 121, 121, 121, 121, 122, 122, 122, 122, 122, 122, 123, 123, 123, 123, 123, 123, 124, 124, 124, 124, 124, 124, 125, 125, 125, 125, 125, 125, 126, 126, 126, 126, 126, 126, 126, 127, 127, 127, 127, 127, 127, 128, 128, 128, 128, 128, 128, 129, 129, 129, 129, 129, 129, 129, 130, 130, 130, 130, 130, 130, 131, 131, 131, 131, 131, 131, 131, 132, 132, 132, 132, 132, 132, 133, 133, 133, 133, 133, 133, 133, 134, 134, 134, 134, 134, 134, 134, 135, 135, 135, 135, 135, 135, 135, 136, 136, 136, 136, 136, 136, 137, 137, 137, 137, 137, 137, 137, 138, 138, 138, 138, 138, 138, 138, 139, 139, 139, 139, 139, 139, 139, 140, 140, 140, 140, 140, 140, 140, 141, 141, 141, 141, 141, 141, 141, 142, 142, 142, 142, 142, 142, 142, 143, 143, 143, 143, 143, 143, 143, 144, 144, 144, 144, 144, 144, 144, 144, 145, 145, 145, 145, 145, 145, 145, 146, 146, 146, 146, 146, 146, 146, 147, 147, 147, 147, 147, 147, 147, 148, 148, 148, 148, 148, 148, 148, 148, 149, 149, 149, 149, 149, 149, 149, 150, 150, 150, 150, 150, 150, 150, 150, 151, 151, 151, 151, 151, 151, 151, 152, 152, 152, 152, 152, 152, 152, 152, 153, 153, 153, 153, 153, 153, 153, 153, 154, 154, 154, 154, 154, 154, 154, 155, 155, 155, 155, 155, 155, 155, 155, 156, 156, 156, 156, 156, 156, 156, 156, 157, 157, 157, 157, 157, 157, 157, 157, 158, 158, 158, 158, 158, 158, 158, 158, 159, 159, 159, 159, 159, 159, 159, 159, 160, 160, 160, 160, 160, 160, 160, 160, 161, 161, 161, 161, 161, 161, 161, 161, 162, 162, 162, 162, 162, 162, 162, 162, 163, 163, 163, 163, 163, 163, 163, 163, 164, 164, 164, 164, 164, 164, 164, 164, 165, 165, 165, 165, 165, 165, 165, 165, 166, 166, 166, 166, 166, 166, 166, 166, 166, 167, 167, 167, 167, 167, 167, 167, 167, 168, 168, 168, 168, 168, 168, 168, 168, 169, 169, 169, 169, 169, 169, 169, 169, 169, 170, 170, 170, 170, 170, 170, 170, 170, 171, 171, 171, 171, 171, 171, 171, 171, 171, 172, 172, 172, 172, 172, 172, 172, 172, 173, 173, 173, 173, 173, 173, 173, 173, 173, 174, 174, 174, 174, 174, 174, 174, 174, 174, 175, 175, 175, 175, 175, 175, 175, 175, 175, 176, 176, 176, 176, 176, 176, 176, 176, 177, 177, 177, 177, 177, 177, 177, 177, 177, 178, 178, 178, 178, 178, 178, 178, 178, 178, 179, 179, 179, 179, 179, 179, 179, 179, 179, 180, 180, 180, 180, 180, 180, 180, 180, 180, 181, 181, 181, 181, 181, 181, 181, 181, 181, 182, 182, 182, 182, 182, 182, 182, 182, 182, 183, 183, 183, 183, 183, 183, 183, 183, 183, 184, 184, 184, 184, 184, 184, 184, 184, 184, 184, 185, 185, 185, 185, 185, 185, 185, 185, 185, 186, 186, 186, 186, 186, 186, 186, 186, 186, 187, 187, 187, 187, 187, 187, 187, 187, 187, 188, 188, 188, 188, 188, 188, 188, 188, 188, 188, 189, 189, 189, 189, 189, 189, 189, 189, 189, 190, 190, 190, 190, 190, 190, 190, 190, 190, 190, 191, 191, 191, 191, 191, 191, 191, 191, 191, 192, 192, 192, 192, 192, 192, 192, 192, 192, 192, 193, 193, 193, 193, 193, 193, 193, 193, 193, 193, 194, 194, 194, 194, 194, 194, 194, 194, 194, 195, 195, 195, 195, 195, 195, 195, 195, 195, 195, 196, 196, 196, 196, 196, 196, 196, 196, 196, 196, 197, 197, 197, 197, 197, 197, 197, 197, 197, 197, 198, 198, 198, 198, 198, 198, 198, 198, 198, 198, 199, 199, 199, 199, 199, 199, 199, 199, 199, 199, 200, 200, 200, 200, 200, 200, 200, 200, 200, 200, 201, 201, 201, 201, 201, 201, 201, 201, 201, 201, 202, 202, 202, 202, 202, 202, 202, 202, 202, 202, 203, 203, 203, 203, 203, 203, 203, 203, 203, 203, 204, 204, 204, 204, 204, 204, 204, 204, 204, 204, 205, 205, 205, 205, 205, 205, 205, 205, 205, 205, 206, 206, 206, 206, 206, 206, 206, 206, 206, 206, 206, 207, 207, 207, 207, 207, 207, 207, 207, 207, 207, 208, 208, 208, 208, 208, 208, 208, 208, 208, 208, 209, 209, 209, 209, 209, 209, 209, 209, 209, 209, 209, 210, 210, 210, 210, 210, 210, 210, 210, 210, 210, 211, 211, 211, 211, 211, 211, 211, 211, 211, 211, 211, 212, 212, 212, 212, 212, 212, 212, 212, 212, 212, 213, 213, 213, 213, 213, 213, 213, 213, 213, 213, 213, 214, 214, 214, 214, 214, 214, 214, 214, 214, 214, 214, 215, 215, 215, 215, 215, 215, 215, 215, 215, 215, 215, 216, 216, 216, 216, 216, 216, 216, 216, 216, 216, 217, 217, 217, 217, 217, 217, 217, 217, 217, 217, 217, 218, 218, 218, 218, 218, 218, 218, 218, 218, 218, 218, 219, 219, 219, 219, 219, 219, 219, 219, 219, 219, 219, 220, 220, 220, 220, 220, 220, 220, 220, 220, 220, 220, 221, 221, 221, 221, 221, 221, 221, 221, 221, 221, 221, 222, 222, 222, 222, 222, 222, 222, 222, 222, 222, 222, 223, 223, 223, 223, 223, 223, 223, 223, 223, 223, 223, 224, 224, 224, 224, 224, 224, 224, 224, 224, 224, 224, 224, 225, 225, 225, 225, 225, 225, 225, 225, 225, 225, 225, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 226, 227, 227, 227, 227, 227, 227, 227, 227, 227, 227, 227, 228, 228, 228, 228, 228, 228, 228, 228, 228, 228, 228, 228, 229, 229, 229, 229, 229, 229, 229, 229, 229, 229, 229, 230, 230, 230, 230, 230, 230, 230, 230, 230, 230, 230, 230, 231, 231, 231, 231, 231, 231, 231, 231, 231, 231, 231, 232, 232, 232, 232, 232, 232, 232, 232, 232, 232, 232, 232, 233, 233, 233, 233, 233, 233, 233, 233, 233, 233, 233, 233, 234, 234, 234, 234, 234, 234, 234, 234, 234, 234, 234, 235, 235, 235, 235, 235, 235, 235, 235, 235, 235, 235, 235, 236, 236, 236, 236, 236, 236, 236, 236, 236, 236, 236, 236, 237, 237, 237, 237, 237, 237, 237, 237, 237, 237, 237, 237, 238, 238, 238, 238, 238, 238, 238, 238, 238, 238, 238, 238, 239, 239, 239, 239, 239, 239, 239, 239, 239, 239, 239, 239, 240, 240, 240, 240, 240, 240, 240, 240, 240, 240, 240, 240, 241, 241, 241, 241, 241, 241, 241, 241, 241, 241, 241, 241, 242, 242, 242, 242, 242, 242, 242, 242, 242, 242, 242, 242, 243, 243, 243, 243, 243, 243, 243, 243, 243, 243, 243, 243, 244, 244, 244, 244, 244, 244, 244, 244, 244, 244, 244, 244, 245, 245, 245, 245, 245, 245, 245, 245, 245, 245, 245, 245, 246, 246, 246, 246, 246, 246, 246, 246, 246, 246, 246, 246, 246, 247, 247, 247, 247, 247, 247, 247, 247, 247, 247, 247, 247, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 248, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 249, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 251, 251, 251, 251, 251, 251, 251, 251, 251, 251, 251, 251, 251, 252, 252, 252, 252, 252, 252, 252, 252, 252, 252, 252, 252, 253, 253, 253, 253, 253, 253, 253, 253, 253, 253, 253, 253, 253, 254, 254, 254, 254, 254, 254, 254, 254, 254, 254, 254, 254, 254, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 257, 257, 257, 257, 257, 257, 257, 257, 257, 257, 257, 257, 257, 258, 258, 258, 258, 258, 258, 258, 258, 258, 258, 258, 258, 258, 259, 259, 259, 259, 259, 259, 259, 259, 259, 259, 259, 259, 259, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 260, 261, 261, 261, 261, 261, 261, 261, 261, 261, 261, 261, 261, 261, 262, 262, 262, 262, 262, 262, 262, 262, 262, 262, 262, 262, 262, 263, 263, 263, 263, 263, 263, 263, 263, 263, 263, 263, 263, 263, 264, 264, 264, 264, 264, 264, 264, 264, 264, 264, 264, 264, 264, 264, 265, 265, 265, 265, 265, 265, 265, 265, 265, 265, 265, 265, 265, 266, 266, 266, 266, 266, 266, 266, 266, 266, 266, 266, 266, 266, 267, 267, 267, 267, 267, 267, 267, 267, 267, 267, 267, 267, 267, 268, 268, 268, 268, 268, 268, 268, 268, 268, 268, 268, 268, 268, 268, 269, 269, 269, 269, 269, 269, 269, 269, 269, 269, 269, 269, 269, 270, 270, 270, 270, 270, 270, 270, 270, 270, 270, 270, 270, 270, 270, 271, 271, 271, 271, 271, 271, 271, 271, 271, 271, 271, 271, 271, 272, 272, 272, 272, 272, 272, 272, 272, 272, 272, 272, 272, 272, 272, 273, 273, 273, 273, 273, 273, 273, 273, 273, 273, 273, 273, 273, 273, 274, 274, 274, 274, 274, 274, 274, 274, 274, 274, 274, 274, 274, 275, 275, 275, 275, 275, 275, 275, 275, 275, 275, 275, 275, 275, 275, 276, 276, 276, 276, 276, 276, 276, 276, 276, 276, 276, 276, 276, 276, 277, 277, 277, 277, 277, 277, 277, 277, 277, 277, 277, 277, 277, 277, 278, 278, 278, 278, 278, 278, 278, 278, 278, 278, 278, 278, 278, 278, 279, 279, 279, 279, 279, 279, 279, 279, 279, 279, 279, 279, 279, 279, 280, 280, 280, 280, 280, 280, 280, 280, 280, 280, 280, 280, 280, 280, 281, 281, 281, 281, 281, 281, 281, 281, 281, 281, 281, 281, 281, 281, 282, 282, 282, 282, 282, 282, 282, 282, 282, 282, 282, 282, 282, 282, 283, 283, 283, 283 }; /* Gets the distance between two points on the screen. * implements sqrt((x*x)+(y*y)), but faster... */ long lut_dist_x (long x, long y) { /* Some goofy ideas I had. I was trying to play with a linear series guestimate. But screw it, I didn't feel like doing real math. It would be slower, but would require NO tables. */ /* long ax, ay;*/ /* ax = ABS(x); ay = ABS(y);*/ /* return ((ax+ay)+(ax-ay))/2;*/ /* return ((ax+ay)+(2*ax+ay)+(ax+2*ay))/3;*/ /* return ((ax+ay)+ABS(ay-ax)+ABS(ax-ay))*3;*/ /* return ay + (ax + lut_sin (ax));*/ /* return (ax+ay)*(ax/(ay+1));*/ /* return (ABS(x)+ABS(y));*/ /* return ((ABS(x)+ABS(y))*140)/100;*/ long sqr; if (x < 0) x *= -1; if (y < 0) y *= -1; sqr = (long)SQUARE_TABLE[x] + (long)SQUARE_TABLE[y]; /* 40 is the magic number for table hashing. */ if (sqr >= SQRT_TABLE_LO_SIZE) return (SQRT_TABLE_HI [(int)(sqr/40)]); else return (SQRT_TABLE_LO[sqr]); } #endif #if 0 int VERT_ANGLE [VMAX] = { 0, 0, 0, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 7, 7, 7, 7, 8, 8, 8, 8, 9, 9, 9, 9, 10, 10, 10, 10, 10, 11, 11, 11, 11, 12, 12, 12, 12, 13, 13, 13, 13, 13, 14, 14, 14, 14, 15, 15, 15, 15, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 18, 18, 18, 18, 18, 19, 19, 19, 19, 19, 20, 20, 20, 20, 20, 21, 21, 21, 21, 21, 22, 22, 22, 22, 22, 23, 23 }; int HORIZ_ANGLE [HMAX] = { 63, 63, 63, 63, 62, 62, 62, 61, 61, 60, 60, 60, 59, 59, 58, 58, 58, 57, 57, 56, 56, 56, 55, 55, 54, 54, 54, 53, 53, 52, 52, 52, 51, 51, 51, 50, 50, 50, 49, 49, 48, 48, 48, 47, 47, 47, 46, 46, 46, 45, 45, 45, 44, 44, 44, 44, 43, 43, 43, 42, 42, 42, 41, 41, 41, 41, 40, 40, 40, 39, 39, 39, 39, 38, 38, 38, 38, 37, 37, 37, 37, 36, 36, 36, 36, 35, 35, 35, 35, 34, 34, 34, 34, 33, 33, 33, 33, 33, 32, 32, 32, 32, 32, 31, 31, 31, 31, 31, 30, 30, 30, 30, 30, 30, 29, 29, 29, 29, 29, 28, 28, 28, 28, 28, 28, 27, 27, 27, 27, 27, 27, 27, 26, 26, 26, 26, 26, 26, 26, 25, 25, 25, 25, 25, 25, 25, 24, 24, 24, 24, 24, 24, 24, 24, 23, 23, 23, 23, 23, 23, 23 }; /* You think you're ready? Why didn't I just use some arctan table? * Arctan would still require a double as an argument because arctan * takes the RATIO of y to x. Using the method I have created here the * division can be done last which enables me to make use of integer * division instead of having to convert numbers to double. * This method involves storing the angles at the edges of the screen. * For points inside the edges I use similar triangles to find the * appropriate angles. If I remember, this is about 10X faster. HA HA HA */ /* Get the angle between two points. Assume the angle starts from the * horizontal axis and moves counter-clockwise. */ long lut_angle (long dx, long dy) { int quadrant; long angle; long index; quadrant = 1; if (dx < 0) { quadrant = 2; dx *= -1; if (dy < 0) { quadrant = 3; dy *= -1; } } else if (dy < 0) { quadrant = 4; dy *= -1; } if (dy) /* Don't divide by zero */ index = 100 * dx / dy; else index = 999; /* 160 is half the size of the screen. Damn, I should use constants... */ if (index >= 160) { if (dx) index = 160 * dy / dx; else index = 0; angle = VERT_ANGLE[index]; } else angle = HORIZ_ANGLE [index]; switch (quadrant) { case 1: return angle; case 2: return ANGLE_UNIT_QUART - 1 + (ANGLE_UNIT_QUART - angle); case 3: return ANGLE_UNIT_HALF + angle; case 4: return ANGLE_UNIT_THREE_HALVES - 1+ (ANGLE_UNIT_QUART - angle); } return -ANGLE_UNIT; } #endif /* Many tables are hashed. Some tables are split into two parts to * save space. The works because some functions vary greatly at one * part of the curve, but they remain roughly linear at another part. * rather than have one large highly accurate table I used a small * highly accurate table for areas where the slope of a graph is steep * and I used a small low accuracy table for areas where the slope of * a graph is small. */ /* These coorespond to half of 320 and 200 plus one for fencepost * correction. The MCGA screen is 320x200 with 256 colors, so * many of these functions are optimised strictly for that screen size. * lut_angle() and lut_dist() are two that would require extensive * rewriting to make work on a non 320x200 screen */ #define HMAX 161 #define VMAX 101