#
# Curve - point/line/curve closest/"farest" location/intersection.
#
#					Gershon Elber, May 1993.
#

#
# Curve point distance.
#
crv1 = cbspline( 3,
		 list( ctlpt( E2,  0.5, -1.0 ),
		       ctlpt( E2, -1.5, -0.5 ),
		       ctlpt( E2,  1.5,  0.0 ),
		       ctlpt( E2, -0.5,  1.0 ) ),
		 list( KV_OPEN ) );
color( crv1, green );
attrib( crv1, "width", 0.02 );

pt_param = crvptdst( crv1, point( 0, 0, 0 ), false, -0.001 );
pt_extrem = nil();
for ( i = 1, 1, listsize( pt_param ),
	pt = ceval( crv1, nth( pt_param, i ) ):
	snoc( coerce( pt, vector_type ) , pt_extrem )
    );
interact( list( axes, crv1, pt_extrem ) );

pt_min = ceval( crv1, crvptdst( crv1, point( 0, 0, 0 ), true, 0.001 ) );
interact( list( axes, crv1, coerce( pt_min, vector_type ) ) );
pt_max = ceval( crv1, crvptdst( crv1, point( 0, 0, 0 ), false, 0.001 ) );
interact( list( axes, crv1, coerce( pt_max, vector_type ) ) );

crv2 = cbspline( 3,
		 list( ctlpt( E2, 0,  1 ),
		       ctlpt( E2, 1, -1 ),
		       ctlpt( E2, 0, -1 ),
		       ctlpt( E2, 1,  1 ) ),
		 list( KV_OPEN ) );
color( crv2, green );
attrib( crv2, "width", 0.02 );

pt_param = crvptdst( crv2, point( 0, 0, 0 ), false, -0.001 );
pt_extrem = nil();
for ( i = 1, 1, listsize( pt_param ),
	pt = ceval( crv2, nth( pt_param, i ) ):
	snoc( coerce( pt, vector_type ) , pt_extrem )
    );
interact( list( axes, crv2, pt_extrem ) );

pt_min = ceval( crv2, crvptdst( crv2, point( 0, 0, 0 ), true, 0.001 ) );
interact( list( axes, crv2, coerce( pt_min, vector_type ) ) );
pt_max = ceval( crv2, crvptdst( crv2, point( 0, 0, 0 ), false, 0.001 ) );
interact( list( axes, crv2, coerce( pt_max, vector_type ) ) );

crv3a = cbspline( 3,
		  list( ctlpt( E2, 0.0,  1.0 ),
		        ctlpt( E2, 0.1,  1.0 ),
		        ctlpt( E2, 0.1,  0.1 ),
		        ctlpt( E2, 0.5,  0.1 ),
		        ctlpt( E2, 0.5, -0.1 ),
		        ctlpt( E2, 0.1, -0.1 ),
		        ctlpt( E2, 0.1, -1.0 ),
		        ctlpt( E2, 0.0, -1.0 ) ),
		  list( KV_OPEN ) );
crv3b = crv3a * rotz( 180 );
crv3 = ( crv3a + crv3b ) * trans( vector( 0.1, 0.1, 0.1 ) );
color( crv3, green );
attrib( crv3, "width", 0.02 );

pt_param = crvptdst( crv3, point( 0, 0, 0 ), false, -0.001 );
pt_extrem = nil();
for ( i = 1, 1, listsize( pt_param ),
	pt = ceval( crv3, nth( pt_param, i ) ):
	snoc( coerce( pt, vector_type ) , pt_extrem )
    );
interact( list( axes, crv3, pt_extrem ) );

pt_min = ceval( crv3, crvptdst( crv3, point( 0, 0, 0 ), true, 0.001 ) );
interact( list( axes, crv3, coerce( pt_min, vector_type ) ) );
pt_max = ceval( crv3, crvptdst( crv3, point( 0, 0, 0 ), false, 0.001 ) );
interact( list( axes, crv3, coerce( pt_max, vector_type ) ) );

#
# Curve line distance.
#
line_pt = point( -1, 1.2, 0 );
line_vec = vector( 1, -1, 0 );
line_pt2 = line_pt + ( line_vec * 2 );
line = poly( list( line_pt, line_pt2 ), true );
crv1 = cbspline( 3,
		 list( ctlpt( E2,  0.5, -1.0 ),
		       ctlpt( E2, -1.5, -0.5 ),
		       ctlpt( E2,  1.5,  0.0 ),
		       ctlpt( E2, -0.5,  1.0 ) ),
		 list( KV_OPEN ) );
color( crv1, green );
attrib( crv1, "width", 0.02 );

pt_param = crvlndst( crv1, line_pt, line_vec, false, -0.001 );
pt_extrem = nil();
for ( i = 1, 1, listsize( pt_param ),
	pt = ceval( crv1, nth( pt_param, i ) ):
	snoc( pt, pt_extrem )
    );
interact( list( axes, crv1, line, pt_extrem ) );

pt_min = ceval( crv1, crvlndst( crv1, line_pt, line_vec, true, 0.001 ) );
interact( list( axes, crv1, line, pt_min ) );
pt_max = ceval( crv1, crvlndst( crv1, line_pt, line_vec, false, 0.001 ) );
interact( list( axes, crv1, line, pt_max ) );

line_pt = point( -1, -1.5, 0 );
line_vec = vector( 1, 2, 0 );
line_pt2 = line_pt + ( line_vec * 2 );
line = poly( list( line_pt, line_pt2 ), true );
crv2 = cbspline( 3,
		 list( ctlpt( E2, 0,  1 ),
		       ctlpt( E2, 1, -1 ),
		       ctlpt( E2, 0, -1 ),
		       ctlpt( E2, 1,  1 ) ),
		 list( KV_OPEN ) );
color( crv2, green );
attrib( crv2, "width", 0.02 );

pt_param = crvlndst( crv2, line_pt, line_vec, false, -0.001 );
pt_extrem = nil();
for ( i = 1, 1, listsize( pt_param ),
	pt = ceval( crv2, nth( pt_param, i ) ):
	snoc( pt, pt_extrem )
    );
interact( list( axes, crv2, line, pt_extrem ) );

pt_min = ceval( crv2, crvlndst( crv2, line_pt, line_vec, true, 0.001 ) );
interact( list( axes, crv2, line, pt_min ) );
pt_max = ceval( crv2, crvlndst( crv2, line_pt, line_vec, false, 0.001 ) );
interact( list( axes, crv2, line, pt_max ) );

crv3a = cbspline( 3,
		  list( ctlpt( E2, 0.0,  1.0 ),
		        ctlpt( E2, 0.1,  1.0 ),
		        ctlpt( E2, 0.1,  0.1 ),
		        ctlpt( E2, 0.5,  0.1 ),
		        ctlpt( E2, 0.5, -0.1 ),
		        ctlpt( E2, 0.1, -0.1 ),
		        ctlpt( E2, 0.1, -1.0 ),
		        ctlpt( E2, 0.0, -1.0 ) ),
		  list( KV_OPEN ) );
crv3b = crv3a * rotz( 180 );
crv3 = ( crv3a + crv3b ) * trans( vector( 0.1, 0.1, 0.1 ) );
color( crv3, green );
attrib( crv3, "width", 0.02 );

pt_param = crvlndst( crv3, line_pt, line_vec, false, -0.001 );
pt_extrem = nil();
for ( i = 1, 1, listsize( pt_param ),
	pt = ceval( crv3, nth( pt_param, i ) ):
	snoc( pt, pt_extrem )
    );
interact( list( axes, crv3, line, pt_extrem ) );

pt_min = ceval( crv3, crvlndst( crv3, line_pt, line_vec, true, 0.001 ) );
interact( list( axes, crv3, line, pt_min ) );
pt_max = ceval( crv3, crvlndst( crv3, line_pt, line_vec, false, 0.001 ) );
interact( list( axes, crv3, line, pt_max ) );

#
# Curve curve intersetion/distance square scalar field.
#
crv1 = cbspline( 3,
		 list( ctlpt( E2, 0, 0 ),
		       ctlpt( E2, 0.25, 0.5 ),
		       ctlpt( E2, 0.5, 0.7 ),
		       ctlpt( E2, 1, 1 ) ),
		  list( 0.0, 0.0, 0.0, 0.5, 1.0, 1.0, 1.0 ) );
crv2 = cbspline( 3,
		 list( ctlpt( E2, 1, 0 ),
		       ctlpt( E2, 0.7, 0.25 ),
		       ctlpt( E2, 0.3, 0.5 ),
		       ctlpt( E2, 0, 1 ) ),
		 list( 0.0, 0.0, 0.0, 0.5, 1.0, 1.0, 1.0 ) );
inter_pts = ccinter( crv1, crv2, 0.0001, false );
pt_inters = nil();
pt = nil();
for ( i = 1, 1, listsize( inter_pts ),
	pt = ceval( crv1, coord( nth( inter_pts, i ), 0 ) ):
	snoc( pt, pt_inters ):
	pt = ceval( crv2, coord( nth( inter_pts, i ), 1 ) ):
	snoc( pt, pt_inters )
    );
interact( list( axes, crv1, crv2, pt_inters ) );
dist_srf = ccinter( crv1, crv2, -1, true );
color( dist_srf, magenta );
interact( list( axes, dist_srf ) );

crv1 = cbspline( 3,
		 list( ctlpt( E2, 0, 0 ),
		       ctlpt( E2, 0.25, 2 ),
		       ctlpt( E2, 0.5, -1 ),
		       ctlpt( E2, 1, 1 ) ),
		 list( 0.0, 0.0, 0.0, 0.5, 1.0, 1.0, 1.0 ) );
crv2 = cbspline( 3,
		 list( ctlpt( E2, 1, 0 ),
		       ctlpt( E2, -1, 0.25 ),
		       ctlpt( E2, 2, 0.5 ),
		       ctlpt( E2, 0, 1 ) ),
		 list( 0.0, 0.0, 0.0, 0.5, 1.0, 1.0, 1.0 ) );

inter_pts = ccinter( crv1, crv2, 0.001, false );
pt_inters = nil();
pt = nil();
for ( i = 1, 1, listsize( inter_pts ),
	pt = ceval( crv1, coord( nth( inter_pts, i ), 0 ) ):
	snoc( pt, pt_inters ):
	pt = ceval( crv2, coord( nth( inter_pts, i ), 1 ) ):
	snoc( pt, pt_inters )
    );
interact( list( axes, crv1, crv2, pt_inters ) );
dist_srf = ccinter( crv1, crv2, -1, true );
color( dist_srf, magenta );
interact( list( axes, dist_srf ) );

crv1 = cbspline( 3,
		 list( ctlpt( E2, 0, 1 ),
		       ctlpt( E2, 1, 0 ),
		       ctlpt( E2, 0, 0 ),
		       ctlpt( E2, 1, 1 ) ),
		 list( 0.0, 0.0, 0.0, 0.5, 1.0, 1.0, 1.0 ) );
inter_pts = ccinter( crv1, crv1, 0.001, true );
pt_inters = nil();
pt = nil();
for ( i = 1, 1, listsize( inter_pts ),
	pt = ceval( crv1, coord( nth( inter_pts, i ), 0 ) ):
	snoc( pt, pt_inters )
    );
interact( list( axes, crv1, pt_inters ) );
dist_srf = ccinter( crv1, crv1, -1, true );
color( dist_srf, magenta );
interact( list( axes, dist_srf ) );

free(inter_pts);
free(pt_inters);
free(crv1);
free(crv2);