source: osgVisual/src/manip_Spacemouse/manip_nodeTrackerSpaceMouse.cpp @ 94

Last change on this file since 94 was 88, checked in by Torben Dannhauer, 14 years ago

Moved memory leak detection from source file to headerfile. Its still in the class but at least not in the source file.

The leak detection works, but the false positives are not stopped.
Use Linux/Valgrind? to make your final leak detection beyond the easy first approach in MSVC

File size: 20.4 KB
Line 
1/* -*-c++-*- osgVisual - Copyright (C) 2009-2010 Torben Dannhauer
2 *
3 * This library is based on OpenSceneGraph, open source and may be redistributed and/or modified under
4 * the terms of the OpenSceneGraph Public License (OSGPL) version 0.0 or
5 * (at your option) any later version.  The full license is in LICENSE file
6 * included with this distribution, and on the openscenegraph.org website.
7 *
8 * osgVisual requires for some proprietary modules a license from the correspondig manufacturer.
9 * You have to aquire licenses for all used proprietary modules.
10 *
11 * This library is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14 * OpenSceneGraph Public License for more details.
15*/
16
17#include <manip_nodeTrackerSpaceMouse.h>
18
19
20
21using namespace osg;
22using namespace osgVisual;
23
24
25NodeTrackerSpaceMouse::NodeTrackerSpaceMouse(SpaceMouse* spacemouse) : _spaceMouse(spacemouse)
26{
27    _trackerMode = NODE_CENTER_AND_ROTATION; 
28    _rotationMode = TRACKBALL; 
29
30    _distance = 1.0;
31
32    _thrown = false;
33        _autohoming = false;
34        _ah_init = false;
35        _distanceDependendAV = false;
36
37        TZ=0;
38        RX=0;
39        RY=0;
40        RZ=0;
41
42}
43
44
45NodeTrackerSpaceMouse::~NodeTrackerSpaceMouse()
46{
47}
48
49void NodeTrackerSpaceMouse::setTrackNode(osg::Node* node)
50{
51    if (!node)
52    {
53        osg::notify(osg::NOTICE)<<"NodeTrackerSpaceMouse::setTrackNode(Node*):  Unable to set tracked node due to null Node*"<<std::endl;
54        return;
55    }
56
57    osg::NodePathList parentNodePaths = node->getParentalNodePaths();
58
59    if (!parentNodePaths.empty())
60    {
61        osg::notify(osg::INFO)<<"NodeTrackerSpaceMouse::setTrackNode(Node*): Path set"<<std::endl;
62        setTrackNodePath(parentNodePaths[0]);
63    }
64    else
65    {
66        osg::notify(osg::NOTICE)<<"NodeTrackerSpaceMouse::setTrackNode(Node*): Unable to set tracked node due to empty parental path."<<std::endl;
67    }
68}
69
70osg::Node* NodeTrackerSpaceMouse::getTrackNode()
71{
72    osg::NodePath nodePath;
73    if (_trackNodePath.getNodePath(nodePath)) return nodePath.back();
74    else return 0;
75}
76
77const osg::Node* NodeTrackerSpaceMouse::getTrackNode() const
78{
79    osg::NodePath nodePath;
80    if (_trackNodePath.getNodePath(nodePath)) return nodePath.back();
81    else return 0;
82}
83
84void NodeTrackerSpaceMouse::setTrackerMode(TrackerMode mode)
85{
86    _trackerMode = mode;
87}
88
89void NodeTrackerSpaceMouse::setRotationMode(RotationMode mode)
90{
91    _rotationMode = mode;
92    if (getAutoComputeHomePosition()) computeHomePosition();   
93}
94
95void NodeTrackerSpaceMouse::setNode(osg::Node* node)
96{
97    _node = node;
98   
99    if (_node.get())
100    {
101        const osg::BoundingSphere& boundingSphere=_node->getBound();
102        const float minimumDistanceScale = 0.001f;
103        _minimumDistance = osg::clampBetween(
104            float(boundingSphere._radius) * minimumDistanceScale,
105            0.00001f,1.0f);
106           
107        osg::notify(osg::INFO)<<"Setting Tracker manipulator _minimumDistance to "<<_minimumDistance<<std::endl;
108    }
109    if (getAutoComputeHomePosition()) computeHomePosition();   
110}
111
112const osg::Node* NodeTrackerSpaceMouse::getNode() const
113{
114    return _node.get();
115}
116
117
118osg::Node* NodeTrackerSpaceMouse::getNode()
119{
120    return _node.get();
121}
122
123
124void NodeTrackerSpaceMouse::home(const GUIEventAdapter& ,GUIActionAdapter& us)
125{
126    if (getAutoComputeHomePosition()) computeHomePosition();
127
128    computePosition(_homeEye, _homeCenter, _homeUp);
129    us.requestRedraw();
130}
131
132void NodeTrackerSpaceMouse::computeHomePosition()
133{
134    osg::NodePath nodePath;
135    _trackNodePath.getNodePath(nodePath);
136
137    osg::Node* node = nodePath.empty() ? getNode() : nodePath.back();
138   
139    if(node)
140    {
141        const osg::BoundingSphere& boundingSphere=node->getBound();
142
143        setHomePosition(boundingSphere._center+osg::Vec3( 0.0,-3.5f * boundingSphere._radius,0.0f),
144                        boundingSphere._center,
145                        osg::Vec3(0.0f,0.0f,1.0f),
146                        _autoComputeHomePosition);
147    }
148}
149
150
151
152void NodeTrackerSpaceMouse::init(const GUIEventAdapter& ,GUIActionAdapter& )
153{
154    flushMouseEventStack();
155}
156
157
158void NodeTrackerSpaceMouse::getUsage(osg::ApplicationUsage& usage) const
159{
160    usage.addKeyboardMouseBinding("SN Tracker: Space","Toggle camera auto homing");
161        usage.addKeyboardMouseBinding("SN Tracker: m","Toggle Nodetracker Mode");
162        usage.addKeyboardMouseBinding("SN Tracker: n","Toggle distance dependend angular velocity");
163    usage.addKeyboardMouseBinding("SN Tracker: +","When in stereo, increase the fusion distance");
164    usage.addKeyboardMouseBinding("SN Tracker: -","When in stereo, reduce the fusion distance");
165}
166
167bool NodeTrackerSpaceMouse::handle(const GUIEventAdapter& ea,GUIActionAdapter& us)
168{
169    switch(ea.getEventType())
170    {
171        case(GUIEventAdapter::PUSH):
172        {
173            flushMouseEventStack();
174            addMouseEvent(ea);
175            if (calcMovement()) us.requestRedraw();
176            us.requestContinuousUpdate(false);
177            _thrown = false;
178            return true;
179        }
180
181        case(GUIEventAdapter::RELEASE):
182        {
183            if (ea.getButtonMask()==0)
184            {
185
186                double timeSinceLastRecordEvent = _ga_t0.valid() ? (ea.getTime() - _ga_t0->getTime()) : DBL_MAX;
187                if (timeSinceLastRecordEvent>0.02) flushMouseEventStack();
188
189                if (isMouseMoving())
190                {
191                    if (calcMovement())
192                    {
193                        us.requestRedraw();
194                        us.requestContinuousUpdate(true);
195                        _thrown = true;
196                    }
197                }
198                else
199                {
200                    flushMouseEventStack();
201                    addMouseEvent(ea);
202                    if (calcMovement()) us.requestRedraw();
203                    us.requestContinuousUpdate(false);
204                    _thrown = false;
205                }
206
207            }
208            else
209            {
210                flushMouseEventStack();
211                addMouseEvent(ea);
212                if (calcMovement()) us.requestRedraw();
213                us.requestContinuousUpdate(false);
214                _thrown = false;
215            }
216            return true;
217        }
218
219        case(GUIEventAdapter::DRAG):
220        {
221            addMouseEvent(ea);
222            if (calcMovement()) us.requestRedraw();
223            us.requestContinuousUpdate(false);
224            _thrown = false;
225            return true;
226        }
227
228        case(GUIEventAdapter::MOVE):
229        {
230            return false;
231        }
232
233        case(GUIEventAdapter::KEYDOWN):
234            if (ea.getKey()==' ')
235            {
236                                if (_autohoming)
237                                        _autohoming=false;
238                                else
239                                        _autohoming=true;
240                return true;
241            }
242                        if (ea.getKey()=='m')
243            {
244                                switch(_trackerMode)
245                                {
246                                        case NODE_CENTER: _trackerMode=NODE_CENTER_AND_AZIM;
247                                                break;
248                                        case NODE_CENTER_AND_AZIM: _trackerMode=NODE_CENTER_AND_ROTATION;
249                                                break;
250                                        case  NODE_CENTER_AND_ROTATION: _trackerMode=NODE_CENTER;
251                                                break;
252                                        default: _trackerMode = NODE_CENTER;
253                                                break;
254                                };
255                return true;
256            }
257                        if (ea.getKey()=='n')
258            {
259                                if (_distanceDependendAV)
260                                        _distanceDependendAV=false;
261                                else
262                                        _distanceDependendAV=true;
263                return true;
264            }
265            return false;
266        case(GUIEventAdapter::FRAME):
267            if (_thrown)
268            {
269                if (calcMovement()) us.requestRedraw();
270            }
271
272                        computeHomePosition();
273                        computePosition(_homeEye, _homeCenter, _homeUp);
274                        if (calcMovementSpaceMouse())
275                                us.requestRedraw();
276
277            return false;
278        default:
279            return false;
280    }
281}
282
283
284bool NodeTrackerSpaceMouse::isMouseMoving()
285{
286    if (_ga_t0.get()==NULL || _ga_t1.get()==NULL) return false;
287
288    static const float velocity = 0.1f;
289
290    float dx = _ga_t0->getXnormalized()-_ga_t1->getXnormalized();
291    float dy = _ga_t0->getYnormalized()-_ga_t1->getYnormalized();
292    float len = sqrtf(dx*dx+dy*dy);
293    float dt = _ga_t0->getTime()-_ga_t1->getTime();
294
295    return (len>dt*velocity);
296}
297
298
299void NodeTrackerSpaceMouse::flushMouseEventStack()
300{
301    _ga_t1 = NULL;
302    _ga_t0 = NULL;
303}
304
305
306void NodeTrackerSpaceMouse::addMouseEvent(const GUIEventAdapter& ea)
307{
308    _ga_t1 = _ga_t0;
309    _ga_t0 = &ea;
310}
311
312void NodeTrackerSpaceMouse::setByMatrix(const osg::Matrixd& matrix)
313{
314    osg::Vec3d eye,center,up;
315    matrix.getLookAt(eye,center,up,_distance);
316    computePosition(eye,center,up);
317}
318
319void NodeTrackerSpaceMouse::computeNodeCenterAndRotation(osg::Vec3d& nodeCenter, osg::Quat& nodeRotation) const
320{
321    osg::Matrixd localToWorld, worldToLocal;
322    osg::NodePath nodePath;
323    if (_trackNodePath.getNodePath(nodePath))
324    {
325        worldToLocal = osg::computeWorldToLocal(nodePath);
326        localToWorld = osg::computeLocalToWorld(nodePath);
327        nodeCenter = osg::Vec3d(nodePath.back()->getBound().center())*localToWorld;
328    }
329    else
330    {
331        nodeCenter = osg::Vec3d(0.0f,0.0f,0.0f)*localToWorld;
332    }
333
334
335    switch(_trackerMode)
336    {
337        case(NODE_CENTER_AND_AZIM):
338        {
339            CoordinateFrame coordinateFrame = getCoordinateFrame(nodeCenter);
340            osg::Matrixd localToFrame(localToWorld*osg::Matrixd::inverse(coordinateFrame));
341
342            double azim = atan2(-localToFrame(0,1),localToFrame(0,0));
343            osg::Quat nodeRotationRelToFrame, rotationOfFrame;
344            nodeRotationRelToFrame.makeRotate(-azim,0.0,0.0,1.0);
345            rotationOfFrame = coordinateFrame.getRotate();
346            nodeRotation = nodeRotationRelToFrame*rotationOfFrame;
347            break;
348        }
349        case(NODE_CENTER_AND_ROTATION):
350        {
351            // scale the matrix to get rid of any scales before we extract the rotation.
352            double sx = 1.0/sqrt(localToWorld(0,0)*localToWorld(0,0) + localToWorld(1,0)*localToWorld(1,0) + localToWorld(2,0)*localToWorld(2,0));
353            double sy = 1.0/sqrt(localToWorld(0,1)*localToWorld(0,1) + localToWorld(1,1)*localToWorld(1,1) + localToWorld(2,1)*localToWorld(2,1));
354            double sz = 1.0/sqrt(localToWorld(0,2)*localToWorld(0,2) + localToWorld(1,2)*localToWorld(1,2) + localToWorld(2,2)*localToWorld(2,2));
355            localToWorld = localToWorld*osg::Matrixd::scale(sx,sy,sz);
356
357            nodeRotation = localToWorld.getRotate();
358            break;
359        }
360        case(NODE_CENTER):
361        default:
362        {
363            CoordinateFrame coordinateFrame = getCoordinateFrame(nodeCenter);
364            nodeRotation = coordinateFrame.getRotate();
365            break;
366        }
367    }
368
369}
370
371
372osg::Matrixd NodeTrackerSpaceMouse::getMatrix() const
373{
374    osg::Vec3d nodeCenter;
375    osg::Quat nodeRotation;
376    computeNodeCenterAndRotation(nodeCenter,nodeRotation);
377    return osg::Matrixd::translate(0.0,0.0,_distance)*osg::Matrixd::rotate(_rotation)*osg::Matrixd::rotate(nodeRotation)*osg::Matrix::translate(nodeCenter);
378}
379
380osg::Matrixd NodeTrackerSpaceMouse::getInverseMatrix() const
381{
382    osg::Vec3d nodeCenter;
383    osg::Quat nodeRotation;
384    computeNodeCenterAndRotation(nodeCenter,nodeRotation);
385    return osg::Matrixd::translate(-nodeCenter)*osg::Matrixd::rotate(nodeRotation.inverse())*osg::Matrixd::rotate(_rotation.inverse())*osg::Matrixd::translate(0.0,0.0,-_distance);
386}
387
388void NodeTrackerSpaceMouse::computePosition(const osg::Vec3d& eye,const osg::Vec3d& center,const osg::Vec3d& up)
389{
390    if (!_node) return;
391
392    // compute rotation matrix
393    osg::Vec3 lv(center-eye);
394    _distance = lv.length();
395   
396    osg::Matrixd lookat;
397    lookat.makeLookAt(eye,center,up);
398   
399    _rotation = lookat.getRotate().inverse();
400}
401
402bool NodeTrackerSpaceMouse::calcMovement()
403{
404    // return if less then two events have been added.
405    if (_ga_t0.get()==NULL || _ga_t1.get()==NULL) return false;
406
407    double dx = _ga_t0->getXnormalized()-_ga_t1->getXnormalized();
408    double dy = _ga_t0->getYnormalized()-_ga_t1->getYnormalized();
409
410
411    float distance = sqrtf(dx*dx + dy*dy);
412    // return if movement is too fast, indicating an error in event values or change in screen.
413    if (distance>0.5)
414    {
415        return false;
416    }
417   
418    // return if there is no movement.
419    if (distance==0.0f)
420    {
421        return false;
422    }
423
424    osg::Vec3d nodeCenter;
425    osg::Quat nodeRotation;
426    computeNodeCenterAndRotation(nodeCenter, nodeRotation);
427
428    unsigned int buttonMask = _ga_t1->getButtonMask();
429
430    if (buttonMask==GUIEventAdapter::LEFT_MOUSE_BUTTON)
431    {
432
433        if (_rotationMode==TRACKBALL)
434        {
435
436            // rotate camera.
437            osg::Vec3 axis;
438            double angle;
439
440            double px0 = _ga_t0->getXnormalized();
441            double py0 = _ga_t0->getYnormalized();
442
443            double px1 = _ga_t1->getXnormalized();
444            double py1 = _ga_t1->getYnormalized();
445
446
447            trackball(axis,angle,px1,py1,px0,py0);
448
449            osg::Quat new_rotate;
450            new_rotate.makeRotate(angle,axis);
451
452            _rotation = _rotation*new_rotate;
453        }
454        else
455        {
456            osg::Matrix rotation_matrix;
457            rotation_matrix.makeRotate(_rotation);
458
459            osg::Vec3d lookVector = -getUpVector(rotation_matrix);
460            osg::Vec3d sideVector = getSideVector(rotation_matrix);
461            osg::Vec3d upVector = getFrontVector(rotation_matrix);
462           
463            osg::Vec3d localUp(0.0f,0.0f,1.0f);
464
465            osg::Vec3d forwardVector = localUp^sideVector;
466            sideVector = forwardVector^localUp;
467
468            forwardVector.normalize();
469            sideVector.normalize();
470           
471            osg::Quat rotate_elevation;
472            rotate_elevation.makeRotate(dy,sideVector);
473
474            osg::Quat rotate_azim;
475            rotate_azim.makeRotate(-dx,localUp);
476           
477            _rotation = _rotation * rotate_elevation * rotate_azim;
478           
479        }
480       
481        return true;
482
483    }
484    else if (buttonMask==GUIEventAdapter::MIDDLE_MOUSE_BUTTON ||
485        buttonMask==(GUIEventAdapter::LEFT_MOUSE_BUTTON|GUIEventAdapter::RIGHT_MOUSE_BUTTON))
486    {
487        return true;
488    }
489    else if (buttonMask==GUIEventAdapter::RIGHT_MOUSE_BUTTON)
490    {
491
492        // zoom model.
493
494        double fd = _distance;
495        double scale = 1.0f+dy;
496        if (fd*scale>_minimumDistance)
497        {
498
499            _distance *= scale;
500
501        } else
502        {
503            _distance = _minimumDistance;
504        }
505       
506        return true;
507
508    }
509
510    return false;
511}
512
513bool NodeTrackerSpaceMouse::calcMovementSpaceMouse()
514{
515        // TRACKBALL mode is not possible, because it uses the mouseevent, which does not happen with Spacemouse.
516
517   double dTX, dTY, dTZ;        // Only TZ is used later.
518    _spaceMouse->getTranslations(dTX, dTY, dTZ);
519
520    double dRX, dRY, dRZ;
521    _spaceMouse->getRotations(dRX, dRY, dRZ);
522
523        //OSG_NOTIFY( osg::INFO ) << "Recieved Spacemousevalues: dRX:" << dRX << ", dRY:" << dRY << ", dRZ:" << dRZ << ", dTX:" << dTX << std::endl;
524
525        if(!_ah_init)
526        {
527                TZ=dTZ;
528                RX=dRX;
529                RY=dRY;
530                RZ=dRZ;
531                _ah_init=true;
532        }
533
534        if (_autohoming)        // The factors are required to allow macroscopic movements.
535        {
536                TZ=dTZ*5;
537                if (!_distanceDependendAV)
538                {
539                        RX=dRX*100;
540                        RY=dRY*100;
541                        RZ=dRZ*100;
542                }
543        }
544        else            // NOT Autohoming
545        {
546                TZ+=dTZ;
547                //OSG_NOTIFY( osg::INFO ) << "Stored Spacemousevalues: RX:" << RX << ", RY:" << RY << ", RZ:" << RZ << ", TX:" << TX << std::endl;
548                if (_distanceDependendAV)
549                {
550                        RX+=(dRX*800/TZ);
551                        RY+=(dRY*800/TZ);
552                        RZ+=(dRZ*800/TZ);
553                }
554                else
555                {
556                        RX+=dRX;
557                        RY+=dRY;
558                        RZ+=dRZ;
559                }
560        }
561
562    osg::Vec3d nodeCenter;
563    osg::Quat nodeRotation;
564    computeNodeCenterAndRotation(nodeCenter, nodeRotation);
565
566    // ROTATION PART
567    if (_rotationMode==TRACKBALL)
568    {
569                // return if less then two events have been added.
570                if (_ga_t0.get()==NULL || _ga_t1.get()==NULL) return false;
571
572        // rotate camera.
573        osg::Vec3 axis;
574        double angle;
575
576        double px0 = _ga_t0->getXnormalized();
577        double py0 = _ga_t0->getYnormalized();
578
579        double px1 = _ga_t1->getXnormalized();
580        double py1 = _ga_t1->getYnormalized();
581
582
583        trackball(axis,angle,px1,py1,px0,py0);
584
585        osg::Quat new_rotate;
586        new_rotate.makeRotate(angle,axis);
587
588        _rotation = _rotation*new_rotate;
589    }
590    else
591    {
592        osg::Matrix rotation_matrix;
593        rotation_matrix.makeRotate(_rotation);
594
595        osg::Vec3d lookVector = -getUpVector(rotation_matrix);
596        osg::Vec3d sideVector = getSideVector(rotation_matrix);
597        osg::Vec3d upVector = getFrontVector(rotation_matrix);
598       
599        osg::Vec3d localUp(0.0f,0.0f,1.0f);
600
601        osg::Vec3d forwardVector = localUp^sideVector;
602        sideVector = forwardVector^localUp;
603
604        forwardVector.normalize();
605        sideVector.normalize();
606       
607        osg::Quat rotate_elevation;
608                rotate_elevation.makeRotate(osg::DegreesToRadians(RX*90), sideVector);
609
610        osg::Quat rotate_azim;
611                rotate_azim.makeRotate(osg::DegreesToRadians(RY*90), localUp);
612       
613        _rotation = _rotation * rotate_elevation * rotate_azim;
614       
615    }
616
617    // TRANSLATION PART
618    double scale = 1.0f+TZ/100;
619        if (_distance*scale>_minimumDistance)
620    {
621
622        _distance *= scale;
623
624    } else
625    {
626                TZ = (_minimumDistance/_distance - 1)*100;      // Reset TZ to the value f(_minimalDistance, _distance, scale-function).
627        _distance = _minimumDistance;
628    }
629    return true;
630}
631
632void NodeTrackerSpaceMouse::clampOrientation()
633{
634}
635
636
637/*
638 * This size should really be based on the distance from the center of
639 * rotation to the point on the object underneath the mouse.  That
640 * point would then track the mouse as closely as possible.  This is a
641 * simple example, though, so that is left as an Exercise for the
642 * Programmer.
643 */
644const float TRACKBALLSIZE = 0.8f;
645
646/*
647 * Ok, simulate a track-ball.  Project the points onto the virtual
648 * trackball, then figure out the axis of rotation, which is the cross
649 * product of P1 P2 and O P1 (O is the center of the ball, 0,0,0)
650 * Note:  This is a deformed trackball-- is a trackball in the center,
651 * but is deformed into a hyperbolic sheet of rotation away from the
652 * center.  This particular function was chosen after trying out
653 * several variations.
654 *
655 * It is assumed that the arguments to this routine are in the range
656 * (-1.0 ... 1.0)
657 */
658void NodeTrackerSpaceMouse::trackball(osg::Vec3& axis,double & angle, double  p1x, double  p1y, double  p2x, double  p2y)
659{
660    /*
661     * First, figure out z-coordinates for projection of P1 and P2 to
662     * deformed sphere
663     */
664
665    osg::Matrix rotation_matrix(_rotation);
666
667
668    osg::Vec3d uv = osg::Vec3d(0.0,1.0,0.0)*rotation_matrix;
669    osg::Vec3d sv = osg::Vec3d(1.0,0.0,0.0)*rotation_matrix;
670    osg::Vec3d lv = osg::Vec3d(0.0,0.0,-1.0)*rotation_matrix;
671
672    osg::Vec3d p1 = sv*p1x+uv*p1y-lv*tb_project_to_sphere(TRACKBALLSIZE,p1x,p1y);
673    osg::Vec3d p2 = sv*p2x+uv*p2y-lv*tb_project_to_sphere(TRACKBALLSIZE,p2x,p2y);
674
675    /*
676     *  Now, we want the cross product of P1 and P2
677     */
678
679// Robert,
680//
681// This was the quick 'n' dirty  fix to get the trackball doing the right
682// thing after fixing the Quat rotations to be right-handed.  You may want
683// to do something more elegant.
684//   axis = p1^p2;
685axis = p2^p1;
686    axis.normalize();
687
688    /*
689     *  Figure out how much to rotate around that axis.
690     */
691    double t = (p2-p1).length() / (2.0*TRACKBALLSIZE);
692
693    /*
694     * Avoid problems with out-of-control values...
695     */
696    if (t > 1.0) t = 1.0;
697    if (t < -1.0) t = -1.0;
698    angle = inRadians(asin(t));
699
700}
701
702
703/*
704 * Project an x,y pair onto a sphere of radius r OR a hyperbolic sheet
705 * if we are away from the center of the sphere.
706 */
707double NodeTrackerSpaceMouse::tb_project_to_sphere(double  r, double  x, double  y)
708{
709    float d, t, z;
710
711    d = sqrt(x*x + y*y);
712                                 /* Inside sphere */
713    if (d < r * 0.70710678118654752440)
714    {
715        z = sqrt(r*r - d*d);
716    }                            /* On hyperbola */
717    else
718    {
719        t = r / 1.41421356237309504880;
720        z = t*t / d;
721    }
722    return z;
723}
Note: See TracBrowser for help on using the repository browser.