F܈      !"#$%&'()*+, - . / 0 1 2 3 4 5 6 7 8 9 : ; < = > ? @ A B C D E FGHIJKLMN O P Q R S T U V W X Y Z [ \ ] ^ _ ` a b c d e f g h i j k l m n o p q r s t u v w x y z { | } ~portable provisionalEdward Kmett <ekmett@gmail.com>None portable provisionalEdward Kmett <ekmett@gmail.com>None9A vector is an additive group with additional structure. The zero vector Compute the sum of two vectors V2 1 2 ^+^ V2 3 4V2 4 6+Compute the difference between two vectors V2 4 5 - V2 3 1V2 1 4*Linearly interpolate between two vectors. Apply a function to merge the  'non-zero'= components of two vectors, unioning the rest of the values. * For a dense vector this is equivalent to . + For a sparse vector this is equivalent to . 3Apply a function to the components of two vectors. * For a dense vector this is equivalent to . + For a sparse vector this is equivalent to . !Compute the negation of a vector negated (V2 2 4) V2 (-2) (-4)  Compute the left scalar product  2 *^ V2 3 4V2 6 8 !Compute the right scalar product  V2 3 4 ^* 2V2 6 8 +Compute division by a scalar on the right. BProduce a default basis for a vector space. If the dimensionality 2 of the vector space is not statically known, see  . :Produce a default basis for a vector space from which the  argument is drawn. )Produce a diagonal matrix from a vector. &Outer (tensor) product of two vectors #    portable provisionalEdward Kmett <ekmett@gmail.com> Safe-InferredEProvides a fairly subjective test to see if a quantity is near zero. nearZero (1e-11 :: Double)FalsenearZero (1e-17 :: Double)TruenearZero (1e-5 :: Float)FalsenearZero (1e-7 :: Float)True&Determine if a quantity is near zero.  a  1e-12  a  1e-6  a  1e-12  a  1e-6 non-portable experimentalEdward Kmett <ekmett@gmail.com>NoneFree and sparse inner product/metric spaces. ;Compute the inner product of two vectors or (equivalently)  convert a vector f a into a covector f a -> a. V2 1 2 `dot` V2 3 4119Compute the squared norm. The name quadrance arises from  Norman J. Wildberger's rational trigonometry. (Compute the quadrance of the difference ;Compute the distance between two vectors in a metric space /Compute the norm of a vector in a metric space *Convert a non-zero vector to unit vector.  Normalize a  functor to have unit . This function $ does not change the functor if its  is 0 or 1.  non-portable experimentalEdward Kmett <ekmett@gmail.com> Safe-InferredA  f+ is corepresentable if it is isomorphic to (x -> a) E for some x. Nearly all such functors can be represented by choosing x to be ? the set of lenses that are polymorphic in the contents of the ,  that is to say x = Rep f is a valid choice of x for (nearly) every   Representable .  non-portable experimentalEdward Kmett <ekmett@gmail.com>NoneA 0-dimensional vector pure 1 :: V0 IntV0V0 + V0V0 non-portable experimentalEdward Kmett <ekmett@gmail.com>None6A space that distinguishes 2 orthogonal basis vectors  and  , but may have more.  V2 1 2 ^._x1V2 1 2 & _x .~ 3V2 3 2    :: Lens' (t a) a   V2 1 2 ^._y2V2 1 2 & _y .~ 3V2 1 3     :: Lens' (t a) a !   ! :: Lens' (t a) (" a) "A 2-dimensional vector pure 1 :: V2 IntV2 1 1V2 1 2 + V2 3 4V2 4 6V2 1 2 * V2 3 4V2 3 8 sum (V2 1 2)3$+the counter-clockwise perpendicular vector perp $ V2 10 20 V2 (-20) 10 !"#$ !"#$"# !$ !"#$ non-portable experimentalEdward Kmett <ekmett@gmail.com>None%7A space that distinguishes 3 orthogonal basis vectors: ,  , and &. (It may have more) &   & :: Lens' (t a) a '   ' :: Lens' (t a) (( a) (A 3-dimensional vector *cross product +scalar triple product %&'()*+  !%&'()*+ ()*+ !%&'%&'()*+  non-portable experimentalEdward Kmett <ekmett@gmail.com>None,4A space that distinguishes orthogonal basis vectors ,  , &, -. (It may have more.) -   - :: Lens' (t a) a .   . :: Lens' (t a) (/ a) /A 4-dimensional vector. 1OConvert a 3-dimensional affine vector into a 4-dimensional homogeneous vector. 2NConvert a 3-dimensional affine point into a 4-dimensional homogeneous vector. ,-./012 !%&',-./012/012 !%&',-.,-./012  non-portable experimentalEdward Kmett <ekmett@gmail.com>None38Plcker coordinates for lines in a 3-dimensional space. 5[Given a pair of points represented by homogeneous coordinates generate Plcker coordinates  for the line through them. 6OThese elements form a basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   6 :: Lens' (3 a) a  7 :: Lens' (3 a) a  8 :: Lens' (3 a) a  9 :: Lens' (3 a) a  : :: Lens' (3 a) a  ; :: Lens' (3 a) a 7OThese elements form a basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   6 :: Lens' (3 a) a  7 :: Lens' (3 a) a  8 :: Lens' (3 a) a  9 :: Lens' (3 a) a  : :: Lens' (3 a) a  ; :: Lens' (3 a) a 8OThese elements form a basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   6 :: Lens' (3 a) a  7 :: Lens' (3 a) a  8 :: Lens' (3 a) a  9 :: Lens' (3 a) a  : :: Lens' (3 a) a  ; :: Lens' (3 a) a 9OThese elements form a basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   6 :: Lens' (3 a) a  7 :: Lens' (3 a) a  8 :: Lens' (3 a) a  9 :: Lens' (3 a) a  : :: Lens' (3 a) a  ; :: Lens' (3 a) a :OThese elements form a basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   6 :: Lens' (3 a) a  7 :: Lens' (3 a) a  8 :: Lens' (3 a) a  9 :: Lens' (3 a) a  : :: Lens' (3 a) a  ; :: Lens' (3 a) a ;OThese elements form a basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   6 :: Lens' (3 a) a  7 :: Lens' (3 a) a  8 :: Lens' (3 a) a  9 :: Lens' (3 a) a  : :: Lens' (3 a) a  ; :: Lens' (3 a) a <ZThese elements form an alternate basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   < ::  a => Lens' (3 a) a  = ::  a => Lens' (3 a) a  > ::  a => Lens' (3 a) a  ? ::  a => Lens' (3 a) a  @ ::  a => Lens' (3 a) a  A ::  a => Lens' (3 a) a =ZThese elements form an alternate basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   < ::  a => Lens' (3 a) a  = ::  a => Lens' (3 a) a  > ::  a => Lens' (3 a) a  ? ::  a => Lens' (3 a) a  @ ::  a => Lens' (3 a) a  A ::  a => Lens' (3 a) a >ZThese elements form an alternate basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   < ::  a => Lens' (3 a) a  = ::  a => Lens' (3 a) a  > ::  a => Lens' (3 a) a  ? ::  a => Lens' (3 a) a  @ ::  a => Lens' (3 a) a  A ::  a => Lens' (3 a) a ?ZThese elements form an alternate basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   < ::  a => Lens' (3 a) a  = ::  a => Lens' (3 a) a  > ::  a => Lens' (3 a) a  ? ::  a => Lens' (3 a) a  @ ::  a => Lens' (3 a) a  A ::  a => Lens' (3 a) a @ZThese elements form an alternate basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   < ::  a => Lens' (3 a) a  = ::  a => Lens' (3 a) a  > ::  a => Lens' (3 a) a  ? ::  a => Lens' (3 a) a  @ ::  a => Lens' (3 a) a  A ::  a => Lens' (3 a) a AZThese elements form an alternate basis for the Plcker space, or the Grassmanian manifold Gr(2,V4).   < ::  a => Lens' (3 a) a  = ::  a => Lens' (3 a) a  > ::  a => Lens' (3 a) a  ? ::  a => Lens' (3 a) a  @ ::  a => Lens' (3 a) a  A ::  a => Lens' (3 a) a BValid Plcker coordinates p will have B p  0 PThat said, floating point makes a mockery of this claim, so you may want to use . CThis isn'Yt th actual metric because this bilinear form gives rise to an isotropic quadratic space DrChecks if the line is near-isotropic (isotropic vectors in this quadratic space represent lines in real 3d space) E:Checks if the two vectors intersect (or nearly intersect) &3456789:;<=>?@ABCDE     3456789:;<=>?@ABCDE34BDC5E678<;@=A9>:?%3456789:;<=>?@ABCDE     NoneFGHIJKLM !"#$%&' FGHIJKLM FGJHIKLMFGHIJKLM !"#$%&'  non-portable experimentalEdward Kmett <ekmett@gmail.com> Safe-InferredNAn involutive ring O<Conjugate a value. This defaults to the trivial involution. conjugate (1 :+ 2) 1.0 :+ (-2.0) conjugate 11NO()*+,-./01234567NONONO()*+,-./01234567  non-portable experimentalEdward Kmett <ekmett@gmail.com>NoneP0A vector space that includes the basis elements U, V, Q and R Q   Q :: Lens' (t a) a R   R :: Lens' (t a) a S   S :: Lens' (t a) (V3 a) T0A vector space that includes the basis elements U and V U   U :: Lens' (t a) a V   V :: Lens' (t a) a W Quaternions 8%quadrance of the imaginary component Y norm of the imaginary component Zraise a W to a scalar power 91Helper for calculating with specific branch cuts :1Helper for calculating with specific branch cuts [; with a specified branch cut. \< with a specified branch cut. ]= with a specified branch cut. ^> with a specified branch cut. _? with a specified branch cut. `@ with a specified branch cut. a8Spherical linear interpolation between two quaternions. bApply a rotation to a vector. cc axis theta builds a W representing a  rotation of theta radians about axis. .PQRSTUVWXAB8YZ9:[\]^_`abcCDEFGHIJKLMNOPQRSTUVWPQRSTUVWXYZ[\]^_`abcWXTUVPQRSa[\]^_`YZbc(PQRSTUVWXAB8YZ9:[\]^_`abcCDEFGHIJKLMNOPQRSTUVW  non-portable experimentalEdward Kmett <ekmett@gmail.com>Noned+A 4x3 matrix with row-major representation e+A 4x4 matrix with row-major representation f+A 3x3 matrix with row-major representation g+A 2x2 matrix with row-major representation hdMatrix product. This can compute mixed dense-dense, sparse-dense and sparse-sparse matrix products. :V2 (V3 1 2 3) (V3 4 5 6) !*! V3 (V2 1 2) (V2 3 4) (V2 4 5)V2 (V2 19 25) (V2 43 58)SV2 (fromList [(1,2)]) (fromList [(2,3)]) !*! fromList [(1,V3 0 0 1), (2, V3 0 0 5)]V2 (V3 0 0 2) (V3 0 0 15)iEntry-wise matrix addition. 5V2 (V3 1 2 3) (V3 4 5 6) !+! V2 (V3 7 8 9) (V3 1 2 3)V2 (V3 8 10 12) (V3 5 7 9)jEntry-wise matrix subtraction. 5V2 (V3 1 2 3) (V3 4 5 6) !-! V2 (V3 7 8 9) (V3 1 2 3)!V2 (V3 (-6) (-6) (-6)) (V3 3 3 3)kMatrix * column vector $V2 (V3 1 2 3) (V3 4 5 6) !* V3 7 8 9 V2 50 122lRow vector * matrix "V2 1 2 *! V2 (V3 3 4 5) (V3 6 7 8) V3 15 18 21mScalar-matrix product 5 *!! V2 (V2 1 2) (V2 3 4)V2 (V2 5 10) (V2 15 20)nMatrix-scalar product V2 (V2 1 2) (V2 3 4) !!* 5V2 (V2 5 10) (V2 15 20)o+Hermitian conjugate or conjugate transpose :adjoint (V2 (V2 (1 :+ 2) (3 :+ 4)) (V2 (5 :+ 6) (7 :+ 8)))LV2 (V2 (1.0 :+ (-2.0)) (5.0 :+ (-6.0))) (V2 (3.0 :+ (-4.0)) (7.0 :+ (-8.0)))pCompute the trace of a matrix trace (V2 (V2 a b) (V2 c d))a + dq$Build a rotation matrix from a unit W. r;Build a transformation matrix from a rotation matrix and a  translation vector. s=Build a transformation matrix from a rotation expressed as a  W and a translation vector. tAConvert from a 4x3 matrix to a 4x4 matrix, extending it with the  [ 0 0 0 1 ] column vector u8Convert a 3x3 matrix to a 4x4 matrix extending it with 0's in the new row and column. v2x2 identity matrix. eye2V2 (V2 1 0) (V2 0 1)w3x3 identity matrix. eye3#V3 (V3 1 0 0) (V3 0 1 0) (V3 0 0 1)x4x4 identity matrix. eye46V4 (V4 1 0 0 0) (V4 0 1 0 0) (V4 0 0 1 0) (V4 0 0 0 1)y@Extract the translation vector (first three entries of the last " column) from a 3x4 or 4x4 matrix z2x2 matrix determinant. det22 (V2 (V2 a b) (V2 c d)) a * d - b * c{3x3 matrix determinant. +det33 (V3 (V3 a b c) (V3 d e f) (V3 g h i))?a * (e * i - f * h) - d * (b * i - c * h) + g * (b * f - c * e)|2x2 matrix inverse. inv22 $ V2 (V2 1 2) (V2 3 4))Just (V2 (V2 (-2.0) 1.0) (V2 1.5 (-0.5)))}3x3 matrix inverse. +inv33 $ V3 (V3 1 2 4) (V3 4 2 2) (V3 1 1 1)JJust (V3 (V3 0.0 0.5 (-1.0)) (V3 (-0.5) (-0.75) 3.5) (V3 0.5 0.25 (-1.5)))defghijklmnopqrstuvwxyz{|}defghijklmnopqrstuvwxyz{|}hijklnmogfedutz{|}vwxpyqsrdefghijklmnopqrstuvwxyz{|}portable provisionalEdward Kmett <ekmett@gmail.com>None~?A handy wrapper to help distinguish points from vectors at the  type level ;An affine space is roughly a vector space in which we have = forgotten or at least pretend to have forgotten the origin.  a .+^ (b .-. a) = b@ & (a .+^ u) .+^ v = a .+^ (u ^+^ v)@ & (a .-. b) ^+^ v = (a .+^ v) .-. q@ :Get the difference between two points as a vector offset.  Add a vector offset to a point. 'Subtract a vector offset from a point. ECompute the quadrance of the difference (the square of the distance) /Distance between two points in an affine space Vector spaces have origins. ~XYZ[\]^_`abcdefghi ~ ~~XYZ[\]^_`abcdefghi non-portable experimentalEdward Kmett <ekmett@gmail.com>Noneu  !"#$%&'()*+,-./0123456789:;<=>?@ABCDENOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{|}j !"#$%&'()*+,-.//012344567899:; < = > ? ? @ A B B C D E F G H I J K L M N O P Q R STUVWXYZ[ \ ] ^ _ ` a b c d e e f g h i j k l m n o p q r s t u v w x y z { | } ~                                  ! " # $ % & '(T)*+,-./0123456789:; < = > ? @ A B C D E F G H I J K L M NOPOQOROSOTOU V W X Y Z [ \ ] ^ _ ` a b c d e f g h i j k lmnopqrstuvwxyz{|}~ linear-1.1.2Linear.V Linear.VectorLinear.Epsilon Linear.Metric Linear.Core Linear.V0 Linear.V2 Linear.V3 Linear.V4Linear.PluckerLinear.ConjugateLinear.Quaternion Linear.Matrix Linear.AffineLinear.InstancesLinearreflection-1.3.2Data.ReflectionintAdditivezero^+^^-^lerpliftU2liftI2negated*^^*^/basisbasisFor kroneckerouterEpsilonnearZeroMetricdot quadranceqddistancenormsignorm normalizeCorecoreV0R2_x_y_xyV2perpR3_z_xyzV3crosstripleR4_w_xyzwV4vectorpointPluckerpluckerp01p02p03p23p31p12p10p20p30p32p13p21 squaredError>< isotropic intersectsVtoVectorDim reflectDimdimreifyDim reifyVector fromVector Conjugate conjugate Hamiltonian_j_k_ijk Complicated_e_i Quaternionabsipowasinqacosqatanqasinhqacoshqatanhqslerprotate axisAngleM43M44M33M22!*!!+!!-!!**!*!!!!*adjointtracefromQuaternionmkTransformationMatmkTransformation m43_to_m44 m33_to_m44eye2eye3eye4 translationdet22det33inv22inv33PointPAffineDiff.-..+^.-^qdA distanceAorigin$fTraversable1Complex$fFoldable1Complex$fTraversableComplex$fFoldableComplex$fMonadComplex $fBindComplex$fApplicativeComplex$fApplyComplex$fFunctorComplex $fBindHashMap$fApplyHashMapbaseControl.ApplicativeliftA2containers-0.5.0.0 Data.Map.Base unionWithintersectionWith GAdditivegzerogliftU2gliftI2 setElement$fAdditiveIdentity$fAdditiveComplex$fAdditive(->)$fAdditiveHashMap $fAdditiveMap$fAdditiveIntMap $fAdditive[]$fAdditiveMaybe$fAdditiveVector$fAdditiveZipList$fGAdditivePar1 $fGAdditiveM1$fGAdditiveRec1$fGAdditive:*: $fGAdditiveU1$fEpsilonCDoubleGHC.Numabsghc-prim GHC.Classes<=$fEpsilonCFloat$fEpsilonDouble$fEpsilonFloat$fMetricVector$fMetricIdentityGHC.BaseFunctor $fStorableV0 $fEpsilonV0$fDistributiveV0$fCoreV0 $fMetricV0$fFractionalV0$fNumV0 $fMonadV0$fBindV0 $fAdditiveV0$fApplicativeV0 $fApplyV0$fTraversableV0 $fFoldableV0 $fFunctorV0$fIxV2 $fStorableV2 $fEpsilonV2$fDistributiveV2$fCoreV2$fR2V2 $fMetricV2$fFractionalV2$fNumV2 $fMonadV2$fBindV2 $fAdditiveV2$fApplicativeV2 $fApplyV2$fTraversable1V2 $fFoldable1V2$fTraversableV2 $fFoldableV2 $fFunctorV2$fIxV3 $fEpsilonV3 $fStorableV3$fCoreV3$fR3V3$fR2V3$fDistributiveV3 $fMetricV3$fFractionalV3$fNumV3 $fMonadV3$fBindV3 $fAdditiveV3$fApplicativeV3 $fApplyV3$fTraversable1V3 $fFoldable1V3$fTraversableV3 $fFoldableV3 $fFunctorV3$fIxV4 $fEpsilonV4 $fStorableV4$fCoreV4$fR4V4$fR3V4$fR2V4$fDistributiveV4 $fMetricV4$fFractionalV4$fNumV4 $fMonadV4$fBindV4 $fAdditiveV4 $fApplyV4$fApplicativeV4$fTraversable1V4 $fFoldable1V4$fTraversableV4 $fFoldableV4 $fFunctorV4Num==anti$fEpsilonPlucker$fMetricPlucker$fStorablePlucker$fFractionalPlucker $fNumPlucker $fIxPlucker$fTraversable1Plucker$fFoldable1Plucker$fTraversablePlucker$fFoldablePlucker $fCorePlucker$fDistributivePlucker$fMonadPlucker $fBindPlucker$fAdditivePlucker$fApplicativePlucker$fApplyPlucker$fFunctorPlucker ReifiedDimretagDim $fMetricV $fEpsilonV $fStorableV$fDistributiveV$fCoreV $fFractionalV$fNumV $fAdditiveV$fMonadV$fBindV$fApplicativeV$fApplyV$fTraversableV $fFoldableV $fFunctorV$fDim*V$fDim*ReifiedDim $fDimNatn$fConjugateComplex$fConjugateCDouble$fConjugateCFloat$fConjugateFloat$fConjugateDouble$fConjugateWord8$fConjugateWord16$fConjugateWord32$fConjugateWord64$fConjugateWord$fConjugateInt8$fConjugateInt16$fConjugateInt32$fConjugateInt64$fConjugateInt$fConjugateIntegerqicutcutWith GHC.FloatasinacosatanasinhacoshatanhqNaN reimagine$fEpsilonQuaternion$fFloatingQuaternion$fConjugateQuaternion$fDistributiveQuaternion$fHamiltonianQuaternion$fComplicatedQuaternion$fComplicatedComplex$fMetricQuaternion$fFractionalQuaternion$fNumQuaternion$fStorableQuaternion$fTraversableQuaternion$fFoldableQuaternion$fCoreQuaternion$fIxQuaternion$fMonadQuaternion$fBindQuaternion$fAdditiveQuaternion$fApplicativeQuaternion$fApplyQuaternion$fFunctorQuaternion $fAffinePoint $fAffineV$fAffineHashMap $fAffineMap $fAffine(->)$fAffineQuaternion$fAffinePlucker $fAffineV4 $fAffineV3 $fAffineV2 $fAffineV0$fAffineVector$fAffineIdentity$fAffineIntMap $fAffineMaybe$fAffineZipList$fAffineComplex $fAffine[]