module Diagrams.ThreeD.Shapes
(
Ellipsoid(..)
, sphere
) where
import Graphics.Rendering.Diagrams
import Diagrams.ThreeD.Types
import Diagrams.Util
import Data.Monoid
data Ellipsoid = Ellipsoid T3
type instance V Ellipsoid = R3
instance Transformable Ellipsoid where
transform t1 (Ellipsoid t2) = Ellipsoid (t1 <> t2)
sphere :: (Backend b R3, Renderable Ellipsoid b) => Diagram b R3
sphere = mkAD (Prim $ Ellipsoid mempty)
(Bounds sphereBounds)
mempty
(Query sphereQuery)
where sphereBounds (x,y,z) = 1 / sqrt(x*x + y*y + z*z)
sphereQuery (P (x,y,z)) = Any $ x*x + y*y + z*z <= 1