Merge pull request #130 from MultiSimOLab/sqrt-dev

Stable `sqrt`

Author: miguelmaso

Project.toml

@@ -34,7 +34,7 @@ AbstractTrees = "0.4"
 BlockArrays = "1.9"
 DrWatson = "2.19"
 ForwardDiff = "1.3"
-Gridap = "0.19.7"
+Gridap = "0.19.8"
 GridapDistributed = "0.4"
 GridapGmsh = "0.7"
 GridapPETSc = "0.5"

src/TensorAlgebra/Functions.jl

@@ -1,15 +1,4 @@
 
-"""
-    sqrt(A::TensorValue{3})::TensorValue{3}
-
-Compute the square root of a 3x3 matrix by means of eigen decomposition.
-"""
-function sqrt(A::TensorValue{3})
-  λ, Q = eigen(A)
-  λ = sqrt.(λ)
-  TensorValue{3}(λ[1]*Q[1:3]*Q[1:3]' + λ[2]*Q[4:6]*Q[4:6]' + λ[3]*Q[7:9]*Q[7:9]')
-end
-
 
 """
     cof(A::TensorValue)::TensorValue

src/TensorAlgebra/TensorAlgebra.jl

@@ -6,7 +6,6 @@ using StaticArrays
 using LinearAlgebra
 import Base: *
 import Base: +
-import Base: sqrt
 
 export (*)
 export (×ᵢ⁴)
@@ -32,7 +31,6 @@ export δᵢₗδⱼₖ2D
 export δᵢⱼδₖₗ3D
 export δᵢₖδⱼₗ3D
 export δᵢₗδⱼₖ3D
-export sqrt
 export cof
 export contraction_IP_JPKL
 export contraction_IP_PJKL
@@ -42,9 +40,9 @@ export Ellipsoid
 
 # outer ⊗ \otimes
 # inner ⊙ \odot
-# cross × \times
-# sum +
-# dot ⋅ * 
+# cross ×  \times
+# sum   +  +
+# dot   ⋅  * 
 
 include("FlatIndexing.jl")