/ MacWorld 1997 September / Macworld (1997-09).dmg / Serious Software / Cherwell Scientific Demos / pro Fit / pro Fit 5.0 demo (68k).sea / pro Fit 5.0 demo (68k) / External Modules / External modules sources / C / Contour plotting / ContourPlotter.c

< prev

next >

view JSON data

view as text

open on a Mac

open on a PC

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 00 10 43 6f 6e 74 6f 75 | 72 50 6c 6f 74 74 65 72 |..Contou|rPlotter|
|00000010| 2e 63 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.c......|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 54 45 58 54 43 57 49 | 45 01 00 00 00 00 00 00 |.TEXTCWI|E.......|
|00000050| 00 00 00 00 00 2e 31 00 | 00 01 9a ad a0 43 cf ad |......1.|.....C..|
|00000060| a0 43 cf 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.C......|........|
|00000070| 00 00 00 00 00 00 00 00 | 00 00 81 81 60 a8 00 00 |........|....`...|
|00000080| 23 69 6e 63 6c 75 64 65 | 20 22 70 72 6f 46 69 74 |#include| "proFit|
|00000090| 5f 69 6e 74 65 72 66 61 | 63 65 2e 68 22 0d 0d 23 |_interfa|ce.h"..#|
|000000a0| 69 66 6e 64 65 66 20 5f | 5f 43 4f 4e 54 4f 55 52 |ifndef _|_CONTOUR|
|000000b0| 5f 50 4c 4f 54 54 45 52 | 5f 5f 0d 23 69 6e 63 6c |_PLOTTER|__.#incl|
|000000c0| 75 64 65 20 22 43 6f 6e | 74 6f 75 72 50 6c 6f 74 |ude "Con|tourPlot|
|000000d0| 74 65 72 2e 68 22 0d 23 | 65 6e 64 69 66 0d 0d 0d |ter.h".#|endif...|
|000000e0| 74 79 70 65 64 65 66 20 | 73 74 72 75 63 74 20 4d |typedef |struct M|
|000000f0| 79 50 72 69 76 61 74 65 | 50 6c 6f 74 74 69 6e 67 |yPrivate|Plotting|
|00000100| 44 61 74 61 0d 7b 0d 09 | 6c 6f 6e 67 09 09 09 6e |Data.{..|long...n|
|00000110| 72 52 6f 77 73 2c 20 6e | 72 43 6f 6c 73 3b 09 09 |rRows, n|rCols;..|
|00000120| 09 2f 2f 20 74 68 65 20 | 6e 75 6d 62 65 72 20 6f |.// the |number o|
|00000130| 66 20 72 6f 77 73 20 61 | 6e 64 20 63 6f 6c 75 6d |f rows a|nd colum|
|00000140| 6e 73 20 74 6f 20 73 63 | 61 6e 0d 09 76 6f 69 64 |ns to sc|an..void|
|00000150| 2a 2a 09 09 09 6d 65 73 | 68 3b 09 09 09 09 09 2f |**...mes|h;...../|
|00000160| 2f 20 70 6f 69 6e 74 65 | 72 20 74 6f 20 61 20 32 |/ pointe|r to a 2|
|00000170| 2d 64 20 61 72 72 61 79 | 20 6f 66 20 6d 65 73 68 |-d array| of mesh|
|00000180| 20 70 6f 69 6e 74 73 2c | 20 73 65 65 20 62 65 6c | points,| see bel|
|00000190| 6f 77 0d 09 42 6f 6f 6c | 65 61 6e 09 09 09 65 78 |ow..Bool|ean...ex|
|000001a0| 69 74 3b 09 09 09 09 09 | 2f 2f 20 73 65 74 20 74 |it;.....|// set t|
|000001b0| 6f 20 74 72 75 65 20 69 | 66 20 77 65 20 6d 75 73 |o true i|f we mus|
|000001c0| 74 20 65 78 69 74 0d 09 | 64 6f 75 62 6c 65 09 09 |t exit..|double..|
|000001d0| 09 63 4d 69 6e 2c 20 63 | 4d 61 78 3b 09 09 09 09 |.cMin, c|Max;....|
|000001e0| 2f 2f 20 76 61 6c 75 65 | 73 20 63 6f 72 72 65 73 |// value|s corres|
|000001f0| 70 6f 6e 64 69 6e 67 20 | 74 6f 20 66 69 72 73 74 |ponding |to first|
|00000200| 20 61 6e 64 20 6c 61 73 | 74 20 63 6f 6c 75 6d 6e | and las|t column|
|00000210| 0d 09 64 6f 75 62 6c 65 | 09 09 09 72 4d 69 6e 2c |..double|...rMin,|
|00000220| 20 72 4d 61 78 3b 09 09 | 09 09 2f 2f 20 76 61 6c | rMax;..|..// val|
|00000230| 75 65 73 20 63 6f 72 72 | 65 73 70 6f 6e 64 69 6e |ues corr|espondin|
|00000240| 67 20 74 6f 20 66 69 72 | 73 74 20 61 6e 64 20 6c |g to fir|st and l|
|00000250| 61 73 74 20 72 6f 77 0d | 09 64 6f 75 62 6c 65 09 |ast row.|.double.|
|00000260| 09 09 72 55 6e 69 74 2c | 20 63 55 6e 69 74 3b 09 |..rUnit,| cUnit;.|
|00000270| 09 09 2f 2f 20 28 78 69 | 4d 61 78 2d 78 69 4d 69 |..// (xi|Max-xiMi|
|00000280| 6e 29 2f 28 6e 72 52 6f | 77 73 28 6e 72 43 6f 6c |n)/(nrRo|ws(nrCol|
|00000290| 73 29 2d 31 29 0d 0d 7d | 4d 79 50 6c 6f 74 74 69 |s)-1)..}|MyPlotti|
|000002a0| 6e 67 44 61 74 61 3b 0d | 0d 0d 2f 2f 20 74 68 65 |ngData;.|..// the|
|000002b0| 20 66 6f 6c 6c 6f 77 69 | 6e 67 20 73 74 72 75 63 | followi|ng struc|
|000002c0| 74 20 67 69 76 65 73 20 | 69 6e 66 6f 20 61 62 6f |t gives |info abo|
|000002d0| 75 74 20 61 20 73 69 6e | 67 6c 65 20 6d 65 73 68 |ut a sin|gle mesh|
|000002e0| 20 70 6f 69 6e 74 0d 74 | 79 70 65 64 65 66 20 73 | point.t|ypedef s|
|000002f0| 74 72 75 63 74 0d 7b 0d | 09 64 6f 75 62 6c 65 09 |truct.{.|.double.|
|00000300| 76 61 6c 75 65 3b 09 09 | 09 09 09 2f 2a 20 74 68 |value;..|.../* th|
|00000310| 65 20 76 61 6c 75 65 20 | 61 74 20 74 68 65 20 67 |e value |at the g|
|00000320| 69 76 65 6e 20 70 6f 69 | 6e 74 20 2a 2f 0d 09 42 |iven poi|nt */..B|
|00000330| 6f 6f 6c 65 61 6e 09 76 | 61 6c 69 64 3b 09 09 09 |oolean.v|alid;...|
|00000340| 09 09 2f 2a 20 74 72 75 | 65 20 69 66 20 74 68 69 |../* tru|e if thi|
|00000350| 73 20 70 6f 69 6e 74 20 | 68 61 73 20 76 61 6c 69 |s point |has vali|
|00000360| 64 20 64 61 74 61 20 2a | 2f 0d 09 42 6f 6f 6c 65 |d data *|/..Boole|
|00000370| 61 6e 09 61 74 52 69 67 | 68 74 3b 09 09 09 09 2f |an.atRig|ht;..../|
|00000380| 2a 20 74 72 75 65 20 69 | 66 20 61 20 63 6f 6e 74 |* true i|f a cont|
|00000390| 6f 75 72 20 70 61 73 73 | 65 73 20 74 68 65 20 6d |our pass|es the m|
|000003a0| 65 73 68 20 74 6f 20 74 | 68 65 20 72 69 67 68 74 |esh to t|he right|
|000003b0| 20 6f 66 20 74 68 65 20 | 70 6f 69 6e 74 20 2a 2f | of the |point */|
|000003c0| 0d 09 42 6f 6f 6c 65 61 | 6e 09 61 74 42 6f 74 74 |..Boolea|n.atBott|
|000003d0| 6f 6d 3b 09 09 09 09 2f | 2a 20 74 72 75 65 20 69 |om;..../|* true i|
|000003e0| 66 20 61 20 63 6f 6e 74 | 6f 75 72 20 70 61 73 73 |f a cont|our pass|
|000003f0| 65 73 20 74 68 65 20 6d | 65 73 73 20 74 6f 20 74 |es the m|ess to t|
|00000400| 68 65 20 62 6f 74 74 6f | 6d 20 6f 66 20 74 68 65 |he botto|m of the|
|00000410| 20 70 6f 69 6e 74 20 2a | 2f 0d 7d 20 4d 65 73 68 | point *|/.} Mesh|
|00000420| 50 6f 69 6e 74 3b 0d 0d | 0d 0d 0d 2f 2a 20 74 68 |Point;..|.../* th|
|00000430| 65 20 66 6f 6c 6c 6f 77 | 69 6e 67 20 69 73 20 74 |e follow|ing is t|
|00000440| 68 65 20 73 69 7a 65 20 | 6f 66 20 61 6e 20 65 6e |he size |of an en|
|00000450| 74 72 79 20 69 6e 20 74 | 68 65 20 6d 65 73 68 20 |try in t|he mesh |
|00000460| 61 72 72 61 79 20 2d 20 | 6e 6f 74 65 20 74 68 61 |array - |note tha|
|00000470| 74 20 77 65 20 6d 61 6b | 65 20 73 75 72 65 20 2a |t we mak|e sure *|
|00000480| 2f 0d 2f 2a 20 74 68 61 | 74 20 74 68 65 20 73 69 |/./* tha|t the si|
|00000490| 7a 65 20 69 73 20 61 20 | 6d 75 6c 74 69 70 6c 65 |ze is a |multiple|
|000004a0| 20 6f 66 20 38 20 66 6f | 72 20 69 6e 63 72 65 61 | of 8 fo|r increa|
|000004b0| 73 69 6e 67 20 73 70 65 | 65 64 20 6f 6e 20 73 79 |sing spe|ed on sy|
|000004c0| 73 74 65 6d 73 20 77 68 | 65 72 65 20 64 6f 75 62 |stems wh|ere doub|
|000004d0| 6c 65 73 20 73 68 6f 75 | 6c 64 20 6e 6f 74 20 63 |les shou|ld not c|
|000004e0| 72 6f 73 73 20 2a 2f 0d | 2f 2a 20 32 2d 20 6f 72 |ross */.|/* 2- or|
|000004f0| 20 34 20 62 79 74 65 20 | 62 6f 75 6e 64 61 72 69 | 4 byte |boundari|
|00000500| 65 73 20 2a 2f 0d 23 64 | 65 66 69 6e 65 20 6d 65 |es */.#d|efine me|
|00000510| 73 68 50 6f 69 6e 74 53 | 69 7a 65 20 28 28 6c 6f |shPointS|ize ((lo|
|00000520| 6e 67 29 28 28 73 69 7a | 65 6f 66 28 4d 65 73 68 |ng)((siz|eof(Mesh|
|00000530| 50 6f 69 6e 74 29 2b 37 | 29 26 30 78 46 46 46 46 |Point)+7|)&0xFFFF|
|00000540| 46 46 38 29 29 0d 0d 23 | 64 65 66 69 6e 65 20 41 |FF8))..#|define A|
|00000550| 42 53 28 78 29 20 20 28 | 28 78 29 3c 30 3f 20 28 |BS(x)  (|(x)<0? (|
|00000560| 2d 28 78 29 29 3a 20 28 | 78 29 29 0d 0d 0d 73 74 |-(x)): (|x))...st|
|00000570| 61 74 69 63 20 4d 65 73 | 68 50 6f 69 6e 74 2a 20 |atic Mes|hPoint* |
|00000580| 47 65 74 4d 65 73 68 50 | 6f 69 6e 74 20 28 6c 6f |GetMeshP|oint (lo|
|00000590| 6e 67 20 72 6f 77 2c 20 | 6c 6f 6e 67 20 63 6f 6c |ng row, |long col|
|000005a0| 75 6d 6e 2c 20 4d 79 50 | 6c 6f 74 74 69 6e 67 44 |umn, MyP|lottingD|
|000005b0| 61 74 61 2a 20 63 6f 6e | 73 74 20 76 29 0d 09 2f |ata* con|st v)../|
|000005c0| 2f 20 72 65 74 75 72 6e | 73 20 74 68 65 20 61 64 |/ return|s the ad|
|000005d0| 64 72 65 73 73 20 6f 66 | 20 74 68 65 20 6d 65 73 |dress of| the mes|
|000005e0| 68 20 70 6f 69 6e 74 20 | 61 74 20 74 68 65 20 67 |h point |at the g|
|000005f0| 69 76 65 6e 20 72 6f 77 | 2f 63 6f 6c 75 6d 6e 0d |iven row|/column.|
|00000600| 09 2f 2f 20 72 6f 77 3a | 20 30 2e 2e 76 2d 3e 6e |.// row:| 0..v->n|
|00000610| 72 52 6f 77 73 2d 31 2c | 20 63 6f 6c 3a 20 30 2e |rRows-1,| col: 0.|
|00000620| 2e 76 2d 3e 6e 72 43 6f | 6c 73 2d 31 0d 7b 0d 09 |.v->nrCo|ls-1.{..|
|00000630| 6c 6f 6e 67 09 6f 66 66 | 73 65 74 20 3d 20 6d 65 |long.off|set = me|
|00000640| 73 68 50 6f 69 6e 74 53 | 69 7a 65 20 2a 20 28 72 |shPointS|ize * (r|
|00000650| 6f 77 2a 76 2d 3e 6e 72 | 43 6f 6c 73 20 2b 20 63 |ow*v->nr|Cols + c|
|00000660| 6f 6c 75 6d 6e 29 3b 0d | 0d 09 72 65 74 75 72 6e |olumn);.|..return|
|00000670| 20 28 4d 65 73 68 50 6f | 69 6e 74 2a 29 28 2a 28 | (MeshPo|int*)(*(|
|00000680| 63 68 61 72 2a 2a 29 76 | 2d 3e 6d 65 73 68 20 2b |char**)v|->mesh +|
|00000690| 20 6f 66 66 73 65 74 29 | 3b 0d 7d 0d 0d 0d 73 74 | offset)|;.}...st|
|000006a0| 61 74 69 63 20 42 6f 6f | 6c 65 61 6e 20 54 65 73 |atic Boo|lean Tes|
|000006b0| 74 49 66 45 78 69 74 28 | 4d 79 50 6c 6f 74 74 69 |tIfExit(|MyPlotti|
|000006c0| 6e 67 44 61 74 61 2a 20 | 63 6f 6e 73 74 20 76 29 |ngData* |const v)|
|000006d0| 0d 09 2f 2a 20 72 65 74 | 75 72 6e 73 20 74 72 75 |../* ret|urns tru|
|000006e0| 65 20 69 66 20 76 2d 3e | 65 78 69 74 20 69 73 20 |e if v->|exit is |
|000006f0| 74 72 75 65 20 6f 72 20 | 69 66 20 75 73 65 72 20 |true or |if user |
|00000700| 70 72 65 73 73 65 64 20 | 63 6d 64 2d 2e 20 2a 2f |pressed |cmd-. */|
|00000710| 0d 7b 0d 09 69 66 20 28 | 54 65 73 74 53 74 6f 70 |.{..if (|TestStop|
|00000720| 28 29 29 20 76 2d 3e 65 | 78 69 74 20 3d 20 74 72 |()) v->e|xit = tr|
|00000730| 75 65 3b 0d 09 72 65 74 | 75 72 6e 20 76 2d 3e 65 |ue;..ret|urn v->e|
|00000740| 78 69 74 3b 0d 7d 0d 0d | 73 74 61 74 69 63 20 76 |xit;.}..|static v|
|00000750| 6f 69 64 20 41 64 64 54 | 6f 41 72 72 61 79 28 64 |oid AddT|oArray(d|
|00000760| 6f 75 62 6c 65 2a 2a 20 | 61 72 72 61 79 2c 20 64 |ouble** |array, d|
|00000770| 6f 75 62 6c 65 20 76 61 | 6c 75 65 2c 20 6c 6f 6e |ouble va|lue, lon|
|00000780| 67 20 6e 72 56 61 6c 75 | 65 73 29 0d 09 2f 2a 20 |g nrValu|es)../* |
|00000790| 61 64 64 20 76 61 6c 75 | 65 20 74 6f 20 61 72 72 |add valu|e to arr|
|000007a0| 61 79 2c 20 61 72 72 61 | 79 20 70 72 65 73 65 6e |ay, arra|y presen|
|000007b0| 74 6c 79 20 68 6f 6c 64 | 69 6e 67 20 6e 72 56 61 |tly hold|ing nrVa|
|000007c0| 6c 75 65 73 20 76 61 6c | 75 65 73 20 2a 2f 0d 09 |lues val|ues */..|
|000007d0| 2f 2a 20 72 65 73 69 7a | 67 69 6e 67 20 74 68 65 |/* resiz|ging the|
|000007e0| 20 61 72 72 61 79 20 69 | 66 20 6e 65 63 65 73 73 | array i|f necess|
|000007f0| 61 72 79 20 2a 2f 0d 7b | 0d 09 6c 6f 6e 67 09 72 |ary */.{|..long.r|
|00000800| 65 71 75 69 72 65 64 53 | 69 7a 65 20 3d 20 28 6e |equiredS|ize = (n|
|00000810| 72 56 61 6c 75 65 73 2b | 31 29 2a 73 69 7a 65 6f |rValues+|1)*sizeo|
|00000820| 66 28 64 6f 75 62 6c 65 | 29 3b 0d 09 69 66 20 28 |f(double|);..if (|
|00000830| 47 65 74 48 61 6e 64 6c | 65 53 69 7a 65 28 28 48 |GetHandl|eSize((H|
|00000840| 61 6e 64 6c 65 29 61 72 | 72 61 79 29 20 3c 20 72 |andle)ar|ray) < r|
|00000850| 65 71 75 69 72 65 64 53 | 69 7a 65 29 0d 09 7b 09 |equiredS|ize)..{.|
|00000860| 53 65 74 48 61 6e 64 6c | 65 53 69 7a 65 28 28 48 |SetHandl|eSize((H|
|00000870| 61 6e 64 6c 65 29 61 72 | 72 61 79 2c 20 72 65 71 |andle)ar|ray, req|
|00000880| 75 69 72 65 64 53 69 7a | 65 2b 31 30 32 34 29 3b |uiredSiz|e+1024);|
|00000890| 0d 09 09 69 66 20 28 4d | 65 6d 45 72 72 6f 72 28 |...if (M|emError(|
|000008a0| 29 20 21 3d 20 6e 6f 45 | 72 72 29 20 72 65 74 75 |) != noE|rr) retu|
|000008b0| 72 6e 3b 0d 09 7d 0d 09 | 28 2a 61 72 72 61 79 29 |rn;..}..|(*array)|
|000008c0| 5b 6e 72 56 61 6c 75 65 | 73 5d 20 3d 20 76 61 6c |[nrValue|s] = val|
|000008d0| 75 65 3b 0d 7d 0d 0d 73 | 74 61 74 69 63 20 76 6f |ue;.}..s|tatic vo|
|000008e0| 69 64 20 53 63 61 6c 65 | 28 64 6f 75 62 6c 65 20 |id Scale|(double |
|000008f0| 72 6f 77 2c 20 64 6f 75 | 62 6c 65 20 63 6f 6c 2c |row, dou|ble col,|
|00000900| 20 64 6f 75 62 6c 65 2a | 20 63 6f 6e 73 74 20 78 | double*| const x|
|00000910| 31 2c 20 64 6f 75 62 6c | 65 2a 20 63 6f 6e 73 74 |1, doubl|e* const|
|00000920| 20 78 32 2c 20 4d 79 50 | 6c 6f 74 74 69 6e 67 44 | x2, MyP|lottingD|
|00000930| 61 74 61 2a 20 63 6f 6e | 73 74 20 76 29 0d 09 2f |ata* con|st v)../|
|00000940| 2f 20 63 6f 6e 76 65 72 | 74 73 20 75 6e 69 74 73 |/ conver|ts units|
|00000950| 20 69 6e 20 72 6f 77 73 | 2f 63 6f 6c 75 6d 6e 73 | in rows|/columns|
|00000960| 20 74 6f 20 75 6e 69 74 | 73 20 69 6e 20 78 31 2c | to unit|s in x1,|
|00000970| 20 78 32 0d 7b 0d 09 2a | 78 31 20 3d 20 76 2d 3e | x2.{..*|x1 = v->|
|00000980| 63 4d 69 6e 20 2b 20 63 | 6f 6c 20 2a 20 76 2d 3e |cMin + c|ol * v->|
|00000990| 63 55 6e 69 74 3b 0d 09 | 2a 78 32 20 3d 20 76 2d |cUnit;..|*x2 = v-|
|000009a0| 3e 72 4d 69 6e 20 2b 20 | 72 6f 77 20 2a 20 76 2d |>rMin + |row * v-|
|000009b0| 3e 72 55 6e 69 74 3b 0d | 7d 0d 0d 73 74 61 74 69 |>rUnit;.|}..stati|
|000009c0| 63 20 76 6f 69 64 20 44 | 72 61 77 4f 6e 65 43 6f |c void D|rawOneCo|
|000009d0| 6e 74 6f 75 72 28 64 6f | 75 62 6c 65 20 6c 65 76 |ntour(do|uble lev|
|000009e0| 65 6c 2c 20 6c 6f 6e 67 | 20 73 74 61 72 74 52 2c |el, long| startR,|
|000009f0| 20 6c 6f 6e 67 20 73 74 | 61 72 74 43 2c 20 4d 79 | long st|artC, My|
|00000a00| 50 6c 6f 74 74 69 6e 67 | 44 61 74 61 2a 20 63 6f |Plotting|Data* co|
|00000a10| 6e 73 74 20 76 29 0d 09 | 2f 2a 20 64 72 61 77 20 |nst v)..|/* draw |
|00000a20| 6f 6e 65 20 63 6f 6e 74 | 6f 75 72 20 73 74 61 72 |one cont|our star|
|00000a30| 74 69 6e 67 20 61 74 20 | 72 6f 77 2f 63 6f 6c 20 |ting at |row/col |
|00000a40| 73 74 61 72 74 52 2f 73 | 74 61 72 74 43 20 2a 2f |startR/s|tartC */|
|00000a50| 0d 7b 0d 09 6c 6f 6e 67 | 09 09 72 2c 20 63 3b 0d |.{..long|..r, c;.|
|00000a60| 09 42 6f 6f 6c 65 61 6e | 09 09 64 6f 6e 65 20 3d |.Boolean|..done =|
|00000a70| 20 66 61 6c 73 65 3b 0d | 09 69 6e 74 09 09 09 77 | false;.|.int...w|
|00000a80| 68 65 72 65 3b 09 09 09 | 09 09 09 09 09 2f 2a 20 |here;...|...../* |
|00000a90| 30 20 69 66 20 70 72 65 | 73 65 6e 74 20 70 6f 69 |0 if pre|sent poi|
|00000aa0| 6e 74 20 69 73 20 61 74 | 20 72 69 67 68 74 2c 20 |nt is at| right, |
|00000ab0| 31 20 69 66 20 61 74 20 | 62 75 74 74 6f 6e 20 2a |1 if at |button *|
|00000ac0| 2f 0d 09 4d 65 73 68 50 | 6f 69 6e 74 2a 09 69 6e |/..MeshP|oint*.in|
|00000ad0| 66 6f 20 3d 20 47 65 74 | 4d 65 73 68 50 6f 69 6e |fo = Get|MeshPoin|
|00000ae0| 74 28 73 74 61 72 74 52 | 2c 20 73 74 61 72 74 43 |t(startR|, startC|
|00000af0| 2c 20 76 29 3b 09 09 2f | 2a 20 74 68 65 20 70 6f |, v);../|* the po|
|00000b00| 69 6e 74 20 74 68 61 74 | 20 77 65 20 61 72 65 20 |int that| we are |
|00000b10| 63 75 72 72 65 6e 74 6c | 79 20 6c 6f 6f 6b 69 6e |currentl|y lookin|
|00000b20| 67 20 61 74 20 2a 2f 0d | 09 4d 65 73 68 50 6f 69 |g at */.|.MeshPoi|
|00000b30| 6e 74 2a 09 6f 72 69 67 | 49 6e 66 6f 20 3d 20 69 |nt*.orig|Info = i|
|00000b40| 6e 66 6f 3b 09 09 09 09 | 09 09 09 2f 2a 20 74 68 |nfo;....|.../* th|
|00000b50| 65 20 70 6f 69 6e 74 20 | 77 68 65 72 65 20 77 65 |e point |where we|
|00000b60| 20 73 74 61 72 74 65 64 | 20 2a 2f 0d 09 64 6f 75 | started| */..dou|
|00000b70| 62 6c 65 09 09 2a 2a 78 | 56 61 6c 75 65 73 31 3d |ble..**x|Values1=|
|00000b80| 6e 69 6c 2c 20 2a 2a 79 | 56 61 6c 75 65 73 31 3d |nil, **y|Values1=|
|00000b90| 6e 69 6c 3b 09 09 2f 2a | 20 78 2d 20 61 6e 64 20 |nil;../*| x- and |
|00000ba0| 79 2d 76 61 6c 75 65 73 | 20 61 74 20 66 69 72 73 |y-values| at firs|
|00000bb0| 74 20 62 72 61 6e 63 68 | 20 2a 2f 0d 09 64 6f 75 |t branch| */..dou|
|00000bc0| 62 6c 65 09 09 2a 2a 78 | 56 61 6c 75 65 73 32 3d |ble..**x|Values2=|
|00000bd0| 6e 69 6c 2c 20 2a 2a 79 | 56 61 6c 75 65 73 32 3d |nil, **y|Values2=|
|00000be0| 6e 69 6c 3b 09 09 2f 2a | 20 78 2d 20 61 6e 64 20 |nil;../*| x- and |
|00000bf0| 79 2d 76 61 6c 75 65 73 | 20 61 74 20 73 65 63 6f |y-values| at seco|
|00000c00| 6e 64 20 62 72 61 6e 63 | 68 20 2a 2f 0d 09 64 6f |nd branc|h */..do|
|00000c10| 75 62 6c 65 09 09 2a 2a | 78 56 61 6c 75 65 73 2c |uble..**|xValues,|
|00000c20| 20 2a 2a 79 56 61 6c 75 | 65 73 3b 0d 09 6c 6f 6e | **yValu|es;..lon|
|00000c30| 67 09 09 6e 72 56 61 6c | 75 65 73 31 3d 30 2c 20 |g..nrVal|ues1=0, |
|00000c40| 6e 72 56 61 6c 75 65 73 | 32 3d 30 3b 09 09 09 2f |nrValues|2=0;.../|
|00000c50| 2a 20 6e 75 6d 62 65 72 | 20 6f 66 20 76 61 6c 75 |* number| of valu|
|00000c60| 65 73 20 69 6e 20 66 69 | 72 73 74 20 61 6e 64 20 |es in fi|rst and |
|00000c70| 73 65 63 6f 6e 64 20 62 | 72 61 6e 63 68 20 2a 2f |second b|ranch */|
|00000c80| 0d 09 6c 6f 6e 67 09 09 | 2a 6e 72 56 61 6c 75 65 |..long..|*nrValue|
|00000c90| 73 3b 0d 09 69 6e 74 09 | 09 09 62 72 61 6e 63 68 |s;..int.|..branch|
|00000ca0| 3d 30 3b 09 09 09 09 09 | 09 09 2f 2a 20 30 3a 20 |=0;.....|../* 0: |
|00000cb0| 72 75 6e 6e 69 6e 67 20 | 66 69 72 73 74 20 62 72 |running |first br|
|00000cc0| 61 6e 63 68 2c 20 31 3a | 20 73 65 63 6f 6e 64 20 |anch, 1:| second |
|00000cd0| 62 72 61 6e 63 68 20 2a | 2f 0d 09 69 6e 74 09 09 |branch *|/..int..|
|00000ce0| 09 66 69 72 73 74 57 68 | 65 72 65 3b 09 09 09 09 |.firstWh|ere;....|
|00000cf0| 09 09 09 2f 2a 20 74 68 | 65 20 76 61 6c 75 65 20 |.../* th|e value |
|00000d00| 6f 66 20 22 77 68 65 72 | 65 22 20 61 74 20 74 68 |of "wher|e" at th|
|00000d10| 65 20 66 69 72 73 74 20 | 70 6f 69 6e 74 20 2a 2f |e first |point */|
|00000d20| 0d 0d 09 78 56 61 6c 75 | 65 73 31 20 3d 20 28 64 |...xValu|es1 = (d|
|00000d30| 6f 75 62 6c 65 2a 2a 29 | 4e 65 77 48 61 6e 64 6c |ouble**)|NewHandl|
|00000d40| 65 28 32 35 36 29 3b 0d | 09 78 56 61 6c 75 65 73 |e(256);.|.xValues|
|00000d50| 32 20 3d 20 28 64 6f 75 | 62 6c 65 2a 2a 29 4e 65 |2 = (dou|ble**)Ne|
|00000d60| 77 48 61 6e 64 6c 65 28 | 32 35 36 29 3b 0d 09 79 |wHandle(|256);..y|
|00000d70| 56 61 6c 75 65 73 31 20 | 3d 20 28 64 6f 75 62 6c |Values1 |= (doubl|
|00000d80| 65 2a 2a 29 4e 65 77 48 | 61 6e 64 6c 65 28 32 35 |e**)NewH|andle(25|
|00000d90| 36 29 3b 0d 09 79 56 61 | 6c 75 65 73 32 20 3d 20 |6);..yVa|lues2 = |
|00000da0| 28 64 6f 75 62 6c 65 2a | 2a 29 4e 65 77 48 61 6e |(double*|*)NewHan|
|00000db0| 64 6c 65 28 32 35 36 29 | 3b 0d 09 69 66 20 28 78 |dle(256)|;..if (x|
|00000dc0| 56 61 6c 75 65 73 31 3d | 3d 6e 69 6c 20 7c 7c 20 |Values1=|=nil || |
|00000dd0| 78 56 61 6c 75 65 73 32 | 3d 3d 6e 69 6c 20 7c 7c |xValues2|==nil |||
|00000de0| 20 79 56 61 6c 75 65 73 | 31 3d 3d 6e 69 6c 20 7c | yValues|1==nil ||
|00000df0| 7c 20 79 56 61 6c 75 65 | 73 32 3d 3d 6e 69 6c 29 || yValue|s2==nil)|
|00000e00| 20 67 6f 74 6f 20 64 6f | 6e 65 3b 0d 0d 0d 09 66 | goto do|ne;....f|
|00000e10| 69 72 73 74 57 68 65 72 | 65 20 3d 20 21 69 6e 66 |irstWher|e = !inf|
|00000e20| 6f 2d 3e 61 74 52 69 67 | 68 74 3b 0d 09 77 68 69 |o->atRig|ht;..whi|
|00000e30| 6c 65 20 28 21 54 65 73 | 74 49 66 45 78 69 74 28 |le (!Tes|tIfExit(|
|00000e40| 76 29 20 26 26 20 62 72 | 61 6e 63 68 20 3c 20 32 |v) && br|anch < 2|
|00000e50| 29 09 09 09 2f 2a 20 77 | 68 69 6c 65 20 74 68 65 |).../* w|hile the|
|00000e60| 72 65 27 73 20 73 6f 6d | 65 74 68 69 6e 67 20 74 |re's som|ething t|
|00000e70| 6f 20 74 72 61 63 65 20 | 2a 2f 0d 09 7b 0d 09 09 |o trace |*/..{...|
|00000e80| 69 66 20 28 62 72 61 6e | 63 68 20 3d 3d 20 30 29 |if (bran|ch == 0)|
|00000e90| 0d 09 09 7b 09 77 68 65 | 72 65 20 3d 20 66 69 72 |...{.whe|re = fir|
|00000ea0| 73 74 57 68 65 72 65 3b | 0d 09 09 09 78 56 61 6c |stWhere;|....xVal|
|00000eb0| 75 65 73 20 3d 20 78 56 | 61 6c 75 65 73 31 3b 0d |ues = xV|alues1;.|
|00000ec0| 09 09 09 79 56 61 6c 75 | 65 73 20 3d 20 79 56 61 |...yValu|es = yVa|
|00000ed0| 6c 75 65 73 31 3b 0d 09 | 09 09 6e 72 56 61 6c 75 |lues1;..|..nrValu|
|00000ee0| 65 73 20 3d 20 26 6e 72 | 56 61 6c 75 65 73 31 3b |es = &nr|Values1;|
|00000ef0| 0d 09 09 7d 0d 09 09 65 | 6c 73 65 09 09 09 09 09 |...}...e|lse.....|
|00000f00| 09 09 09 09 09 2f 2a 20 | 69 66 20 73 65 63 6f 6e |...../* |if secon|
|00000f10| 64 20 62 72 61 6e 63 68 | 20 2a 2f 0d 09 09 7b 09 |d branch| */...{.|
|00000f20| 77 68 65 72 65 20 3d 20 | 66 69 72 73 74 57 68 65 |where = |firstWhe|
|00000f30| 72 65 3b 0d 09 09 09 78 | 56 61 6c 75 65 73 20 3d |re;....x|Values =|
|00000f40| 20 78 56 61 6c 75 65 73 | 32 3b 0d 09 09 09 79 56 | xValues|2;....yV|
|00000f50| 61 6c 75 65 73 20 3d 20 | 79 56 61 6c 75 65 73 32 |alues = |yValues2|
|00000f60| 3b 0d 09 09 09 6e 72 56 | 61 6c 75 65 73 20 3d 20 |;....nrV|alues = |
|00000f70| 26 6e 72 56 61 6c 75 65 | 73 32 3b 0d 09 09 7d 0d |&nrValue|s2;...}.|
|00000f80| 09 09 72 3d 73 74 61 72 | 74 52 2c 20 63 3d 73 74 |..r=star|tR, c=st|
|00000f90| 61 72 74 43 3b 0d 09 09 | 69 6e 66 6f 20 3d 20 47 |artC;...|info = G|
|00000fa0| 65 74 4d 65 73 68 50 6f | 69 6e 74 28 72 2c 20 63 |etMeshPo|int(r, c|
|00000fb0| 2c 20 76 29 3b 0d 09 09 | 62 72 61 6e 63 68 20 2b |, v);...|branch +|
|00000fc0| 3d 20 31 3b 0d 09 09 77 | 68 69 6c 65 20 28 21 54 |= 1;...w|hile (!T|
|00000fd0| 65 73 74 49 66 45 78 69 | 74 28 76 29 29 09 09 09 |estIfExi|t(v))...|
|00000fe0| 09 09 09 2f 2a 20 74 72 | 61 63 65 20 6f 6e 65 20 |.../* tr|ace one |
|00000ff0| 62 72 61 6e 63 68 20 2a | 2f 0d 09 09 7b 0d 09 09 |branch *|/...{...|
|00001000| 09 4d 65 73 68 50 6f 69 | 6e 74 2a 09 6e 65 78 74 |.MeshPoi|nt*.next|
|00001010| 3b 0d 09 09 09 64 6f 75 | 62 6c 65 09 09 64 3b 0d |;....dou|ble..d;.|
|00001020| 09 0d 09 09 09 69 66 20 | 28 62 72 61 6e 63 68 20 |.....if |(branch |
|00001030| 3d 3d 20 31 20 26 26 20 | 2a 6e 72 56 61 6c 75 65 |== 1 && |*nrValue|
|00001040| 73 20 3d 3d 20 32 29 09 | 09 2f 2a 20 61 66 74 65 |s == 2).|./* afte|
|00001050| 72 20 74 77 6f 20 73 74 | 65 70 73 20 2a 2f 0d 09 |r two st|eps */..|
|00001060| 09 09 7b 09 69 66 20 28 | 66 69 72 73 74 57 68 65 |..{.if (|firstWhe|
|00001070| 72 65 20 3d 3d 20 30 29 | 09 09 09 09 2f 2a 20 77 |re == 0)|..../* w|
|00001080| 65 20 6d 75 73 74 20 73 | 65 74 20 6f 72 69 67 69 |e must s|et origi|
|00001090| 6e 61 6c 20 70 6f 69 6e | 74 20 62 61 63 6b 20 74 |nal poin|t back t|
|000010a0| 6f 20 74 72 75 65 20 28 | 69 66 20 6c 6f 6f 70 20 |o true (|if loop |
|000010b0| 69 73 20 63 6c 6f 73 65 | 64 2c 20 6f 72 69 67 69 |is close|d, origi|
|000010c0| 6e 61 6c 20 70 6f 69 6e | 74 20 69 73 20 75 73 65 |nal poin|t is use|
|000010d0| 64 20 74 77 69 63 65 29 | 20 2a 2f 0d 09 09 09 09 |d twice)| */.....|
|000010e0| 09 6f 72 69 67 49 6e 66 | 6f 2d 3e 61 74 52 69 67 |.origInf|o->atRig|
|000010f0| 68 74 20 3d 20 74 72 75 | 65 3b 0d 09 09 09 09 65 |ht = tru|e;.....e|
|00001100| 6c 73 65 0d 09 09 09 09 | 09 6f 72 69 67 49 6e 66 |lse.....|.origInf|
|00001110| 6f 2d 3e 61 74 42 6f 74 | 74 6f 6d 20 3d 20 74 72 |o->atBot|tom = tr|
|00001120| 75 65 3b 0d 09 09 09 7d | 0d 0d 09 09 09 69 66 20 |ue;....}|.....if |
|00001130| 28 77 68 65 72 65 3d 3d | 30 29 09 09 09 09 09 09 |(where==|0)......|
|00001140| 09 2f 2a 20 69 66 20 77 | 65 20 61 72 65 20 74 6f |./* if w|e are to|
|00001150| 20 74 68 65 20 72 69 67 | 68 74 20 6f 66 20 74 68 | the rig|ht of th|
|00001160| 65 20 70 72 65 73 65 6e | 74 20 6d 65 73 68 20 70 |e presen|t mesh p|
|00001170| 6f 69 6e 74 20 2a 2f 0d | 09 09 09 7b 0d 09 09 09 |oint */.|...{....|
|00001180| 09 69 6e 66 6f 2d 3e 61 | 74 52 69 67 68 74 20 3d |.info->a|tRight =|
|00001190| 20 66 61 6c 73 65 3b 09 | 09 09 09 2f 2a 20 64 65 | false;.|.../* de|
|000011a0| 6c 65 74 65 20 74 68 69 | 73 20 70 6f 69 6e 74 20 |lete thi|s point |
|000011b0| 2d 20 69 74 20 69 73 20 | 74 72 61 63 65 64 20 6e |- it is |traced n|
|000011c0| 6f 77 20 2a 2f 0d 09 09 | 09 09 6e 65 78 74 20 3d |ow */...|..next =|
|000011d0| 20 47 65 74 4d 65 73 68 | 50 6f 69 6e 74 28 72 2c | GetMesh|Point(r,|
|000011e0| 20 63 2b 31 2c 20 76 29 | 3b 0d 09 09 09 09 69 66 | c+1, v)|;.....if|
|000011f0| 20 28 6e 65 78 74 2d 3e | 76 61 6c 69 64 3d 3d 66 | (next->|valid==f|
|00001200| 61 6c 73 65 29 0d 09 09 | 09 09 09 62 72 65 61 6b |alse)...|...break|
|00001210| 3b 0d 09 0d 09 09 09 09 | 64 20 3d 20 41 42 53 28 |;.......|d = ABS(|
|00001220| 6e 65 78 74 2d 3e 76 61 | 6c 75 65 20 2d 20 69 6e |next->va|lue - in|
|00001230| 66 6f 2d 3e 76 61 6c 75 | 65 29 3b 0d 09 09 09 09 |fo->valu|e);.....|
|00001240| 69 66 20 28 64 20 3d 3d | 20 30 29 0d 09 09 09 09 |if (d ==| 0).....|
|00001250| 09 62 72 65 61 6b 3b 0d | 09 09 09 09 0d 09 09 09 |.break;.|........|
|00001260| 09 41 64 64 54 6f 41 72 | 72 61 79 28 78 56 61 6c |.AddToAr|ray(xVal|
|00001270| 75 65 73 2c 20 72 2c 20 | 2a 6e 72 56 61 6c 75 65 |ues, r, |*nrValue|
|00001280| 73 29 3b 0d 09 09 09 09 | 41 64 64 54 6f 41 72 72 |s);.....|AddToArr|
|00001290| 61 79 28 79 56 61 6c 75 | 65 73 2c 20 63 20 2b 20 |ay(yValu|es, c + |
|000012a0| 41 42 53 28 69 6e 66 6f | 2d 3e 76 61 6c 75 65 2d |ABS(info|->value-|
|000012b0| 6c 65 76 65 6c 29 2f 64 | 2c 20 28 2a 6e 72 56 61 |level)/d|, (*nrVa|
|000012c0| 6c 75 65 73 29 2b 2b 29 | 3b 0d 09 0d 09 09 09 09 |lues)++)|;.......|
|000012d0| 69 66 20 28 69 6e 66 6f | 2d 3e 61 74 42 6f 74 74 |if (info|->atBott|
|000012e0| 6f 6d 29 0d 09 09 09 09 | 7b 09 77 68 65 72 65 20 |om).....|{.where |
|000012f0| 3d 20 31 3b 0d 09 09 09 | 09 09 63 6f 6e 74 69 6e |= 1;....|..contin|
|00001300| 75 65 3b 0d 09 09 09 09 | 7d 0d 0d 09 09 09 09 69 |ue;.....|}......i|
|00001310| 66 20 28 63 2b 31 20 21 | 3d 20 76 2d 3e 6e 72 43 |f (c+1 !|= v->nrC|
|00001320| 6f 6c 73 29 09 09 09 09 | 09 2f 2a 20 69 66 20 6e |ols)....|./* if n|
|00001330| 6f 74 20 79 65 74 20 61 | 74 20 72 69 67 68 74 20 |ot yet a|t right |
|00001340| 6f 66 20 67 72 69 64 20 | 2a 2f 0d 09 09 09 09 7b |of grid |*/.....{|
|00001350| 09 6e 65 78 74 20 3d 20 | 47 65 74 4d 65 73 68 50 |.next = |GetMeshP|
|00001360| 6f 69 6e 74 28 72 2c 20 | 63 2b 31 2c 20 76 29 3b |oint(r, |c+1, v);|
|00001370| 0d 09 09 09 09 09 69 66 | 20 28 6e 65 78 74 2d 3e |......if| (next->|
|00001380| 76 61 6c 69 64 20 26 26 | 20 6e 65 78 74 2d 3e 61 |valid &&| next->a|
|00001390| 74 42 6f 74 74 6f 6d 29 | 0d 09 09 09 09 09 7b 09 |tBottom)|......{.|
|000013a0| 63 20 2b 3d 20 31 3b 0d | 09 09 09 09 09 09 77 68 |c += 1;.|......wh|
|000013b0| 65 72 65 20 3d 20 31 3b | 0d 09 09 09 09 09 09 69 |ere = 1;|.......i|
|000013c0| 6e 66 6f 20 3d 20 6e 65 | 78 74 3b 0d 09 09 09 09 |nfo = ne|xt;.....|
|000013d0| 09 09 63 6f 6e 74 69 6e | 75 65 3b 0d 09 09 09 09 |..contin|ue;.....|
|000013e0| 09 7d 0d 09 09 09 09 7d | 0d 09 09 09 09 69 66 20 |.}.....}|.....if |
|000013f0| 28 72 2b 31 20 21 3d 20 | 76 2d 3e 6e 72 52 6f 77 |(r+1 != |v->nrRow|
|00001400| 73 29 09 09 09 09 09 2f | 2a 20 69 66 20 6e 6f 74 |s)...../|* if not|
|00001410| 20 61 74 20 62 6f 74 74 | 6f 6d 20 6f 66 20 67 72 | at bott|om of gr|
|00001420| 69 64 20 2a 2f 0d 09 09 | 09 09 7b 09 6e 65 78 74 |id */...|..{.next|
|00001430| 20 3d 20 47 65 74 4d 65 | 73 68 50 6f 69 6e 74 28 | = GetMe|shPoint(|
|00001440| 72 2b 31 2c 20 63 2c 20 | 76 29 3b 0d 09 09 09 09 |r+1, c, |v);.....|
|00001450| 09 69 66 20 28 6e 65 78 | 74 2d 3e 76 61 6c 69 64 |.if (nex|t->valid|
|00001460| 20 26 26 20 6e 65 78 74 | 2d 3e 61 74 52 69 67 68 | && next|->atRigh|
|00001470| 74 29 0d 09 09 09 09 09 | 7b 09 72 20 2b 3d 20 31 |t)......|{.r += 1|
|00001480| 3b 0d 09 09 09 09 09 09 | 77 68 65 72 65 20 3d 20 |;.......|where = |
|00001490| 30 3b 0d 09 09 09 09 09 | 09 69 6e 66 6f 20 3d 20 |0;......|.info = |
|000014a0| 6e 65 78 74 3b 0d 09 09 | 09 09 09 09 63 6f 6e 74 |next;...|....cont|
|000014b0| 69 6e 75 65 3b 0d 09 09 | 09 09 09 7d 0d 09 09 09 |inue;...|...}....|
|000014c0| 09 7d 0d 09 09 09 09 69 | 66 20 28 72 20 21 3d 20 |.}.....i|f (r != |
|000014d0| 30 29 09 09 09 09 09 09 | 09 09 2f 2a 20 69 66 20 |0)......|../* if |
|000014e0| 6e 6f 74 20 61 74 20 74 | 6f 70 20 6f 66 20 67 72 |not at t|op of gr|
|000014f0| 69 64 20 2a 2f 0d 09 09 | 09 09 7b 0d 09 09 09 09 |id */...|..{.....|
|00001500| 09 6e 65 78 74 20 3d 20 | 47 65 74 4d 65 73 68 50 |.next = |GetMeshP|
|00001510| 6f 69 6e 74 28 72 2d 31 | 2c 20 63 2c 20 76 29 3b |oint(r-1|, c, v);|
|00001520| 0d 09 09 09 09 09 69 66 | 20 28 6e 65 78 74 2d 3e |......if| (next->|
|00001530| 76 61 6c 69 64 20 26 26 | 20 28 6e 65 78 74 2d 3e |valid &&| (next->|
|00001540| 61 74 52 69 67 68 74 20 | 7c 7c 20 6e 65 78 74 2d |atRight ||| next-|
|00001550| 3e 61 74 42 6f 74 74 6f | 6d 29 29 0d 09 09 09 09 |>atBotto|m)).....|
|00001560| 09 7b 09 69 6e 66 6f 20 | 3d 20 6e 65 78 74 3b 0d |.{.info |= next;.|
|00001570| 09 09 09 09 09 09 72 20 | 2d 3d 20 31 3b 0d 09 09 |......r |-= 1;...|
|00001580| 09 09 09 09 77 68 65 72 | 65 20 3d 20 21 69 6e 66 |....wher|e = !inf|
|00001590| 6f 2d 3e 61 74 52 69 67 | 68 74 3b 0d 09 09 09 09 |o->atRig|ht;.....|
|000015a0| 09 09 63 6f 6e 74 69 6e | 75 65 3b 0d 09 09 09 09 |..contin|ue;.....|
|000015b0| 09 7d 0d 09 09 09 09 09 | 69 66 20 28 63 2b 31 20 |.}......|if (c+1 |
|000015c0| 21 3d 20 76 2d 3e 6e 72 | 43 6f 6c 73 29 09 09 09 |!= v->nr|Cols)...|
|000015d0| 09 2f 2a 20 69 66 20 6e | 6f 74 20 79 65 74 20 61 |./* if n|ot yet a|
|000015e0| 74 20 72 69 67 68 74 20 | 6f 66 20 67 72 69 64 20 |t right |of grid |
|000015f0| 2a 2f 0d 09 09 09 09 09 | 7b 09 6e 65 78 74 20 3d |*/......|{.next =|
|00001600| 20 47 65 74 4d 65 73 68 | 50 6f 69 6e 74 28 72 2d | GetMesh|Point(r-|
|00001610| 31 2c 20 63 2b 31 2c 20 | 76 29 3b 0d 09 09 09 09 |1, c+1, |v);.....|
|00001620| 09 09 69 66 20 28 6e 65 | 78 74 2d 3e 76 61 6c 69 |..if (ne|xt->vali|
|00001630| 64 20 26 26 20 6e 65 78 | 74 2d 3e 61 74 42 6f 74 |d && nex|t->atBot|
|00001640| 74 6f 6d 29 0d 09 09 09 | 09 09 09 7b 09 72 20 2d |tom)....|...{.r -|
|00001650| 3d 20 31 3b 0d 09 09 09 | 09 09 09 09 63 20 2b 3d |= 1;....|....c +=|
|00001660| 20 31 3b 0d 09 09 09 09 | 09 09 09 77 68 65 72 65 | 1;.....|...where|
|00001670| 20 3d 20 31 3b 0d 09 09 | 09 09 09 09 09 69 6e 66 | = 1;...|.....inf|
|00001680| 6f 20 3d 20 6e 65 78 74 | 3b 0d 09 09 09 09 09 09 |o = next|;.......|
|00001690| 09 63 6f 6e 74 69 6e 75 | 65 3b 0d 09 09 09 09 09 |.continu|e;......|
|000016a0| 09 7d 0d 09 09 09 09 09 | 7d 0d 09 09 09 09 7d 0d |.}......|}.....}.|
|000016b0| 09 09 09 09 62 72 65 61 | 6b 3b 09 09 09 09 09 09 |....brea|k;......|
|000016c0| 09 09 09 2f 2a 20 64 69 | 64 20 6e 6f 74 20 66 69 |.../* di|d not fi|
|000016d0| 6e 64 20 63 6f 6e 74 69 | 6e 75 61 74 69 6f 6e 20 |nd conti|nuation |
|000016e0| 2a 2f 0d 09 09 09 7d 2f | 2f 69 66 20 77 68 65 72 |*/....}/|/if wher|
|000016f0| 65 3d 3d 30 0d 09 09 09 | 65 6c 73 65 0d 09 09 09 |e==0....|else....|
|00001700| 7b 0d 09 09 09 09 69 6e | 66 6f 2d 3e 61 74 42 6f |{.....in|fo->atBo|
|00001710| 74 74 6f 6d 20 3d 20 66 | 61 6c 73 65 3b 09 09 09 |ttom = f|alse;...|
|00001720| 09 09 2f 2a 20 64 65 6c | 65 74 65 20 74 68 69 73 |../* del|ete this|
|00001730| 20 70 6f 69 6e 74 20 2d | 20 69 74 20 69 73 20 74 | point -| it is t|
|00001740| 72 61 63 65 64 20 6e 6f | 77 20 2a 2f 0d 09 09 09 |raced no|w */....|
|00001750| 09 6e 65 78 74 20 3d 20 | 47 65 74 4d 65 73 68 50 |.next = |GetMeshP|
|00001760| 6f 69 6e 74 28 72 2b 31 | 2c 20 63 2c 20 76 29 3b |oint(r+1|, c, v);|
|00001770| 0d 09 09 09 09 69 66 20 | 28 6e 65 78 74 2d 3e 76 |.....if |(next->v|
|00001780| 61 6c 69 64 3d 3d 66 61 | 6c 73 65 29 0d 09 09 09 |alid==fa|lse)....|
|00001790| 09 09 62 72 65 61 6b 3b | 0d 09 0d 09 09 09 09 64 |..break;|.......d|
|000017a0| 20 3d 20 41 42 53 28 6e | 65 78 74 2d 3e 76 61 6c | = ABS(n|ext->val|
|000017b0| 75 65 20 2d 20 69 6e 66 | 6f 2d 3e 76 61 6c 75 65 |ue - inf|o->value|
|000017c0| 29 3b 0d 09 09 09 09 69 | 66 20 28 64 20 3d 3d 20 |);.....i|f (d == |
|000017d0| 30 29 0d 09 09 09 09 09 | 62 72 65 61 6b 3b 0d 09 |0)......|break;..|
|000017e0| 09 09 09 0d 09 09 09 09 | 41 64 64 54 6f 41 72 72 |........|AddToArr|
|000017f0| 61 79 28 78 56 61 6c 75 | 65 73 2c 20 72 20 2b 20 |ay(xValu|es, r + |
|00001800| 41 42 53 28 69 6e 66 6f | 2d 3e 76 61 6c 75 65 2d |ABS(info|->value-|
|00001810| 6c 65 76 65 6c 29 2f 64 | 2c 20 28 2a 6e 72 56 61 |level)/d|, (*nrVa|
|00001820| 6c 75 65 73 29 29 3b 0d | 09 09 09 09 41 64 64 54 |lues));.|....AddT|
|00001830| 6f 41 72 72 61 79 28 79 | 56 61 6c 75 65 73 2c 20 |oArray(y|Values, |
|00001840| 63 2c 20 28 2a 6e 72 56 | 61 6c 75 65 73 29 2b 2b |c, (*nrV|alues)++|
|00001850| 29 3b 0d 09 0d 09 09 09 | 09 69 66 20 28 69 6e 66 |);......|.if (inf|
|00001860| 6f 2d 3e 61 74 52 69 67 | 68 74 29 0d 09 09 09 09 |o->atRig|ht).....|
|00001870| 7b 09 77 68 65 72 65 20 | 3d 20 30 3b 0d 09 09 09 |{.where |= 0;....|
|00001880| 09 09 63 6f 6e 74 69 6e | 75 65 3b 0d 09 09 09 09 |..contin|ue;.....|
|00001890| 7d 0d 0d 09 09 09 09 69 | 66 20 28 72 2b 31 20 21 |}......i|f (r+1 !|
|000018a0| 3d 20 76 2d 3e 6e 72 52 | 6f 77 73 29 09 09 09 09 |= v->nrR|ows)....|
|000018b0| 09 2f 2a 20 69 66 20 6e | 6f 74 20 79 65 74 20 61 |./* if n|ot yet a|
|000018c0| 74 20 62 6f 74 74 6f 6d | 20 6f 66 20 67 72 69 64 |t bottom| of grid|
|000018d0| 20 2a 2f 0d 09 09 09 09 | 7b 09 6e 65 78 74 20 3d | */.....|{.next =|
|000018e0| 20 47 65 74 4d 65 73 68 | 50 6f 69 6e 74 28 72 2b | GetMesh|Point(r+|
|000018f0| 31 2c 20 63 2c 20 76 29 | 3b 0d 09 09 09 09 09 69 |1, c, v)|;......i|
|00001900| 66 20 28 6e 65 78 74 2d | 3e 76 61 6c 69 64 20 26 |f (next-|>valid &|
|00001910| 26 20 6e 65 78 74 2d 3e | 61 74 52 69 67 68 74 29 |& next->|atRight)|
|00001920| 0d 09 09 09 09 09 7b 09 | 09 72 20 2b 3d 20 31 3b |......{.|.r += 1;|
|00001930| 0d 09 09 09 09 09 09 09 | 77 68 65 72 65 20 3d 20 |........|where = |
|00001940| 30 3b 0d 09 09 09 09 09 | 09 09 69 6e 66 6f 20 3d |0;......|..info =|
|00001950| 20 6e 65 78 74 3b 0d 09 | 09 09 09 09 09 09 63 6f | next;..|......co|
|00001960| 6e 74 69 6e 75 65 3b 0d | 09 09 09 09 09 7d 0d 09 |ntinue;.|.....}..|
|00001970| 09 09 09 7d 0d 09 09 09 | 09 69 66 20 28 63 2b 31 |...}....|.if (c+1|
|00001980| 20 21 3d 20 76 2d 3e 6e | 72 43 6f 6c 73 29 09 09 | != v->n|rCols)..|
|00001990| 09 09 09 2f 2a 20 69 66 | 20 6e 6f 74 20 79 65 74 |.../* if| not yet|
|000019a0| 20 61 74 20 72 69 67 68 | 74 20 6f 66 20 67 72 69 | at righ|t of gri|
|000019b0| 64 20 2a 2f 0d 09 09 09 | 09 7b 09 6e 65 78 74 20 |d */....|.{.next |
|000019c0| 3d 20 47 65 74 4d 65 73 | 68 50 6f 69 6e 74 28 72 |= GetMes|hPoint(r|
|000019d0| 2c 20 63 2b 31 2c 20 76 | 29 3b 0d 09 09 09 09 09 |, c+1, v|);......|
|000019e0| 69 66 20 28 6e 65 78 74 | 2d 3e 76 61 6c 69 64 20 |if (next|->valid |
|000019f0| 26 26 20 6e 65 78 74 2d | 3e 61 74 42 6f 74 74 6f |&& next-|>atBotto|
|00001a00| 6d 29 0d 09 09 09 09 09 | 7b 09 63 20 2b 3d 20 31 |m)......|{.c += 1|
|00001a10| 3b 0d 09 09 09 09 09 09 | 77 68 65 72 65 20 3d 20 |;.......|where = |
|00001a20| 31 3b 0d 09 09 09 09 09 | 09 69 6e 66 6f 20 3d 20 |1;......|.info = |
|00001a30| 6e 65 78 74 3b 0d 09 09 | 09 09 09 09 63 6f 6e 74 |next;...|....cont|
|00001a40| 69 6e 75 65 3b 0d 09 09 | 09 09 09 7d 0d 09 09 09 |inue;...|...}....|
|00001a50| 09 7d 0d 09 09 09 09 69 | 66 20 28 63 20 21 3d 20 |.}.....i|f (c != |
|00001a60| 30 29 09 09 09 09 09 09 | 09 09 2f 2f 20 69 66 20 |0)......|..// if |
|00001a70| 6e 6f 74 20 61 74 20 6c | 65 66 74 20 6f 66 20 67 |not at l|eft of g|
|00001a80| 72 69 64 0d 09 09 09 09 | 7b 09 6e 65 78 74 20 3d |rid.....|{.next =|
|00001a90| 20 47 65 74 4d 65 73 68 | 50 6f 69 6e 74 28 72 2c | GetMesh|Point(r,|
|00001aa0| 20 63 2d 31 2c 20 76 29 | 3b 0d 09 09 09 09 09 69 | c-1, v)|;......i|
|00001ab0| 66 20 28 6e 65 78 74 2d | 3e 76 61 6c 69 64 20 26 |f (next-|>valid &|
|00001ac0| 26 20 28 6e 65 78 74 2d | 3e 61 74 52 69 67 68 74 |& (next-|>atRight|
|00001ad0| 20 7c 7c 20 6e 65 78 74 | 2d 3e 61 74 42 6f 74 74 | || next|->atBott|
|00001ae0| 6f 6d 29 29 0d 09 09 09 | 09 09 7b 09 69 6e 66 6f |om))....|..{.info|
|00001af0| 20 3d 20 6e 65 78 74 3b | 0d 09 09 09 09 09 09 63 | = next;|.......c|
|00001b00| 20 2d 3d 20 31 3b 0d 09 | 09 09 09 09 09 77 68 65 | -= 1;..|.....whe|
|00001b10| 72 65 20 3d 20 21 69 6e | 66 6f 2d 3e 61 74 52 69 |re = !in|fo->atRi|
|00001b20| 67 68 74 3b 0d 09 09 09 | 09 09 09 63 6f 6e 74 69 |ght;....|...conti|
|00001b30| 6e 75 65 3b 0d 09 09 09 | 09 09 7d 0d 09 09 09 09 |nue;....|..}.....|
|00001b40| 09 69 66 20 28 72 2b 31 | 20 21 3d 20 76 2d 3e 6e |.if (r+1| != v->n|
|00001b50| 72 52 6f 77 73 29 09 09 | 09 09 2f 2a 20 69 66 20 |rRows)..|../* if |
|00001b60| 6e 6f 74 20 79 65 74 20 | 61 74 20 62 6f 74 74 6f |not yet |at botto|
|00001b70| 6d 20 6f 66 20 67 72 69 | 64 20 2a 2f 0d 09 09 09 |m of gri|d */....|
|00001b80| 09 09 7b 09 6e 65 78 74 | 20 3d 20 47 65 74 4d 65 |..{.next| = GetMe|
|00001b90| 73 68 50 6f 69 6e 74 28 | 72 2b 31 2c 20 63 2d 31 |shPoint(|r+1, c-1|
|00001ba0| 2c 20 76 29 3b 0d 09 09 | 09 09 09 09 69 66 20 28 |, v);...|....if (|
|00001bb0| 6e 65 78 74 2d 3e 76 61 | 6c 69 64 20 26 26 20 6e |next->va|lid && n|
|00001bc0| 65 78 74 2d 3e 61 74 52 | 69 67 68 74 29 0d 09 09 |ext->atR|ight)...|
|00001bd0| 09 09 09 09 7b 09 63 20 | 2d 3d 20 31 3b 0d 09 09 |....{.c |-= 1;...|
|00001be0| 09 09 09 09 09 72 20 2b | 3d 20 31 3b 0d 09 09 09 |.....r +|= 1;....|
|00001bf0| 09 09 09 09 77 68 65 72 | 65 20 3d 20 30 3b 0d 09 |....wher|e = 0;..|
|00001c00| 09 09 09 09 09 09 69 6e | 66 6f 20 3d 20 6e 65 78 |......in|fo = nex|
|00001c10| 74 3b 0d 09 09 09 09 09 | 09 09 63 6f 6e 74 69 6e |t;......|..contin|
|00001c20| 75 65 3b 0d 09 09 09 09 | 09 09 7d 0d 09 09 09 09 |ue;.....|..}.....|
|00001c30| 09 7d 0d 09 09 09 09 7d | 0d 09 09 09 09 62 72 65 |.}.....}|.....bre|
|00001c40| 61 6b 3b 09 09 09 09 09 | 09 09 09 09 2f 2a 20 64 |ak;.....|..../* d|
|00001c50| 69 64 20 6e 6f 74 20 66 | 69 6e 64 20 63 6f 6e 74 |id not f|ind cont|
|00001c60| 69 6e 75 61 74 69 6f 6e | 20 2a 2f 0d 09 09 09 7d |inuation| */....}|
|00001c70| 2f 2f 69 66 20 77 68 65 | 72 65 3d 3d 30 0d 09 09 |//if whe|re==0...|
|00001c80| 7d 2f 2f 77 68 69 6c 65 | 0d 09 7d 2f 2f 77 68 69 |}//while|..}//whi|
|00001c90| 6c 65 0d 09 69 66 20 28 | 54 65 73 74 49 66 45 78 |le..if (|TestIfEx|
|00001ca0| 69 74 28 76 29 29 20 67 | 6f 74 6f 20 64 6f 6e 65 |it(v)) g|oto done|
|00001cb0| 3b 0d 09 7b 0d 09 09 42 | 6f 6f 6c 65 61 6e 09 66 |;..{...B|oolean.f|
|00001cc0| 69 72 73 74 3d 74 72 75 | 65 3b 0d 09 09 6c 6f 6e |irst=tru|e;...lon|
|00001cd0| 67 09 69 3b 0d 09 09 64 | 6f 75 62 6c 65 09 78 2c |g.i;...d|ouble.x,|
|00001ce0| 20 79 2c 20 78 4c 61 73 | 74 2c 20 79 4c 61 73 74 | y, xLas|t, yLast|
|00001cf0| 3b 0d 0d 09 09 69 66 20 | 28 6e 72 56 61 6c 75 65 |;....if |(nrValue|
|00001d00| 73 31 3e 30 29 0d 09 09 | 7b 09 66 6f 72 20 28 69 |s1>0)...|{.for (i|
|00001d10| 3d 6e 72 56 61 6c 75 65 | 73 31 2d 31 3b 20 69 3e |=nrValue|s1-1; i>|
|00001d20| 3d 30 3b 20 69 2d 2d 29 | 0d 09 09 09 7b 09 53 63 |=0; i--)|....{.Sc|
|00001d30| 61 6c 65 28 28 2a 78 56 | 61 6c 75 65 73 31 29 5b |ale((*xV|alues1)[|
|00001d40| 69 5d 2c 20 28 2a 79 56 | 61 6c 75 65 73 31 29 5b |i], (*yV|alues1)[|
|00001d50| 69 5d 2c 20 26 78 2c 20 | 26 79 2c 20 76 29 3b 0d |i], &x, |&y, v);.|
|00001d60| 09 09 09 09 69 66 20 28 | 66 69 72 73 74 29 20 47 |....if (|first) G|
|00001d70| 72 4d 6f 76 65 54 6f 28 | 78 2c 20 79 29 3b 0d 09 |rMoveTo(|x, y);..|
|00001d80| 09 09 09 65 6c 73 65 20 | 47 72 4c 69 6e 65 54 6f |...else |GrLineTo|
|00001d90| 28 78 2c 20 79 29 3b 0d | 09 09 09 09 66 69 72 73 |(x, y);.|....firs|
|00001da0| 74 20 3d 20 66 61 6c 73 | 65 3b 0d 09 09 09 7d 0d |t = fals|e;....}.|
|00001db0| 09 09 09 78 4c 61 73 74 | 20 3d 20 78 3b 0d 09 09 |...xLast| = x;...|
|00001dc0| 09 79 4c 61 73 74 20 3d | 20 79 3b 0d 09 09 7d 0d |.yLast =| y;...}.|
|00001dd0| 09 09 66 6f 72 20 28 69 | 3d 30 3b 20 69 20 3c 20 |..for (i|=0; i < |
|00001de0| 6e 72 56 61 6c 75 65 73 | 32 3b 20 69 2b 2b 29 0d |nrValues|2; i++).|
|00001df0| 09 09 7b 09 53 63 61 6c | 65 28 28 2a 78 56 61 6c |..{.Scal|e((*xVal|
|00001e00| 75 65 73 32 29 5b 69 5d | 2c 20 28 2a 79 56 61 6c |ues2)[i]|, (*yVal|
|00001e10| 75 65 73 32 29 5b 69 5d | 2c 20 26 78 2c 20 26 79 |ues2)[i]|, &x, &y|
|00001e20| 2c 20 76 29 3b 0d 09 09 | 09 69 66 20 28 6e 72 56 |, v);...|.if (nrV|
|00001e30| 61 6c 75 65 73 31 20 3e | 20 30 20 26 26 20 69 20 |alues1 >| 0 && i |
|00001e40| 3d 3d 20 30 29 09 09 09 | 09 2f 2a 20 69 66 20 74 |== 0)...|./* if t|
|00001e50| 68 65 20 66 69 72 73 74 | 20 70 6f 69 6e 74 20 6f |he first| point o|
|00001e60| 66 20 74 68 69 73 20 62 | 72 61 6e 63 68 20 2a 2f |f this b|ranch */|
|00001e70| 0d 09 09 09 7b 09 69 66 | 20 28 78 3d 3d 78 4c 61 |....{.if| (x==xLa|
|00001e80| 73 74 20 26 26 20 79 3d | 3d 79 4c 61 73 74 29 09 |st && y=|=yLast).|
|00001e90| 09 09 09 2f 2a 20 69 66 | 20 6c 61 73 74 20 62 72 |.../* if| last br|
|00001ea0| 61 6e 63 68 20 65 6e 64 | 65 64 20 77 69 74 68 20 |anch end|ed with |
|00001eb0| 73 61 6d 65 20 70 6f 69 | 6e 74 20 61 73 20 2a 2f |same poi|nt as */|
|00001ec0| 0d 09 09 09 09 09 63 6f | 6e 74 69 6e 75 65 3b 09 |......co|ntinue;.|
|00001ed0| 09 09 09 09 09 09 2f 2a | 20 74 68 69 73 20 70 6f |....../*| this po|
|00001ee0| 69 6e 74 20 73 74 61 72 | 74 73 20 77 69 74 68 20 |int star|ts with |
|00001ef0| 2a 2f 0d 09 09 09 7d 0d | 09 09 09 69 66 20 28 66 |*/....}.|...if (f|
|00001f00| 69 72 73 74 29 20 47 72 | 4d 6f 76 65 54 6f 28 78 |irst) Gr|MoveTo(x|
|00001f10| 2c 20 79 29 3b 0d 09 09 | 09 65 6c 73 65 20 47 72 |, y);...|.else Gr|
|00001f20| 4c 69 6e 65 54 6f 28 78 | 2c 20 79 29 3b 0d 09 09 |LineTo(x|, y);...|
|00001f30| 09 66 69 72 73 74 20 3d | 20 66 61 6c 73 65 3b 0d |.first =| false;.|
|00001f40| 09 09 7d 0d 09 7d 0d 0d | 64 6f 6e 65 3a 0d 09 69 |..}..}..|done:..i|
|00001f50| 66 20 28 78 56 61 6c 75 | 65 73 31 20 21 3d 20 6e |f (xValu|es1 != n|
|00001f60| 69 6c 29 20 44 69 73 70 | 6f 73 65 48 61 6e 64 6c |il) Disp|oseHandl|
|00001f70| 65 28 28 48 61 6e 64 6c | 65 29 78 56 61 6c 75 65 |e((Handl|e)xValue|
|00001f80| 73 31 29 3b 0d 09 69 66 | 20 28 78 56 61 6c 75 65 |s1);..if| (xValue|
|00001f90| 73 32 20 21 3d 20 6e 69 | 6c 29 20 44 69 73 70 6f |s2 != ni|l) Dispo|
|00001fa0| 73 65 48 61 6e 64 6c 65 | 28 28 48 61 6e 64 6c 65 |seHandle|((Handle|
|00001fb0| 29 78 56 61 6c 75 65 73 | 32 29 3b 0d 09 69 66 20 |)xValues|2);..if |
|00001fc0| 28 79 56 61 6c 75 65 73 | 31 20 21 3d 20 6e 69 6c |(yValues|1 != nil|
|00001fd0| 29 20 44 69 73 70 6f 73 | 65 48 61 6e 64 6c 65 28 |) Dispos|eHandle(|
|00001fe0| 28 48 61 6e 64 6c 65 29 | 79 56 61 6c 75 65 73 31 |(Handle)|yValues1|
|00001ff0| 29 3b 0d 09 69 66 20 28 | 79 56 61 6c 75 65 73 32 |);..if (|yValues2|
|00002000| 20 21 3d 20 6e 69 6c 29 | 20 44 69 73 70 6f 73 65 | != nil)| Dispose|
|00002010| 48 61 6e 64 6c 65 28 28 | 48 61 6e 64 6c 65 29 79 |Handle((|Handle)y|
|00002020| 56 61 6c 75 65 73 32 29 | 3b 0d 7d 0d 0d 0d 0d 50 |Values2)|;.}....P|
|00002030| 6c 6f 74 74 69 6e 67 44 | 61 74 61 2a 20 49 6e 69 |lottingD|ata* Ini|
|00002040| 74 43 6f 6e 74 6f 75 72 | 50 6c 6f 74 74 65 72 28 |tContour|Plotter(|
|00002050| 6c 6f 6e 67 20 6e 72 52 | 6f 77 73 2c 20 6c 6f 6e |long nrR|ows, lon|
|00002060| 67 20 6e 72 43 6f 6c 73 | 2c 0d 09 09 64 6f 75 62 |g nrCols|,...doub|
|00002070| 6c 65 20 78 31 4d 69 6e | 2c 20 64 6f 75 62 6c 65 |le x1Min|, double|
|00002080| 20 78 31 4d 61 78 2c 20 | 64 6f 75 62 6c 65 20 78 | x1Max, |double x|
|00002090| 32 4d 69 6e 2c 20 64 6f | 75 62 6c 65 20 78 32 4d |2Min, do|uble x2M|
|000020a0| 61 78 29 0d 09 2f 2f 20 | 63 61 6c 6c 20 74 68 69 |ax)..// |call thi|
|000020b0| 73 20 72 6f 75 74 69 6e | 65 20 74 6f 20 69 6e 69 |s routin|e to ini|
|000020c0| 74 69 61 6c 69 7a 65 20 | 74 68 65 20 63 6f 6e 74 |tialize |the cont|
|000020d0| 6f 75 72 20 70 6c 6f 74 | 74 65 72 20 66 6f 72 20 |our plot|ter for |
|000020e0| 0d 09 2f 2f 20 77 6f 72 | 6b 69 6e 67 20 6f 6e 20 |..// wor|king on |
|000020f0| 61 20 6d 61 74 72 69 78 | 20 6f 66 20 73 69 7a 65 |a matrix| of size|
|00002100| 20 6e 72 52 6f 77 73 2c | 20 6e 72 43 6f 6c 73 0d | nrRows,| nrCols.|
|00002110| 09 2f 2f 20 78 31 4d 69 | 6e 2c 20 78 31 4d 61 78 |.// x1Mi|n, x1Max|
|00002120| 2c 20 65 74 63 2e 20 63 | 6f 72 72 65 73 70 6f 6e |, etc. c|orrespon|
|00002130| 64 20 74 6f 20 74 68 65 | 20 6e 75 6d 65 72 69 63 |d to the| numeric|
|00002140| 61 6c 20 76 61 6c 75 65 | 73 20 6f 66 20 74 68 65 |al value|s of the|
|00002150| 20 62 6f 75 6e 64 73 20 | 6f 66 20 74 68 65 20 6d | bounds |of the m|
|00002160| 61 74 72 69 78 0d 09 2f | 2f 20 72 65 74 75 72 6e |atrix../|/ return|
|00002170| 73 20 6e 69 6c 20 69 66 | 20 6e 6f 74 20 65 6e 6f |s nil if| not eno|
|00002180| 75 67 68 20 6d 65 6d 6f | 72 79 0d 09 2f 2f 20 69 |ugh memo|ry..// i|
|00002190| 66 20 49 6e 69 74 43 6f | 6e 74 6f 75 72 20 69 73 |f InitCo|ntour is|
|000021a0| 20 73 75 63 63 65 73 73 | 66 75 6c 2c 20 69 74 20 | success|ful, it |
|000021b0| 72 65 74 75 72 6e 73 20 | 74 68 65 20 61 62 6f 76 |returns |the abov|
|000021c0| 65 20 64 65 66 69 6e 65 | 64 20 73 74 72 75 63 74 |e define|d struct|
|000021d0| 2e 0d 7b 0d 0d 09 4d 79 | 50 6c 6f 74 74 69 6e 67 |..{...My|Plotting|
|000021e0| 44 61 74 61 2a 20 76 20 | 3d 20 28 4d 79 50 6c 6f |Data* v |= (MyPlo|
|000021f0| 74 74 69 6e 67 44 61 74 | 61 2a 29 4e 65 77 50 74 |ttingDat|a*)NewPt|
|00002200| 72 43 6c 65 61 72 28 73 | 69 7a 65 6f 66 28 4d 79 |rClear(s|izeof(My|
|00002210| 50 6c 6f 74 74 69 6e 67 | 44 61 74 61 29 29 3b 0d |Plotting|Data));.|
|00002220| 09 6c 6f 6e 67 09 73 69 | 7a 65 3b 0d 09 4f 53 45 |.long.si|ze;..OSE|
|00002230| 72 72 09 65 72 72 3b 0d | 09 48 61 6e 64 6c 65 09 |rr.err;.|.Handle.|
|00002240| 6d 79 4d 65 6d 6f 72 79 | 3b 0d 0d 09 69 66 20 28 |myMemory|;...if (|
|00002250| 76 3d 3d 6e 69 6c 29 20 | 72 65 74 75 72 6e 20 6e |v==nil) |return n|
|00002260| 69 6c 3b 0d 0d 09 76 2d | 3e 65 78 69 74 20 3d 20 |il;...v-|>exit = |
|00002270| 66 61 6c 73 65 3b 0d 09 | 76 2d 3e 6d 65 73 68 20 |false;..|v->mesh |
|00002280| 3d 20 6e 69 6c 3b 0d 09 | 76 2d 3e 6e 72 52 6f 77 |= nil;..|v->nrRow|
|00002290| 73 20 3d 20 6e 72 52 6f | 77 73 3b 0d 09 76 2d 3e |s = nrRo|ws;..v->|
|000022a0| 6e 72 43 6f 6c 73 20 3d | 20 6e 72 43 6f 6c 73 3b |nrCols =| nrCols;|
|000022b0| 0d 09 76 2d 3e 63 4d 69 | 6e 20 3d 20 78 31 4d 69 |..v->cMi|n = x1Mi|
|000022c0| 6e 3b 20 76 2d 3e 63 4d | 61 78 20 3d 20 78 31 4d |n; v->cM|ax = x1M|
|000022d0| 61 78 3b 0d 09 76 2d 3e | 72 4d 69 6e 20 3d 20 78 |ax;..v->|rMin = x|
|000022e0| 32 4d 69 6e 3b 20 76 2d | 3e 72 4d 61 78 20 3d 20 |2Min; v-|>rMax = |
|000022f0| 78 32 4d 61 78 3b 0d 0d | 09 76 2d 3e 72 55 6e 69 |x2Max;..|.v->rUni|
|00002300| 74 20 3d 20 28 78 32 4d | 61 78 20 2d 20 78 32 4d |t = (x2M|ax - x2M|
|00002310| 69 6e 29 2f 28 6e 72 52 | 6f 77 73 2d 31 29 3b 0d |in)/(nrR|ows-1);.|
|00002320| 09 76 2d 3e 63 55 6e 69 | 74 20 3d 20 28 78 31 4d |.v->cUni|t = (x1M|
|00002330| 61 78 20 2d 20 78 31 4d | 69 6e 29 2f 28 6e 72 43 |ax - x1M|in)/(nrC|
|00002340| 6f 6c 73 2d 31 29 3b 0d | 0d 09 73 69 7a 65 20 3d |ols-1);.|..size =|
|00002350| 20 6e 72 52 6f 77 73 2a | 6e 72 43 6f 6c 73 2a 6d | nrRows*|nrCols*m|
|00002360| 65 73 68 50 6f 69 6e 74 | 53 69 7a 65 3b 0d 09 0d |eshPoint|Size;...|
|00002370| 09 69 66 20 28 54 65 6d | 70 46 72 65 65 4d 65 6d |.if (Tem|pFreeMem|
|00002380| 28 29 20 3e 20 73 69 7a | 65 2b 35 30 30 30 30 29 |() > siz|e+50000)|
|00002390| 09 09 09 09 09 09 09 2f | 2f 20 69 66 20 74 68 65 |......./|/ if the|
|000023a0| 72 65 27 73 20 61 20 6c | 6f 74 20 6f 66 20 74 65 |re's a l|ot of te|
|000023b0| 6d 70 6f 72 61 72 79 20 | 6d 65 6d 0d 09 09 6d 79 |mporary |mem...my|
|000023c0| 4d 65 6d 6f 72 79 20 3d | 20 54 65 6d 70 4e 65 77 |Memory =| TempNew|
|000023d0| 48 61 6e 64 6c 65 28 73 | 69 7a 65 2c 20 26 65 72 |Handle(s|ize, &er|
|000023e0| 72 29 3b 09 09 09 09 2f | 2f 20 74 72 79 20 75 73 |r);..../|/ try us|
|000023f0| 69 6e 67 20 74 65 6d 70 | 6f 72 61 72 79 20 6d 65 |ing temp|orary me|
|00002400| 6d 6f 72 79 0d 09 69 66 | 20 28 6d 79 4d 65 6d 6f |mory..if| (myMemo|
|00002410| 72 79 20 3d 3d 20 6e 69 | 6c 29 0d 09 09 6d 79 4d |ry == ni|l)...myM|
|00002420| 65 6d 6f 72 79 20 3d 20 | 4e 65 77 48 61 6e 64 6c |emory = |NewHandl|
|00002430| 65 28 73 69 7a 65 29 3b | 09 09 09 09 09 09 09 2f |e(size);|......./|
|00002440| 2f 20 74 72 79 20 75 73 | 69 6e 67 20 68 65 61 70 |/ try us|ing heap|
|00002450| 20 6d 65 6d 6f 72 79 0d | 09 69 66 20 28 6d 79 4d | memory.|.if (myM|
|00002460| 65 6d 6f 72 79 20 3d 3d | 20 6e 69 6c 29 09 09 09 |emory ==| nil)...|
|00002470| 09 09 09 09 09 09 2f 2f | 20 69 66 20 61 6c 6c 6f |......//| if allo|
|00002480| 63 61 74 69 6f 6e 20 66 | 61 69 6c 65 64 0d 09 7b |cation f|ailed..{|
|00002490| 09 44 69 73 70 6f 73 65 | 50 74 72 28 28 50 74 72 |.Dispose|Ptr((Ptr|
|000024a0| 29 76 29 3b 0d 09 09 72 | 65 74 75 72 6e 20 6e 69 |)v);...r|eturn ni|
|000024b0| 6c 3b 0d 09 7d 0d 09 48 | 4c 6f 63 6b 28 6d 79 4d |l;..}..H|Lock(myM|
|000024c0| 65 6d 6f 72 79 29 3b 0d | 09 76 2d 3e 6d 65 73 68 |emory);.|.v->mesh|
|000024d0| 20 3d 20 28 76 6f 69 64 | 2a 2a 29 6d 79 4d 65 6d | = (void|**)myMem|
|000024e0| 6f 72 79 3b 0d 09 72 65 | 74 75 72 6e 20 28 50 6c |ory;..re|turn (Pl|
|000024f0| 6f 74 74 69 6e 67 44 61 | 74 61 2a 29 76 3b 0d 7d |ottingDa|ta*)v;.}|
|00002500| 2f 2f 49 6e 69 74 43 6f | 6e 74 6f 75 72 50 6c 6f |//InitCo|ntourPlo|
|00002510| 74 74 65 72 0d 0d 76 6f | 69 64 20 44 69 73 70 6f |tter..vo|id Dispo|
|00002520| 73 65 43 6f 6e 74 6f 75 | 72 50 6c 6f 74 74 65 72 |seContou|rPlotter|
|00002530| 28 50 6c 6f 74 74 69 6e | 67 44 61 74 61 2a 20 70 |(Plottin|gData* p|
|00002540| 29 0d 09 2f 2f 20 64 65 | 61 6c 6c 6f 63 61 74 65 |)..// de|allocate|
|00002550| 73 20 61 6c 6c 20 6d 65 | 6d 6f 72 79 20 61 6c 6c |s all me|mory all|
|00002560| 6f 63 61 74 65 64 20 62 | 79 20 49 6e 69 74 43 6f |ocated b|y InitCo|
|00002570| 6e 74 6f 75 72 50 6c 6f | 74 74 65 72 09 0d 09 2f |ntourPlo|tter.../|
|00002580| 2f 20 64 6f 65 73 20 6e | 6f 74 68 69 6e 67 20 69 |/ does n|othing i|
|00002590| 66 20 76 3d 3d 6e 69 6c | 0d 7b 0d 09 4d 79 50 6c |f v==nil|.{..MyPl|
|000025a0| 6f 74 74 69 6e 67 44 61 | 74 61 2a 20 76 20 3d 20 |ottingDa|ta* v = |
|000025b0| 28 4d 79 50 6c 6f 74 74 | 69 6e 67 44 61 74 61 2a |(MyPlott|ingData*|
|000025c0| 29 70 3b 0d 09 69 66 20 | 28 76 29 0d 09 7b 09 69 |)p;..if |(v)..{.i|
|000025d0| 66 20 28 76 2d 3e 6d 65 | 73 68 29 20 44 69 73 70 |f (v->me|sh) Disp|
|000025e0| 6f 73 65 48 61 6e 64 6c | 65 28 28 48 61 6e 64 6c |oseHandl|e((Handl|
|000025f0| 65 29 76 2d 3e 6d 65 73 | 68 29 3b 0d 09 09 44 69 |e)v->mes|h);...Di|
|00002600| 73 70 6f 73 65 50 74 72 | 28 28 50 74 72 29 76 29 |sposePtr|((Ptr)v)|
|00002610| 3b 0d 09 7d 0d 7d 0d 0d | 0d 42 6f 6f 6c 65 61 6e |;..}.}..|.Boolean|
|00002620| 20 43 61 6c 63 75 6c 61 | 74 65 43 6f 6e 74 6f 75 | Calcula|teContou|
|00002630| 72 4d 61 74 72 69 78 28 | 46 75 6e 63 74 69 6f 6e |rMatrix(|Function|
|00002640| 50 74 72 20 66 75 6e 63 | 2c 20 76 6f 69 64 2a 20 |Ptr func|, void* |
|00002650| 70 61 72 61 6d 2c 20 50 | 6c 6f 74 74 69 6e 67 44 |param, P|lottingD|
|00002660| 61 74 61 2a 20 70 29 0d | 09 2f 2f 20 61 66 74 65 |ata* p).|.// afte|
|00002670| 72 20 68 61 76 69 6e 67 | 20 63 61 6c 6c 65 64 20 |r having| called |
|00002680| 49 6e 69 74 43 6f 6e 74 | 6f 75 72 50 6c 6f 74 74 |InitCont|ourPlott|
|00002690| 65 72 2c 20 63 61 6c 6c | 20 74 68 69 73 20 72 6f |er, call| this ro|
|000026a0| 75 74 69 6e 65 20 74 6f | 20 69 6e 69 74 69 61 6c |utine to| initial|
|000026b0| 69 7a 65 0d 09 2f 2f 20 | 74 68 65 20 70 6c 6f 74 |ize..// |the plot|
|000026c0| 74 69 6e 67 20 6d 61 74 | 72 69 78 2e 0d 09 2f 2f |ting mat|rix...//|
|000026d0| 20 66 75 6e 63 3a 20 61 | 20 66 75 6e 63 74 69 6f | func: a| functio|
|000026e0| 6e 20 74 68 61 74 20 72 | 65 74 75 72 6e 73 20 78 |n that r|eturns x|
|000026f0| 33 20 66 6f 72 20 65 61 | 63 68 20 6d 61 74 72 69 |3 for ea|ch matri|
|00002700| 78 20 70 6f 69 6e 74 20 | 78 31 2c 20 78 32 0d 09 |x point |x1, x2..|
|00002710| 2f 2f 20 70 61 72 61 6d | 20 69 73 20 70 61 73 73 |// param| is pass|
|00002720| 65 64 20 74 6f 20 66 75 | 6e 63 0d 7b 0d 09 64 6f |ed to fu|nc.{..do|
|00002730| 75 62 6c 65 09 78 31 2c | 20 78 32 3b 0d 09 6c 6f |uble.x1,| x2;..lo|
|00002740| 6e 67 09 72 2c 20 63 3b | 0d 09 4d 79 50 6c 6f 74 |ng.r, c;|..MyPlot|
|00002750| 74 69 6e 67 44 61 74 61 | 2a 20 76 20 3d 20 28 4d |tingData|* v = (M|
|00002760| 79 50 6c 6f 74 74 69 6e | 67 44 61 74 61 2a 29 70 |yPlottin|gData*)p|
|00002770| 3b 0d 0d 09 78 32 20 3d | 20 76 2d 3e 72 4d 69 6e |;...x2 =| v->rMin|
|00002780| 3b 0d 09 66 6f 72 20 28 | 72 20 3d 20 30 3b 20 72 |;..for (|r = 0; r|
|00002790| 20 3c 20 76 2d 3e 6e 72 | 52 6f 77 73 3b 20 72 2b | < v->nr|Rows; r+|
|000027a0| 2b 29 09 09 09 2f 2f 20 | 69 6e 69 74 69 61 6c 69 |+)...// |initiali|
|000027b0| 7a 65 20 74 68 65 20 6d | 65 73 68 20 70 6f 69 6e |ze the m|esh poin|
|000027c0| 74 20 61 72 72 61 79 0d | 09 7b 0d 09 09 78 31 20 |t array.|.{...x1 |
|000027d0| 3d 20 76 2d 3e 63 4d 69 | 6e 3b 0d 09 09 66 6f 72 |= v->cMi|n;...for|
|000027e0| 20 28 63 20 3d 20 30 3b | 20 63 20 3c 20 76 2d 3e | (c = 0;| c < v->|
|000027f0| 6e 72 43 6f 6c 73 3b 20 | 63 2b 2b 29 0d 09 09 7b |nrCols; |c++)...{|
|00002800| 09 4d 65 73 68 50 6f 69 | 6e 74 2a 09 6d 65 73 68 |.MeshPoi|nt*.mesh|
|00002810| 20 3d 20 47 65 74 4d 65 | 73 68 50 6f 69 6e 74 28 | = GetMe|shPoint(|
|00002820| 72 2c 20 63 2c 20 76 29 | 3b 0d 09 09 09 6d 65 73 |r, c, v)|;....mes|
|00002830| 68 2d 3e 76 61 6c 69 64 | 20 3d 20 28 2a 66 75 6e |h->valid| = (*fun|
|00002840| 63 29 28 78 31 2c 20 78 | 32 2c 20 26 6d 65 73 68 |c)(x1, x|2, &mesh|
|00002850| 2d 3e 76 61 6c 75 65 2c | 20 70 61 72 61 6d 29 3b |->value,| param);|
|00002860| 0d 09 09 09 69 66 20 28 | 54 65 73 74 49 66 45 78 |....if (|TestIfEx|
|00002870| 69 74 28 76 29 29 20 72 | 65 74 75 72 6e 20 66 61 |it(v)) r|eturn fa|
|00002880| 6c 73 65 3b 0d 09 09 09 | 78 31 20 2b 3d 20 76 2d |lse;....|x1 += v-|
|00002890| 3e 63 55 6e 69 74 3b 0d | 09 09 7d 0d 09 09 78 32 |>cUnit;.|..}...x2|
|000028a0| 20 2b 3d 20 76 2d 3e 72 | 55 6e 69 74 3b 0d 09 7d | += v->r|Unit;..}|
|000028b0| 0d 09 72 65 74 75 72 6e | 20 74 72 75 65 3b 0d 7d |..return| true;.}|
|000028c0| 0d 0d 0d 0d 42 6f 6f 6c | 65 61 6e 20 44 72 61 77 |....Bool|ean Draw|
|000028d0| 43 6f 6e 74 6f 75 72 28 | 64 6f 75 62 6c 65 20 6c |Contour(|double l|
|000028e0| 65 76 65 6c 2c 20 50 6c | 6f 74 74 69 6e 67 44 61 |evel, Pl|ottingDa|
|000028f0| 74 61 2a 20 70 29 0d 09 | 2f 2f 20 61 66 74 65 72 |ta* p)..|// after|
|00002900| 20 68 61 76 69 6e 67 20 | 63 61 6c 6c 65 64 20 43 | having |called C|
|00002910| 61 6c 63 75 6c 61 74 65 | 43 6f 6e 74 6f 75 72 4d |alculate|ContourM|
|00002920| 61 74 72 69 78 2c 20 63 | 61 6c 6c 20 74 68 69 73 |atrix, c|all this|
|00002930| 20 72 6f 75 74 69 6e 65 | 20 66 6f 72 20 61 6c 6c | routine| for all|
|00002940| 20 6c 65 76 65 6c 73 0d | 09 2f 2f 20 28 69 2e 65 | levels.|.// (i.e|
|00002950| 2e 20 78 33 2d 76 61 6c | 75 65 73 29 20 77 68 65 |. x3-val|ues) whe|
|00002960| 72 65 20 79 6f 75 20 77 | 61 6e 74 20 74 6f 20 68 |re you w|ant to h|
|00002970| 61 76 65 20 61 20 63 6f | 6e 74 6f 75 72 20 6c 69 |ave a co|ntour li|
|00002980| 6e 65 0d 09 2f 2f 20 52 | 65 74 75 72 6e 73 20 66 |ne..// R|eturns f|
|00002990| 61 6c 73 65 20 69 66 20 | 65 78 65 63 75 74 69 6f |alse if |executio|
|000029a0| 6e 20 6d 75 73 74 20 62 | 65 20 73 74 6f 70 70 65 |n must b|e stoppe|
|000029b0| 64 0d 0d 7b 0d 09 6c 6f | 6e 67 20 72 2c 20 63 3b |d..{..lo|ng r, c;|
|000029c0| 0d 09 4d 79 50 6c 6f 74 | 74 69 6e 67 44 61 74 61 |..MyPlot|tingData|
|000029d0| 2a 20 76 20 3d 20 28 4d | 79 50 6c 6f 74 74 69 6e |* v = (M|yPlottin|
|000029e0| 67 44 61 74 61 2a 29 70 | 3b 0d 0d 09 2f 2f 20 66 |gData*)p|;...// f|
|000029f0| 69 72 73 74 20 63 68 65 | 63 6b 20 77 68 65 72 65 |irst che|ck where|
|00002a00| 20 77 65 20 68 61 76 65 | 20 63 72 6f 73 73 69 6e | we have| crossin|
|00002a10| 67 20 70 6f 69 6e 74 73 | 20 69 6e 20 74 68 65 20 |g points| in the |
|00002a20| 6d 61 74 72 69 78 3a 0d | 09 66 6f 72 20 28 72 20 |matrix:.|.for (r |
|00002a30| 3d 20 30 3b 20 72 20 3c | 20 76 2d 3e 6e 72 52 6f |= 0; r <| v->nrRo|
|00002a40| 77 73 20 26 26 20 21 54 | 65 73 74 49 66 45 78 69 |ws && !T|estIfExi|
|00002a50| 74 28 76 29 3b 20 72 2b | 2b 29 0d 09 7b 0d 09 09 |t(v); r+|+)..{...|
|00002a60| 66 6f 72 20 28 63 20 3d | 20 30 3b 20 63 20 3c 20 |for (c =| 0; c < |
|00002a70| 76 2d 3e 6e 72 43 6f 6c | 73 20 26 26 20 21 54 65 |v->nrCol|s && !Te|
|00002a80| 73 74 49 66 45 78 69 74 | 28 76 29 3b 20 63 2b 2b |stIfExit|(v); c++|
|00002a90| 29 0d 09 09 7b 0d 09 09 | 09 4d 65 73 68 50 6f 69 |)...{...|.MeshPoi|
|00002aa0| 6e 74 2a 09 6d 65 73 68 | 20 3d 20 47 65 74 4d 65 |nt*.mesh| = GetMe|
|00002ab0| 73 68 50 6f 69 6e 74 28 | 72 2c 20 63 2c 20 76 29 |shPoint(|r, c, v)|
|00002ac0| 3b 0d 09 09 09 69 66 20 | 28 6d 65 73 68 2d 3e 76 |;....if |(mesh->v|
|00002ad0| 61 6c 69 64 29 09 09 09 | 09 09 09 09 09 09 2f 2f |alid)...|......//|
|00002ae0| 20 69 66 20 74 68 65 72 | 65 27 73 20 76 61 6c 69 | if ther|e's vali|
|00002af0| 64 20 64 61 74 61 20 69 | 6e 20 74 68 69 73 20 63 |d data i|n this c|
|00002b00| 65 6c 6c 0d 09 09 09 7b | 0d 09 09 09 09 42 6f 6f |ell....{|.....Boo|
|00002b10| 6c 65 61 6e 09 74 68 69 | 73 49 73 42 65 6c 6f 77 |lean.thi|sIsBelow|
|00002b20| 4c 65 76 65 6c 3b 0d 09 | 09 09 09 6d 65 73 68 2d |Level;..|...mesh-|
|00002b30| 3e 61 74 52 69 67 68 74 | 20 3d 20 6d 65 73 68 2d |>atRight| = mesh-|
|00002b40| 3e 61 74 42 6f 74 74 6f | 6d 20 3d 20 66 61 6c 73 |>atBotto|m = fals|
|00002b50| 65 3b 0d 09 09 09 09 74 | 68 69 73 49 73 42 65 6c |e;.....t|hisIsBel|
|00002b60| 6f 77 4c 65 76 65 6c 20 | 3d 20 6d 65 73 68 2d 3e |owLevel |= mesh->|
|00002b70| 76 61 6c 75 65 20 3c 20 | 6c 65 76 65 6c 3b 0d 09 |value < |level;..|
|00002b80| 09 09 09 69 66 20 28 63 | 20 3c 20 76 2d 3e 6e 72 |...if (c| < v->nr|
|00002b90| 43 6f 6c 73 2d 31 29 0d | 09 09 09 09 7b 09 4d 65 |Cols-1).|....{.Me|
|00002ba0| 73 68 50 6f 69 6e 74 2a | 20 72 69 67 68 74 20 3d |shPoint*| right =|
|00002bb0| 20 47 65 74 4d 65 73 68 | 50 6f 69 6e 74 28 72 2c | GetMesh|Point(r,|
|00002bc0| 20 63 2b 31 2c 20 76 29 | 3b 09 09 2f 2f 20 74 68 | c+1, v)|;..// th|
|00002bd0| 69 73 20 69 73 20 6e 6f | 74 20 76 65 72 79 20 65 |is is no|t very e|
|00002be0| 66 66 69 63 69 65 6e 74 | 2c 20 6d 69 67 68 74 20 |fficient|, might |
|00002bf0| 75 73 65 20 62 75 66 66 | 65 72 69 6e 67 20 68 65 |use buff|ering he|
|00002c00| 72 65 0d 09 09 09 09 09 | 69 66 20 28 72 69 67 68 |re......|if (righ|
|00002c10| 74 2d 3e 76 61 6c 69 64 | 20 26 26 20 28 74 68 69 |t->valid| && (thi|
|00002c20| 73 49 73 42 65 6c 6f 77 | 4c 65 76 65 6c 20 21 3d |sIsBelow|Level !=|
|00002c30| 20 28 72 69 67 68 74 2d | 3e 76 61 6c 75 65 20 3c | (right-|>value <|
|00002c40| 20 6c 65 76 65 6c 29 29 | 29 0d 09 09 09 09 09 09 | level))|).......|
|00002c50| 6d 65 73 68 2d 3e 61 74 | 52 69 67 68 74 20 3d 20 |mesh->at|Right = |
|00002c60| 74 72 75 65 3b 0d 09 09 | 09 09 7d 0d 09 09 09 09 |true;...|..}.....|
|00002c70| 69 66 20 28 72 20 3c 20 | 76 2d 3e 6e 72 52 6f 77 |if (r < |v->nrRow|
|00002c80| 73 2d 31 29 0d 09 09 09 | 09 7b 09 4d 65 73 68 50 |s-1)....|.{.MeshP|
|00002c90| 6f 69 6e 74 2a 20 62 6f | 74 74 6f 6d 20 3d 20 47 |oint* bo|ttom = G|
|00002ca0| 65 74 4d 65 73 68 50 6f | 69 6e 74 28 72 2b 31 2c |etMeshPo|int(r+1,|
|00002cb0| 20 63 2c 20 76 29 3b 0d | 09 09 09 09 09 69 66 20 | c, v);.|.....if |
|00002cc0| 28 62 6f 74 74 6f 6d 2d | 3e 76 61 6c 69 64 20 26 |(bottom-|>valid &|
|00002cd0| 26 20 28 74 68 69 73 49 | 73 42 65 6c 6f 77 4c 65 |& (thisI|sBelowLe|
|00002ce0| 76 65 6c 20 21 3d 20 28 | 62 6f 74 74 6f 6d 2d 3e |vel != (|bottom->|
|00002cf0| 76 61 6c 75 65 20 3c 20 | 6c 65 76 65 6c 29 29 29 |value < |level)))|
|00002d00| 0d 09 09 09 09 09 09 6d | 65 73 68 2d 3e 61 74 42 |.......m|esh->atB|
|00002d10| 6f 74 74 6f 6d 20 3d 20 | 74 72 75 65 3b 0d 09 09 |ottom = |true;...|
|00002d20| 09 09 7d 0d 09 09 09 7d | 0d 09 09 7d 0d 09 7d 2f |..}....}|...}..}/|
|00002d30| 2f 66 6f 72 20 72 0d 0d | 09 2f 2f 20 6e 6f 77 20 |/for r..|.// now |
|00002d40| 66 6f 6c 6c 6f 77 20 61 | 6c 6c 20 63 72 6f 73 73 |follow a|ll cross|
|00002d50| 69 6e 67 20 70 6f 69 6e | 74 73 20 74 6f 20 63 72 |ing poin|ts to cr|
|00002d60| 65 61 74 65 20 74 68 65 | 20 63 6f 6e 74 6f 75 72 |eate the| contour|
|00002d70| 73 3a 0d 09 66 6f 72 20 | 28 72 20 3d 20 30 3b 20 |s:..for |(r = 0; |
|00002d80| 72 20 3c 20 76 2d 3e 6e | 72 52 6f 77 73 2d 31 20 |r < v->n|rRows-1 |
|00002d90| 26 26 20 21 54 65 73 74 | 49 66 45 78 69 74 28 76 |&& !Test|IfExit(v|
|00002da0| 29 3b 20 72 2b 2b 29 0d | 09 7b 09 66 6f 72 20 28 |); r++).|.{.for (|
|00002db0| 63 20 3d 20 30 3b 20 63 | 20 3c 20 76 2d 3e 6e 72 |c = 0; c| < v->nr|
|00002dc0| 43 6f 6c 73 2d 31 20 26 | 26 20 21 54 65 73 74 49 |Cols-1 &|& !TestI|
|00002dd0| 66 45 78 69 74 28 76 29 | 3b 20 63 2b 2b 29 0d 09 |fExit(v)|; c++)..|
|00002de0| 09 7b 0d 09 09 09 4d 65 | 73 68 50 6f 69 6e 74 2a |.{....Me|shPoint*|
|00002df0| 09 69 6e 66 6f 20 3d 20 | 47 65 74 4d 65 73 68 50 |.info = |GetMeshP|
|00002e00| 6f 69 6e 74 28 72 2c 20 | 63 2c 20 76 29 3b 0d 09 |oint(r, |c, v);..|
|00002e10| 09 09 69 66 20 28 69 6e | 66 6f 2d 3e 61 74 52 69 |..if (in|fo->atRi|
|00002e20| 67 68 74 29 20 44 72 61 | 77 4f 6e 65 43 6f 6e 74 |ght) Dra|wOneCont|
|00002e30| 6f 75 72 28 6c 65 76 65 | 6c 2c 20 72 2c 20 63 2c |our(leve|l, r, c,|
|00002e40| 20 76 29 3b 0d 09 09 09 | 69 66 20 28 69 6e 66 6f | v);....|if (info|
|00002e50| 2d 3e 61 74 42 6f 74 74 | 6f 6d 29 20 44 72 61 77 |->atBott|om) Draw|
|00002e60| 4f 6e 65 43 6f 6e 74 6f | 75 72 28 6c 65 76 65 6c |OneConto|ur(level|
|00002e70| 2c 20 72 2c 20 63 2c 20 | 76 29 3b 0d 09 09 7d 20 |, r, c, |v);...} |
|00002e80| 2f 2a 66 6f 72 20 63 2a | 2f 0d 09 7d 20 2f 2a 66 |/*for c*|/..} /*f|
|00002e90| 6f 72 20 72 2a 2f 0d 09 | 72 65 74 75 72 6e 20 21 |or r*/..|return !|
|00002ea0| 54 65 73 74 49 66 45 78 | 69 74 28 76 29 3b 0d 7d |TestIfEx|it(v);.}|
|00002eb0| 0d 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ec0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ed0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ee0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ef0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f00| 00 00 01 00 00 00 01 54 | 00 00 00 54 00 00 00 46 |.......T|...T...F|
|00002f10| 33 00 00 66 66 ff ff 00 | 00 00 00 66 66 cc cc ff |3..ff...|...ff...|
|00002f20| ff 00 00 66 66 cc cc cc | cc 00 00 66 66 cc cc 99 |...ff...|...ff...|
|00002f30| 10 43 6f 6e 74 6f 75 72 | 50 6c 6f 74 74 65 72 2e |.Contour|Plotter.|
|00002f40| 63 74 02 00 00 00 54 45 | 58 54 43 57 49 45 01 00 |ct....TE|XTCWIE..|
|00002f50| 00 00 54 45 58 54 43 57 | 49 45 01 00 00 00 01 80 |..TEXTCW|IE......|
|00002f60| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002f70| 00 00 ad 6e cc 8a 00 00 | 2e 31 00 00 01 9a 66 ff |...n....|.1....f.|
|00002f80| ff 00 00 66 66 66 66 cc | cc 00 00 66 66 66 66 99 |...ffff.|...ffff.|
|00002f90| 99 00 00 66 66 66 66 66 | 66 00 00 66 66 66 66 33 |...fffff|f..ffff3|
|00002fa0| 33 00 00 66 66 66 66 00 | 00 00 00 66 66 33 33 ff |3..ffff.|...ff33.|
|00002fb0| ff 00 00 66 66 33 33 cc | cc 00 00 66 66 33 33 99 |...ff33.|...ff33.|
|00002fc0| 99 00 00 66 66 33 33 66 | 66 00 00 66 66 33 33 33 |...ff33f|f..ff333|
|00002fd0| 33 00 00 66 66 33 33 00 | 00 00 00 66 66 00 00 ff |3..ff33.|...ff...|
|00002fe0| ff 00 00 66 66 00 00 cc | cc 00 00 66 66 00 00 99 |...ff...|...ff...|
|00002ff0| 99 00 00 66 66 00 00 66 | 66 00 00 66 66 00 00 33 |...ff..f|f..ff..3|
|00003000| 00 00 00 48 00 09 4d 6f | 6e 61 63 6f 00 00 00 00 |...H..Mo|naco....|
|00003010| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003020| 00 00 00 00 00 00 00 03 | 00 04 00 57 00 03 03 62 |........|...W...b|
|00003030| 02 7d 00 57 00 03 03 62 | 02 7d ad 9f fd 80 00 00 |.}.W...b|.}......|
|00003040| 2c ad 00 00 2c ad 00 00 | 26 8f 00 00 00 00 00 04 |,...,...|&.......|
|00003050| 00 01 00 01 00 00 01 00 | 00 00 01 54 00 00 00 54 |........|...T...T|
|00003060| 00 00 00 46 02 94 a5 c0 | 24 5c 00 00 00 1c 00 46 |...F....|$\.....F|
|00003070| 00 01 4d 50 53 52 00 00 | 00 12 4d 57 42 42 00 00 |..MPSR..|..MWBB..|
|00003080| 00 1e 03 ed ff ff 00 00 | 00 00 00 00 00 00 03 f0 |........|........|
|00003090| ff ff 00 00 00 4c 00 00 | 00 00 00 00 00 00 00 00 |.....L..|........|
|000030a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000030b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000030c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000030d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000030e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000030f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
+--------+-------------------------+-------------------------+--------+--------+