Merge pull request #92 from WIAS-PDELib/chmerdon/type_fixes

fixed type of result_kernels in some operators (needed for automatic differentiation)

Author: pjaap

.gitignore

@@ -0,0 +1,8 @@
+**.png
+**.html
+**.js
+**.css
+**.json
+**.inv
+docs/build/
+Manifest.toml

CHANGELOG.md

@@ -1,5 +1,10 @@
 # CHANGES
 
+## v1.7.1
+
+### Fixed
+  - fixed type issues in LinearOperator and ItemIntegrator
+
 ## v1.7.0
 
 ### Added

Project.toml

@@ -1,6 +1,6 @@
 name = "ExtendableFEM"
 uuid = "a722555e-65e0-4074-a036-ca7ce79a4aed"
-version = "1.7.0"
+version = "1.7.1"
 authors = ["Christian Merdon <merdon@wias-berlin.de>", "Patrick Jaap <patrick.jaap@wias-berlin.de>"]
 
 [deps]

src/common_operators/item_integrator.jl

@@ -191,7 +191,7 @@ function build_assembler!(O::ItemIntegrator{Tv}, FE_args::Array{<:FEVectorBlock,
         function assembly_loop(b, sol::Array{<:FEVectorBlock, 1}, items, EG::ElementGeometries, QF::QuadratureRule, BE_args::Array{<:FEEvaluator, 1}, L2G::L2GTransformer, QPinfos::QPInfos)
 
             ## prepare parameters
-            result_kernel = zeros(Tv, resultdim)
+            result_kernel = zeros(eltype(b), resultdim)
             input_args = zeros(Tv, op_offsets_args[end])
             ndofs_args::Array{Int, 1} = [size(BE.cvals, 2) for BE in BE_args]
             weights, xref = QF.w, QF.xref

src/common_operators/linear_operator.jl

@@ -372,7 +372,7 @@ function build_assembler!(b, O::LinearOperator{Tv}, FE_test, FE_args; time = 0.0
 
             ## prepare parameters
             input_args = zeros(Tv, op_offsets_args[end])
-            result_kernel = zeros(Tv, op_offsets_test[end])
+            result_kernel = zeros(T, op_offsets_test[end])
             offsets_test = [FE_test[j].offset for j in 1:length(FES_test)]
 
             ndofs_test::Array{Int, 1} = [get_ndofs(AT, FE, EG) for FE in FETypes_test]
@@ -584,7 +584,7 @@ function build_assembler!(b, O::LinearOperator{Tv}, FE_test::Array{<:FEVectorBlo
         function assembly_loop(b::AbstractVector{T}, items, EG::ElementGeometries, QF::QuadratureRule, BE_test::Array{<:FEEvaluator, 1}, BE_test_vals::Array{Array{Tv, 3}, 1}, L2G::L2GTransformer, QPinfos::QPInfos) where {T}
 
             ## prepare parameters
-            result_kernel = zeros(Tv, op_offsets_test[end])
+            result_kernel = zeros(T, op_offsets_test[end])
             offsets_test = [FE_test[j].offset for j in 1:length(FES_test)]
 
             ndofs_test::Array{Int, 1} = [get_ndofs(AT, FE, EG) for FE in FETypes_test]