Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] View of /branches/vis12/test/square-detect/square_score.diderot
ViewVC logotype

View of /branches/vis12/test/square-detect/square_score.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1927 - (download) (annotate)
Sat Jun 23 18:09:16 2012 UTC (7 years, 3 months ago) by jhr
File size: 20333 byte(s)
  converting to use "image" instead of "load" for image nrrd loading
field#2(2)[] T = bspln3 ⊛ image("square-template.nrrd");
field#2(2)[] F = bspln3 ⊛ image("square-rotate.nrrd");

vec2 v00 = [0.1, 0.1];
vec2 v01 = [0.1, 0.2];
vec2 v02 = [0.1, 0.3];
vec2 v03 = [0.1, 0.4];
vec2 v04 = [0.1, 0.5];
vec2 v05 = [0.1, -0.1];
vec2 v06 = [0.1, -0.2];
vec2 v07 = [0.1, -0.3];
vec2 v08 = [0.1, -0.4];
vec2 v09 = [0.1, -0.5];
vec2 v10 = [0.2, 0.1];
vec2 v11 = [0.2, 0.2];
vec2 v12 = [0.2, 0.3];
vec2 v13 = [0.2, 0.4];
vec2 v14 = [0.2, 0.5];
vec2 v15 = [0.2, -0.1];
vec2 v16 = [0.2, -0.2];
vec2 v17 = [0.2, -0.3];
vec2 v18 = [0.2, -0.4];
vec2 v19 = [0.2, -0.5];
vec2 v20 = [0.3, 0.1];
vec2 v21 = [0.3, 0.2];
vec2 v22 = [0.3, 0.3];
vec2 v23 = [0.3, 0.4];
vec2 v24 = [0.3, 0.5];
vec2 v25 = [0.3, -0.1];
vec2 v26 = [0.3, -0.2];
vec2 v27 = [0.3, -0.3];
vec2 v28 = [0.3, -0.4];
vec2 v29 = [0.3, -0.5];
vec2 v30 = [0.4, 0.1];
vec2 v31 = [0.4, 0.2];
vec2 v32 = [0.4, 0.3];
vec2 v33 = [0.4, 0.4];
vec2 v34 = [0.4, 0.5];
vec2 v35 = [0.4, -0.1];
vec2 v36 = [0.4, -0.2];
vec2 v37 = [0.4, -0.3];
vec2 v38 = [0.4, -0.4];
vec2 v39 = [0.4, -0.5];
vec2 v40 = [0.5, 0.1];
vec2 v41 = [0.5, 0.2];
vec2 v42 = [0.5, 0.3];
vec2 v43 = [0.5, 0.4];
vec2 v44 = [0.5, 0.5];
vec2 v45 = [0.5, -0.1];
vec2 v46 = [0.5, -0.2];
vec2 v47 = [0.5, -0.3];
vec2 v48 = [0.5, -0.4];
vec2 v49 = [0.5, -0.5];
vec2 v50 = [-0.1, 0.1];
vec2 v51 = [-0.1, 0.2];
vec2 v52 = [-0.1, 0.3];
vec2 v53 = [-0.1, 0.4];
vec2 v54 = [-0.1, 0.5];
vec2 v55 = [-0.1, -0.1];
vec2 v56 = [-0.1, -0.2];
vec2 v57 = [-0.1, -0.3];
vec2 v58 = [-0.1, -0.4];
vec2 v59 = [-0.1, -0.5];
vec2 v60 = [-0.2, 0.1];
vec2 v61 = [-0.2, 0.2];
vec2 v62 = [-0.2, 0.3];
vec2 v63 = [-0.2, 0.4];
vec2 v64 = [-0.2, 0.5];
vec2 v65 = [-0.2, -0.1];
vec2 v66 = [-0.2, -0.2];
vec2 v67 = [-0.2, -0.3];
vec2 v68 = [-0.2, -0.4];
vec2 v69 = [-0.2, -0.5];
vec2 v70 = [-0.3, 0.1];
vec2 v71 = [-0.3, 0.2];
vec2 v72 = [-0.3, 0.3];
vec2 v73 = [-0.3, 0.4];
vec2 v74 = [-0.3, 0.5];
vec2 v75 = [-0.3, -0.1];
vec2 v76 = [-0.3, -0.2];
vec2 v77 = [-0.3, -0.3];
vec2 v78 = [-0.3, -0.4];
vec2 v79 = [-0.3, -0.5];
vec2 v80 = [-0.4, 0.1];
vec2 v81 = [-0.4, 0.2];
vec2 v82 = [-0.4, 0.3];
vec2 v83 = [-0.4, 0.4];
vec2 v84 = [-0.4, 0.5];
vec2 v85 = [-0.4, -0.1];
vec2 v86 = [-0.4, -0.2];
vec2 v87 = [-0.4, -0.3];
vec2 v88 = [-0.4, -0.4];
vec2 v89 = [-0.4, -0.5];
vec2 v90 = [-0.5, 0.1];
vec2 v91 = [-0.5, 0.2];
vec2 v92 = [-0.5, 0.3];
vec2 v93 = [-0.5, 0.4];
vec2 v94 = [-0.5, 0.5];
vec2 v95 = [-0.5, -0.1];
vec2 v96 = [-0.5, -0.2];
vec2 v97 = [-0.5, -0.3];
vec2 v98 = [-0.5, -0.4];
vec2 v99 = [-0.5, -0.5];

strand R (int i, int j, int k)
{
	vec3 t = [real(i) / 60.0 - 0.5, real(j) / 60.0 - 0.5, real(k) / 60.0 - 0.5];
    output real score = 0.0;

	vec2 r0 = [0.0, 0.0];
	vec2 r1 = [0.0, 0.0];

	vec2 pos = [0.0, 0.0];


	update
	{
		r0 = [1.0 - 2.0 * t[2] * t[2], -2.0*sqrt(1.0 - t[2] * t[2]) * t[2]];
		r1 = [2.0 * sqrt(1.0 - t[2] * t[2]) * t[2], 1.0 - 2.0 * t[2] * t[2]];


        // v00
		pos = [r0 • v00, r1 • v00] + [t[0], t[1]];
		if(inside(pos,F) && inside (v00, T))
        {
        score += (F(pos) - T(v00)) * (F(pos) - T(v00));
        }
        // v01
		pos = [r0 • v01, r1 • v01] + [t[0], t[1]];
		if(inside(pos,F) && inside (v01, T))
        {
		score += (F(pos) - T(v01)) * (F(pos) - T(v01));
        }
        // v02
		pos = [r0 • v02, r1 • v02] + [t[0], t[1]];
		if(inside(pos,F) && inside (v02, T))
        {
		score += (F(pos) - T(v02)) * (F(pos) - T(v02));
        }
        // v03
		pos = [r0 • v03, r1 • v03] + [t[0], t[1]];
		if(inside(pos,F) && inside (v03, T))
        {
		score += (F(pos) - T(v03)) * (F(pos) - T(v03));
        }
        // v04
		pos = [r0 • v04, r1 • v04] + [t[0], t[1]];
		if(inside(pos,F) && inside (v04, T))
        {
		score += (F(pos) - T(v04)) * (F(pos) - T(v04));
        }
        // v05
		pos = [r0 • v05, r1 • v05] + [t[0], t[1]];
		if(inside(pos,F) && inside (v05, T))
        {
		score += (F(pos) - T(v05)) * (F(pos) - T(v05));
        }
        // v06
		pos = [r0 • v06, r1 • v06] + [t[0], t[1]];
		if(inside(pos,F) && inside (v06, T))
        {
		score += (F(pos) - T(v06)) * (F(pos) - T(v06));
        }
        // v07
		pos = [r0 • v07, r1 • v07] + [t[0], t[1]];
		if(inside(pos,F) && inside (v07, T))
        {
		score += (F(pos) - T(v07)) * (F(pos) - T(v07));
        }
        // v08
		pos = [r0 • v08, r1 • v08] + [t[0], t[1]];
		if(inside(pos,F) && inside (v08, T))
        {
		score += (F(pos) - T(v08)) * (F(pos) - T(v08));
        }
        // v09
		pos = [r0 • v09, r1 • v09] + [t[0], t[1]];
		if(inside(pos,F) && inside (v09, T))
        {
		score += (F(pos) - T(v09)) * (F(pos) - T(v09));
        }
        // v10
		pos = [r0 • v10, r1 • v10] + [t[0], t[1]];
		if(inside(pos,F) && inside (v10, T))
        {
		score += (F(pos) - T(v10)) * (F(pos) - T(v10));
        }
        // v11
		pos = [r0 • v11, r1 • v11] + [t[0], t[1]];
		if(inside(pos,F) && inside (v11, T))
        {
		score += (F(pos) - T(v11)) * (F(pos) - T(v11));
        }
        // v12
		pos = [r0 • v12, r1 • v12] + [t[0], t[1]];
		if(inside(pos,F) && inside (v12, T))
        {
		score += (F(pos) - T(v12)) * (F(pos) - T(v12));
        }
        // v13
		pos = [r0 • v13, r1 • v13] + [t[0], t[1]];
		if(inside(pos,F) && inside (v13, T))
        {
		score += (F(pos) - T(v13)) * (F(pos) - T(v13));
        }
        // v14
		pos = [r0 • v14, r1 • v14] + [t[0], t[1]];
		if(inside(pos,F) && inside (v14, T))
        {
		score += (F(pos) - T(v14)) * (F(pos) - T(v14));
        }
        // v15
		pos = [r0 • v15, r1 • v15] + [t[0], t[1]];
		if(inside(pos,F) && inside (v15, T))
        {
		score += (F(pos) - T(v15)) * (F(pos) - T(v15));
        }
        // v16
		pos = [r0 • v16, r1 • v16] + [t[0], t[1]];
		if(inside(pos,F) && inside (v16, T))
        {
		score += (F(pos) - T(v16)) * (F(pos) - T(v16));
        }
        // v17
		pos = [r0 • v17, r1 • v17] + [t[0], t[1]];
		if(inside(pos,F) && inside (v17, T))
        {
		score += (F(pos) - T(v17)) * (F(pos) - T(v17));
        }
        // v18
		pos = [r0 • v18, r1 • v18] + [t[0], t[1]];
		if(inside(pos,F) && inside (v18, T))
        {
		score += (F(pos) - T(v18)) * (F(pos) - T(v18));
        }
        // v19
		pos = [r0 • v19, r1 • v19] + [t[0], t[1]];
		if(inside(pos,F) && inside (v19, T))
        {
		score += (F(pos) - T(v19)) * (F(pos) - T(v19));
        }
        // v20
		pos = [r0 • v20, r1 • v20] + [t[0], t[1]];
		if(inside(pos,F) && inside (v20, T))
        {
		score += (F(pos) - T(v20)) * (F(pos) - T(v20));
        }
        // v21
		pos = [r0 • v21, r1 • v21] + [t[0], t[1]];
		if(inside(pos,F) && inside (v21, T))
        {
		score += (F(pos) - T(v21)) * (F(pos) - T(v21));
        }
        // v22
		pos = [r0 • v22, r1 • v22] + [t[0], t[1]];
		if(inside(pos,F) && inside (v22, T))
        {
		score += (F(pos) - T(v22)) * (F(pos) - T(v22));
        }
        // v23
		pos = [r0 • v23, r1 • v23] + [t[0], t[1]];
		if(inside(pos,F) && inside (v23, T))
        {
		score += (F(pos) - T(v23)) * (F(pos) - T(v23));
        }
        // v24
		pos = [r0 • v24, r1 • v24] + [t[0], t[1]];
		if(inside(pos,F) && inside (v24, T))
        {
		score += (F(pos) - T(v24)) * (F(pos) - T(v24));
        }
        // v25
		pos = [r0 • v25, r1 • v25] + [t[0], t[1]];
		if(inside(pos,F) && inside (v25, T))
        {
		score += (F(pos) - T(v25)) * (F(pos) - T(v25));
        }
        // v26
		pos = [r0 • v26, r1 • v26] + [t[0], t[1]];
		if(inside(pos,F) && inside (v26, T))
        {
		score += (F(pos) - T(v26)) * (F(pos) - T(v26));
        }
        // v27
		pos = [r0 • v27, r1 • v27] + [t[0], t[1]];
		if(inside(pos,F) && inside (v27, T))
        {
		score += (F(pos) - T(v27)) * (F(pos) - T(v27));
        }
        // v28
		pos = [r0 • v28, r1 • v28] + [t[0], t[1]];
		if(inside(pos,F) && inside (v28, T))
        {
		score += (F(pos) - T(v28)) * (F(pos) - T(v28));
        }
        // v29
		pos = [r0 • v29, r1 • v29] + [t[0], t[1]];
		if(inside(pos,F) && inside (v29, T))
        {
		score += (F(pos) - T(v29)) * (F(pos) - T(v29));
        }
        // v30
		pos = [r0 • v30, r1 • v30] + [t[0], t[1]];
		if(inside(pos,F) && inside (v30, T))
        {
		score += (F(pos) - T(v30)) * (F(pos) - T(v30));
        }
        // v31
		pos = [r0 • v31, r1 • v31] + [t[0], t[1]];
		if(inside(pos,F) && inside (v31, T))
        {
		score += (F(pos) - T(v31)) * (F(pos) - T(v31));
        }
        // v32
		pos = [r0 • v32, r1 • v32] + [t[0], t[1]];
		if(inside(pos,F) && inside (v32, T))
        {
		score += (F(pos) - T(v32)) * (F(pos) - T(v32));
        }
        // v33
		pos = [r0 • v33, r1 • v33] + [t[0], t[1]];
		if(inside(pos,F) && inside (v33, T))
        {
		score += (F(pos) - T(v33)) * (F(pos) - T(v33));
        }
        // v34
		pos = [r0 • v34, r1 • v34] + [t[0], t[1]];
		if(inside(pos,F) && inside (v34, T))
        {
		score += (F(pos) - T(v34)) * (F(pos) - T(v34));
        }
        // v35
		pos = [r0 • v35, r1 • v35] + [t[0], t[1]];
		if(inside(pos,F) && inside (v35, T))
        {
		score += (F(pos) - T(v35)) * (F(pos) - T(v35));
        }
        // v36
		pos = [r0 • v36, r1 • v36] + [t[0], t[1]];
		if(inside(pos,F) && inside (v36, T))
        {
		score += (F(pos) - T(v36)) * (F(pos) - T(v36));
        }
        // v37
		pos = [r0 • v37, r1 • v37] + [t[0], t[1]];
		if(inside(pos,F) && inside (v37, T))
        {
		score += (F(pos) - T(v37)) * (F(pos) - T(v37));
        }
        // v38
		pos = [r0 • v38, r1 • v38] + [t[0], t[1]];
		if(inside(pos,F) && inside (v38, T))
        {
		score += (F(pos) - T(v38)) * (F(pos) - T(v38));
        }
        // v39
		pos = [r0 • v39, r1 • v39] + [t[0], t[1]];
		if(inside(pos,F) && inside (v39, T))
        {
		score += (F(pos) - T(v39)) * (F(pos) - T(v39));
        }
        // v40
		pos = [r0 • v40, r1 • v40] + [t[0], t[1]];
		if(inside(pos,F) && inside (v40, T))
        {
		score += (F(pos) - T(v40)) * (F(pos) - T(v40));
        }
        // v41
		pos = [r0 • v41, r1 • v41] + [t[0], t[1]];
		if(inside(pos,F) && inside (v41, T))
        {
		score += (F(pos) - T(v41)) * (F(pos) - T(v41));
        }
        // v42
		pos = [r0 • v42, r1 • v42] + [t[0], t[1]];
		if(inside(pos,F) && inside (v42, T))
        {
		score += (F(pos) - T(v42)) * (F(pos) - T(v42));
        }
        // v43
		pos = [r0 • v43, r1 • v43] + [t[0], t[1]];
		if(inside(pos,F) && inside (v43, T))
        {
		score += (F(pos) - T(v43)) * (F(pos) - T(v43));
        }
        // v44
		pos = [r0 • v44, r1 • v44] + [t[0], t[1]];
		if(inside(pos,F) && inside (v44, T))
        {
		score += (F(pos) - T(v44)) * (F(pos) - T(v44));
        }
        // v45
		pos = [r0 • v45, r1 • v45] + [t[0], t[1]];
		if(inside(pos,F) && inside (v45, T))
        {
		score += (F(pos) - T(v45)) * (F(pos) - T(v45));
        }
        // v46
		pos = [r0 • v46, r1 • v46] + [t[0], t[1]];
		if(inside(pos,F) && inside (v46, T))
        {
		score += (F(pos) - T(v46)) * (F(pos) - T(v46));
        }
        // v47
		pos = [r0 • v47, r1 • v47] + [t[0], t[1]];
		if(inside(pos,F) && inside (v47, T))
        {
		score += (F(pos) - T(v47)) * (F(pos) - T(v47));
        }
        // v48
		pos = [r0 • v48, r1 • v48] + [t[0], t[1]];
		if(inside(pos,F) && inside (v48, T))
        {
		score += (F(pos) - T(v48)) * (F(pos) - T(v48));
        }
        // v49
		pos = [r0 • v49, r1 • v49] + [t[0], t[1]];
		if(inside(pos,F) && inside (v49, T))
        {
		score += (F(pos) - T(v49)) * (F(pos) - T(v49));
        }
        // v50
		pos = [r0 • v50, r1 • v50] + [t[0], t[1]];
		if(inside(pos,F) && inside (v50, T))
        {
		score += (F(pos) - T(v50)) * (F(pos) - T(v50));
        }
        // v51
		pos = [r0 • v51, r1 • v51] + [t[0], t[1]];
		if(inside(pos,F) && inside (v51, T))
        {
		score += (F(pos) - T(v51)) * (F(pos) - T(v51));
        }
        // v52
		pos = [r0 • v52, r1 • v52] + [t[0], t[1]];
		if(inside(pos,F) && inside (v52, T))
        {
		score += (F(pos) - T(v52)) * (F(pos) - T(v52));
        }
        // v53
		pos = [r0 • v53, r1 • v53] + [t[0], t[1]];
		if(inside(pos,F) && inside (v53, T))
        {
		score += (F(pos) - T(v53)) * (F(pos) - T(v53));
        }
        // v54
		pos = [r0 • v54, r1 • v54] + [t[0], t[1]];
		if(inside(pos,F) && inside (v54, T))
        {
		score += (F(pos) - T(v54)) * (F(pos) - T(v54));
        }
        // v55
		pos = [r0 • v55, r1 • v55] + [t[0], t[1]];
		if(inside(pos,F) && inside (v55, T))
        {
		score += (F(pos) - T(v55)) * (F(pos) - T(v55));
        }
        // v56
		pos = [r0 • v56, r1 • v56] + [t[0], t[1]];
		if(inside(pos,F) && inside (v56, T))
        {
		score += (F(pos) - T(v56)) * (F(pos) - T(v56));
        }
        // v57
		pos = [r0 • v57, r1 • v57] + [t[0], t[1]];
		if(inside(pos,F) && inside (v57, T))
        {
		score += (F(pos) - T(v57)) * (F(pos) - T(v57));
        }
        // v58
		pos = [r0 • v58, r1 • v58] + [t[0], t[1]];
		if(inside(pos,F) && inside (v58, T))
        {
		score += (F(pos) - T(v58)) * (F(pos) - T(v58));
        }
        // v59
		pos = [r0 • v59, r1 • v59] + [t[0], t[1]];
		if(inside(pos,F) && inside (v59, T))
        {
		score += (F(pos) - T(v59)) * (F(pos) - T(v59));
        }
        // v60
		pos = [r0 • v60, r1 • v60] + [t[0], t[1]];
		if(inside(pos,F) && inside (v60, T))
        {
		score += (F(pos) - T(v60)) * (F(pos) - T(v60));
        }
        // v61
		pos = [r0 • v61, r1 • v61] + [t[0], t[1]];
		if(inside(pos,F) && inside (v61, T))
        {
		score += (F(pos) - T(v61)) * (F(pos) - T(v61));
        }
        // v62
		pos = [r0 • v62, r1 • v62] + [t[0], t[1]];
		if(inside(pos,F) && inside (v62, T))
        {
		score += (F(pos) - T(v62)) * (F(pos) - T(v62));
        }
        // v63
		pos = [r0 • v63, r1 • v63] + [t[0], t[1]];
		if(inside(pos,F) && inside (v63, T))
        {
		score += (F(pos) - T(v63)) * (F(pos) - T(v63));
        }
        // v64
		pos = [r0 • v64, r1 • v64] + [t[0], t[1]];
		if(inside(pos,F) && inside (v64, T))
        {
		score += (F(pos) - T(v64)) * (F(pos) - T(v64));
        }
        // v65
		pos = [r0 • v65, r1 • v65] + [t[0], t[1]];
		if(inside(pos,F) && inside (v65, T))
        {
		score += (F(pos) - T(v65)) * (F(pos) - T(v65));
        }
        // v66
		pos = [r0 • v66, r1 • v66] + [t[0], t[1]];
		if(inside(pos,F) && inside (v66, T))
        {
		score += (F(pos) - T(v66)) * (F(pos) - T(v66));
        }
        // v67
		pos = [r0 • v67, r1 • v67] + [t[0], t[1]];
		if(inside(pos,F) && inside (v67, T))
        {
		score += (F(pos) - T(v67)) * (F(pos) - T(v67));
        }
        // v68
		pos = [r0 • v68, r1 • v68] + [t[0], t[1]];
		if(inside(pos,F) && inside (v68, T))
        {
		score += (F(pos) - T(v68)) * (F(pos) - T(v68));
        }
        // v69
		pos = [r0 • v69, r1 • v69] + [t[0], t[1]];
		if(inside(pos,F) && inside (v69, T))
        {
		score += (F(pos) - T(v69)) * (F(pos) - T(v69));
        }
        // v70
		pos = [r0 • v70, r1 • v70] + [t[0], t[1]];
		if(inside(pos,F) && inside (v70, T))
        {
		score += (F(pos) - T(v70)) * (F(pos) - T(v70));
        }
        // v71
		pos = [r0 • v71, r1 • v71] + [t[0], t[1]];
		if(inside(pos,F) && inside (v71, T))
        {
		score += (F(pos) - T(v71)) * (F(pos) - T(v71));
        }
        // v72
		pos = [r0 • v72, r1 • v72] + [t[0], t[1]];
		if(inside(pos,F) && inside (v72, T))
        {
		score += (F(pos) - T(v72)) * (F(pos) - T(v72));
        }
        // v73
		pos = [r0 • v73, r1 • v73] + [t[0], t[1]];
		if(inside(pos,F) && inside (v73, T))
        {
		score += (F(pos) - T(v73)) * (F(pos) - T(v73));
        }
        // v74
		pos = [r0 • v74, r1 • v74] + [t[0], t[1]];
		if(inside(pos,F) && inside (v74, T))
        {
		score += (F(pos) - T(v74)) * (F(pos) - T(v74));
        }
        // v75
		pos = [r0 • v75, r1 • v75] + [t[0], t[1]];
		if(inside(pos,F) && inside (v75, T))
        {
		score += (F(pos) - T(v75)) * (F(pos) - T(v75));
        }
        // v76
		pos = [r0 • v76, r1 • v76] + [t[0], t[1]];
		if(inside(pos,F) && inside (v76, T))
        {
		score += (F(pos) - T(v76)) * (F(pos) - T(v76));
        }
        // v77
		pos = [r0 • v77, r1 • v77] + [t[0], t[1]];
		if(inside(pos,F) && inside (v77, T))
        {
		score += (F(pos) - T(v77)) * (F(pos) - T(v77));
        }
        // v78
		pos = [r0 • v78, r1 • v78] + [t[0], t[1]];
		if(inside(pos,F) && inside (v78, T))
        {
		score += (F(pos) - T(v78)) * (F(pos) - T(v78));
        }
        // v79
		pos = [r0 • v79, r1 • v79] + [t[0], t[1]];
		if(inside(pos,F) && inside (v79, T))
        {
		score += (F(pos) - T(v79)) * (F(pos) - T(v79));
        }
        // v80
		pos = [r0 • v80, r1 • v80] + [t[0], t[1]];
		if(inside(pos,F) && inside (v80, T))
        {
		score += (F(pos) - T(v80)) * (F(pos) - T(v80));
        }
        // v81
		pos = [r0 • v81, r1 • v81] + [t[0], t[1]];
		if(inside(pos,F) && inside (v81, T))
        {
		score += (F(pos) - T(v81)) * (F(pos) - T(v81));
        }
        // v82
		pos = [r0 • v82, r1 • v82] + [t[0], t[1]];
		if(inside(pos,F) && inside (v82, T))
        {
		score += (F(pos) - T(v82)) * (F(pos) - T(v82));
        }
        // v83
		pos = [r0 • v83, r1 • v83] + [t[0], t[1]];
		if(inside(pos,F) && inside (v83, T))
        {
		score += (F(pos) - T(v83)) * (F(pos) - T(v83));
        }
        // v84
		pos = [r0 • v84, r1 • v84] + [t[0], t[1]];
		if(inside(pos,F) && inside (v84, T))
        {
		score += (F(pos) - T(v84)) * (F(pos) - T(v84));
        }
        // v85
		pos = [r0 • v85, r1 • v85] + [t[0], t[1]];
		if(inside(pos,F) && inside (v85, T))
        {
		score += (F(pos) - T(v85)) * (F(pos) - T(v85));
        }
        // v86
		pos = [r0 • v86, r1 • v86] + [t[0], t[1]];
		if(inside(pos,F) && inside (v86, T))
        {
		score += (F(pos) - T(v86)) * (F(pos) - T(v86));
        }
        // v87
		pos = [r0 • v87, r1 • v87] + [t[0], t[1]];
		if(inside(pos,F) && inside (v87, T))
        {
		score += (F(pos) - T(v87)) * (F(pos) - T(v87));
        }
        // v88
		pos = [r0 • v88, r1 • v88] + [t[0], t[1]];
		if(inside(pos,F) && inside (v88, T))
        {
		score += (F(pos) - T(v88)) * (F(pos) - T(v88));
        }
        // v89
		pos = [r0 • v89, r1 • v89] + [t[0], t[1]];
		if(inside(pos,F) && inside (v89, T))
        {
		score += (F(pos) - T(v89)) * (F(pos) - T(v89));
        }
        // v90
		pos = [r0 • v90, r1 • v90] + [t[0], t[1]];
		if(inside(pos,F) && inside (v90, T))
        {
		score += (F(pos) - T(v90)) * (F(pos) - T(v90));
        }
        // v91
		pos = [r0 • v91, r1 • v91] + [t[0], t[1]];
		if(inside(pos,F) && inside (v91, T))
        {
		score += (F(pos) - T(v91)) * (F(pos) - T(v91));
        }
        // v92
		pos = [r0 • v92, r1 • v92] + [t[0], t[1]];
		if(inside(pos,F) && inside (v92, T))
        {
		score += (F(pos) - T(v92)) * (F(pos) - T(v92));
        }
        // v93
		pos = [r0 • v93, r1 • v93] + [t[0], t[1]];
		if(inside(pos,F) && inside (v93, T))
        {
		score += (F(pos) - T(v93)) * (F(pos) - T(v93));
        }
        // v94
		pos = [r0 • v94, r1 • v94] + [t[0], t[1]];
		if(inside(pos,F) && inside (v94, T))
        {
		score += (F(pos) - T(v94)) * (F(pos) - T(v94));
        }
        // v95
		pos = [r0 • v95, r1 • v95] + [t[0], t[1]];
		if(inside(pos,F) && inside (v95, T))
        {
		score += (F(pos) - T(v95)) * (F(pos) - T(v95));
        }
        // v96
		pos = [r0 • v96, r1 • v96] + [t[0], t[1]];
		if(inside(pos,F) && inside (v96, T))
        {
		score += (F(pos) - T(v96)) * (F(pos) - T(v96));
        }
        // v97
		pos = [r0 • v97, r1 • v97] + [t[0], t[1]];
		if(inside(pos,F) && inside (v97, T))
        {
		score += (F(pos) - T(v97)) * (F(pos) - T(v97));
        }
        // v98
		pos = [r0 • v98, r1 • v98] + [t[0], t[1]];
		if(inside(pos,F) && inside (v98, T))
        {
		score += (F(pos) - T(v98)) * (F(pos) - T(v98));
        }
        // v99
		pos = [r0 • v99, r1 • v99] + [t[0], t[1]];
		if(inside(pos,F) && inside (v99, T))
        {
		score += (F(pos) - T(v99)) * (F(pos) - T(v99));
        }
        // ********************************************
		stabilize;
	}
}

initially [ R(i, j, k) | i in 0..60, j in 0..60 , k in 0..60];

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0