ftp.nice.ch/pub/next/graphics/3d/geomview.1.4.1.s.tar.gz#/Geomview/src/lib/geometry/hpoint3/hplane3.c

This is hplane3.c in view mode; [Download] [Up]

/* Copyright (c) 1992 The Geometry Center; University of Minnesota
   1300 South Second Street;  Minneapolis, MN  55454, USA;
   
This file is part of geomview/OOGL. geomview/OOGL is free software;
you can redistribute it and/or modify it only under the terms given in
the file COPYING, which you should have received along with this file.
This and other related software may be obtained via anonymous ftp from
geom.umn.edu; email: software@geom.umn.edu. */

/* Authors: Charlie Gunn, Pat Hanrahan, Stuart Levy, Tamara Munzner, Mark Phillips */

#
/*
**	hplane3.c - procedural interface to 3D plane geometry
**
**	pat hanrahan
*/

#include <math.h>
#include "hg4.h"
#include "hpoint3.h"
#include "hplane3.h"
#include "hline3.h"
#include "transform3.h"
#include "tolerance.h"

HPlane3 HPl3Ideal = { 0., 0., 0., 1. };

HPlane3 *
HPl3Create()
{
    return (HPlane3 *) Hg4Create();
}

void
HPl3Delete( pl )
    HPlane3 *pl;
{
    Hg4Delete( (Hg4Tensor1 *)pl );
}

void
HPl3Print( pl )
    HPlane3 *pl;
{
    Hg4Print( (Hg4Tensor1 *)pl );
}

void
HPl3Copy( pl1, pl2 )
    HPlane3 *pl1, *pl2;
{
    Hg4Copy( (Hg4Tensor1 *)pl1, (Hg4Tensor1 *)pl2 );
}

void
HPl3From( pl, a, b, c, d )
    HPlane3 *pl;
    HPl3Coord a, b, c, d;
{
    Hg4From( (Hg4Tensor1 *)pl, a, b, c, d );
}

int
HPl3From3HPt3s( pl, pt1, pt2, pt3 )
    HPlane3 *pl;
    HPoint3 *pt1, *pt2, *pt3;
{
    return Hg4Intersect3( 
	(Hg4Tensor1 *)pt1, (Hg4Tensor1 *)pt2, (Hg4Tensor1 *)pt3,
	(Hg4Tensor1 *)pl, 1 );
}

int
HPl3From2HLn3s( pl, ln1, ln2 )
    HPlane3 *pl;
    HLine3 *ln1, *ln2;
{
    HPoint3 pt;

    return HLn3IntersectHLn3( ln1, ln2, pl, &pt );
}

int
HPl3IntersectHPl3( pl1, pl2, ln )
    HPlane3 *pl1, *pl2;
    HLine3 *ln;
{
    return HLn3From2HPl3s( ln, pl1, pl2 );
}

void
HPl3Pencil( t1, pl1, t2, pl2, pl3 )
    HPl3Coord t1, t2;
    HPlane3 *pl1, *pl2, *pl3;
{
    Hg4Pencil( t1, (Hg4Tensor1 *)pl1, t2, (Hg4Tensor2 *)pl2, 
	(Hg4Tensor1 *)pl3 );
}

float
HPl3DotHPt3( pl, pt )
    HPlane3 *pl;
    HPoint3 *pt;
{
    return Hg4ContractPiQi( (Hg4Tensor1 *)pl, (Hg4Tensor1 *)pt );
}

int
HPl3Undefined( pl )
    HPlane3 *pl;
{
    return Hg4Undefined( (Hg4Tensor1 *)pl );
}

int
HPl3Infinity( pl )
    HPlane3 *pl;
{
    return Hg4Infinity( (Hg4Tensor1 *) pl, 1 );
}

int
HPl3Compare( pl1, pl2 )
    HPlane3 *pl1, *pl2;
{
    return Hg4Compare( (Hg4Tensor1 *)pl1, (Hg4Tensor1 *)pl2 );
}

int
HPl3CoincidentHPt3( pl, pt )
    HPlane3 *pl;
    HPoint3 *pt;
{
    return fzero(HPl3DotHPt3(pl,pt));
}

int
HPl3CoincidentHLn3( pl, ln )
    HPlane3 *pl;
    HLine3 *ln;
{
    HPoint3 pt;

    return HLn3IntersectHPt3( ln, pl, &pt );
}

int
HPl3CoincidentHPl3( pl1, pl2 )
    HPlane3 *pl1, *pl2;
{
    return Hg4Coincident( (Hg4Tensor1 *)pl1, (Hg4Tensor1 *)pl2 );
}

void
HPl3Transform( T, pl1, pl2 )
    Transform3 T;
    HPlane3 *pl1, *pl2;
{
    Hg4Transform( T, (Hg4Tensor1 *)pl1, (Hg4Tensor1 *)pl2 );
}

void
HPl3TransformN( T, pl1, pl2, n )
    Transform3 T;
    HPlane3 *pl1, *pl2;
    int n;
{
    while( n-- )
	Hg4Transform( T, (Hg4Tensor1 *)pl1++, (Hg4Tensor1 *)pl2++ );
}

void
HPl3Dual( pl, pt )
    HPlane3 *pl;
    HPoint3 *pt;
{
    pt->x = pl->a;
    pt->y = pl->b;
    pt->z = pl->c;
    pt->w = pl->d;
}

void
HPl3Perp( pl, pt )
    HPlane3 *pl;
    HPoint3 *pt;
{
/*
    HLine3 ln, lndual;
    HQuadric Q;

    HLn3From2HPl3s( &ln, pl, &HPl3Ideal );
    HLn3CorrelateHLn3( &ln, Q, &lndual );
    HLn3Perp( &lndual, pt ):
*/
    /* if HPl3Ideal == (0,0,0,1) */
    pt->x = pl->a;
    pt->y = pl->b;
    pt->z = pl->c;
    pt->w = 0;
}

HPl3Coord
HPl3Angle( pl1, pl2 )
    HPlane3 *pl1, *pl2;
{
    return 0.0;
}

These are the contents of the former NiCE NeXT User Group NeXTSTEP/OpenStep software archive, currently hosted by Netfuture.ch.