/*-
* #%L
* Fiji distribution of ImageJ for the life sciences.
* %%
* Copyright (C) 2007 - 2017 Fiji developers.
* %%
* This program is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as
* published by the Free Software Foundation, either version 2 of the
* License, or (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public
* License along with this program. If not, see
* <http://www.gnu.org/licenses/gpl-2.0.html>.
* #L%
*/
package spim.vecmath;
/*
* Copyright 1997-2008 Sun Microsystems, Inc. All Rights Reserved.
* DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.
*
* This code is free software; you can redistribute it and/or modify it
* under the terms of the GNU General Public License version 2 only, as
* published by the Free Software Foundation. Sun designates this
* particular file as subject to the "Classpath" exception as provided
* by Sun in the LICENSE file that accompanied this code.
*
* This code is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
* version 2 for more details (a copy is included in the LICENSE file that
* accompanied this code).
*
* You should have received a copy of the GNU General Public License version
* 2 along with this work; if not, write to the Free Software Foundation,
* Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
*
* Please contact Sun Microsystems, Inc., 4150 Network Circle, Santa Clara,
* CA 95054 USA or visit www.sun.com if you need additional information or
* have any questions.
*
*/
/**
* A 4 element vector represented by double precision floating point x,y,z,w
* coordinates.
*
*/
public class Point4d extends Tuple4d implements java.io.Serializable
{
// Compatible with 1.1
static final long serialVersionUID = 1733471895962736949L;
/**
* Constructs and initializes a Point4d from the specified xyzw coordinates.
*
* @param x
* the x coordinate
* @param y
* the y coordinate
* @param z
* the z coordinate
* @param w
* the w coordinate
*/
public Point4d( double x, double y, double z, double w )
{
super( x, y, z, w );
}
/**
* Constructs and initializes a Point4d from the coordinates contained in
* the array.
*
* @param p
* the array of length 4 containing xyzw in order
*/
public Point4d( double[] p )
{
super( p );
}
/**
* Constructs and initializes a Point4d from the specified Point4d.
*
* @param p1
* the Point4d containing the initialization x y z w data
*/
public Point4d( Point4d p1 )
{
super( p1 );
}
/**
* Constructs and initializes a Point4d from the specified Point4f.
*
* @param p1
* the Point4f containing the initialization x y z w data
*/
public Point4d( Point4f p1 )
{
super( p1 );
}
/**
* Constructs and initializes a Point4d from the specified Tuple4f.
*
* @param t1
* the Tuple4f containing the initialization x y z w data
*/
public Point4d( Tuple4f t1 )
{
super( t1 );
}
/**
* Constructs and initializes a Point4d from the specified Tuple4d.
*
* @param t1
* the Tuple4d containing the initialization x y z w data
*/
public Point4d( Tuple4d t1 )
{
super( t1 );
}
/**
* Constructs and initializes a Point4d from the specified Tuple3d. The
* x,y,z components of this point are set to the corresponding components of
* tuple t1. The w component of this point is set to 1.
*
* @param t1
* the tuple to be copied
*
* @since vecmath 1.2
*/
public Point4d( Tuple3d t1 )
{
super( t1.x, t1.y, t1.z, 1.0 );
}
/**
* Constructs and initializes a Point4d to (0,0,0,0).
*/
public Point4d()
{
super();
}
/**
* Sets the x,y,z components of this point to the corresponding components
* of tuple t1. The w component of this point is set to 1.
*
* @param t1
* the tuple to be copied
*
* @since vecmath 1.2
*/
public final void set( Tuple3d t1 )
{
this.x = t1.x;
this.y = t1.y;
this.z = t1.z;
this.w = 1.0;
}
/**
* Returns the square of the distance between this point and point p1.
*
* @param p1
* the first point
* @return the square of distance between this point and point p1
*/
public final double distanceSquared( Point4d p1 )
{
double dx, dy, dz, dw;
dx = this.x - p1.x;
dy = this.y - p1.y;
dz = this.z - p1.z;
dw = this.w - p1.w;
return ( dx * dx + dy * dy + dz * dz + dw * dw );
}
/**
* Returns the distance between this point and point p1.
*
* @param p1
* the first point
* @return the distance between these this point and point p1.
*/
public final double distance( Point4d p1 )
{
double dx, dy, dz, dw;
dx = this.x - p1.x;
dy = this.y - p1.y;
dz = this.z - p1.z;
dw = this.w - p1.w;
return Math.sqrt( dx * dx + dy * dy + dz * dz + dw * dw );
}
/**
* Computes the L-1 (Manhattan) distance between this point and point p1.
* The L-1 distance is equal to: abs(x1-x2) + abs(y1-y2) + abs(z1-z2) +
* abs(w1-w2).
*
* @param p1
* the other point
* @return the L-1 distance
*/
public final double distanceL1( Point4d p1 )
{
return Math.abs( this.x - p1.x ) + Math.abs( this.y - p1.y )
+ Math.abs( this.z - p1.z ) + Math.abs( this.w - p1.w );
}
/**
* Computes the L-infinite distance between this point and point p1. The
* L-infinite distance is equal to MAX[abs(x1-x2), abs(y1-y2), abs(z1-z2),
* abs(w1-w2)].
*
* @param p1
* the other point
* @return the L-infinite distance
*/
public final double distanceLinf( Point4d p1 )
{
double t1, t2;
t1 = Math.max( Math.abs( this.x - p1.x ), Math.abs( this.y - p1.y ) );
t2 = Math.max( Math.abs( this.z - p1.z ), Math.abs( this.w - p1.w ) );
return Math.max( t1, t2 );
}
/**
* Multiplies each of the x,y,z components of the Point4d parameter by 1/w,
* places the projected values into this point, and places a 1 as the w
* parameter of this point.
*
* @param p1
* the source Point4d, which is not modified
*/
public final void project( Point4d p1 )
{
double oneOw;
oneOw = 1 / p1.w;
x = p1.x * oneOw;
y = p1.y * oneOw;
z = p1.z * oneOw;
w = 1.0;
}
}