- Hi all i have a simple question. i really appreciate your help. i want to assemble linear vector of weak form and my testfespace in 2d but it gives me error. do you know how can i solve this?
order = 1
reffe_Disp = ReferenceFE(lagrangian ,VectorValue{2,Float64},order)
V0_Disp = TestFESpace(model,reffe_Disp;conformity =:H1,
dirichlet_tags = ["LeftSupport","RightSupport","LoadLine"],
dirichlet_masks =[(false ,true), (true ,true), (false ,true)])
b_vec = assemble_vector(v->(∫(-v)fem_params.dΓ_Load), fem_params.V0_Disp)
this is error
MethodError: no method matching Float64(::VectorValue{2, Float64})
Closest candidates are:
(::Type{T})(::T) where T<:Number at boot.jl:772
(::Type{T})(::AbstractChar) where T<:Union{AbstractChar, Number} at char.jl:50
(::Type{T})(::Base.TwicePrecision) where T<:Number at twiceprecision.jl:266
...
Stacktrace:
[1] convert(#unused#::Type{Float64}, x::VectorValue{2, Float64})
@ Base .\number.jl:7
[2] setindex!(A::Vector{Float64}, x::VectorValue{2, Float64}, i1::Int64)
@ Base .\array.jl:966
[3] setindex!
@ .\multidimensional.jl:674 [inlined]
[4] add_entry!
@ C:\Users\marya\.julia\packages\Gridap\971dU\src\Algebra\AlgebraInterfaces.jl:98 [inlined]
[5] add_entry!
@ C:\Users\marya\.julia\packages\Gridap\971dU\src\Algebra\AlgebraInterfaces.jl:87 [inlined]
[6] _add_entries!
@ C:\Users\marya\.julia\packages\Gridap\971dU\src\Algebra\AlgebraInterfaces.jl:174 [inlined]
[7] add_entries!
@ C:\Users\marya\.julia\packages\Gridap\971dU\src\Algebra\AlgebraInterfaces.jl:158 [inlined]
[8] evaluate!
@ C:\Users\marya\.julia\packages\Gridap\971dU\src\Arrays\AlgebraMaps.jl:51 [inlined]
[9] _numeric_loop_vector!(vec::Vector{Float64}, caches::Tuple{Nothing, Tuple{Tuple{Nothing, Gridap.Arrays.CachedVector{VectorValue{2, Float64}, Vector{VectorValue{2, Float64}}}, Tuple{Nothing, Nothing, Tuple{Tuple{Tuple{Tuple{Nothing, Nothing, Tuple{Tuple{Tuple{Nothing, Nothing, Tuple{Tuple{Tuple{Nothing, Gridap.Arrays.CachedVector{VectorValue{2, Float64}, Vector{VectorValue{2, Float64}}}, Tuple{Tuple{Tuple{Nothing, Gridap.Arrays.CachedVector{Int32, Vector{Int32}}, Tuple{Nothing}}, Gridap.Arrays.IndexItemPair{Int64, Vector{Int32}}}}}, Gridap.Arrays.IndexItemPair{Int64, Vector{VectorValue{2, Float64}}}}, Nothing}}, Gridap.Arrays.IndexItemPair{Int64, TensorValue{1, 2, Float64, 2}}}}}, Gridap.Arrays.IndexItemPair{Int64, ConstantField{TensorValue{1, 2, Float64, 2}}}}, Gridap.Arrays.CachedVector{TensorValue{1, 2, Float64, 2}, Vector{TensorValue{1, 2, Float64, 2}}}, Tuple{Nothing}}, Gridap.Arrays.IndexItemPair{Int64, Vector{TensorValue{1, 2, Float64, 2}}}}}}, Gridap.Arrays.IndexItemPair{Int64, Vector{VectorValue{2, Float64}}}}, Tuple{Tuple{Nothing, Gridap.Arrays.CachedVector{Int32, Vector{Int32}}, Tuple{Nothing}}, Gridap.Arrays.IndexItemPair{Int64, Vector{Int32}}}}, cell_vals::Gridap.Arrays.LazyArray{FillArrays.Fill{IntegrationMap, 1, Tuple{Base.OneTo{Int64}}}, Vector{VectorValue{2, Float64}}, 1, Tuple{Gridap.Geometry.FaceCompressedVector{Matrix{VectorValue{2, Float64}}, Gridap.Geometry.FaceToCellGlue{Vector{Int64}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Gridap.Arrays.LazyArray{FillArrays.Fill{Reindex{Vector{Int8}}, 1, Tuple{Base.OneTo{Int64}}}, Int8, 1, Tuple{Vector{Int64}}}, Vector{Int8}}}, Gridap.Arrays.CompressedArray{Vector{Float64}, 1, Vector{Vector{Float64}}, Gridap.Arrays.LazyArray{FillArrays.Fill{Reindex{Vector{Int8}}, 1, Tuple{Base.OneTo{Int64}}}, Int8, 1, Tuple{Vector{Int64}}}}, Gridap.Arrays.LazyArray{Gridap.Arrays.LazyArray{FillArrays.Fill{typeof(constant_field), 1, Tuple{Base.OneTo{Int64}}}, ConstantField{TensorValue{1, 2, Float64, 2}}, 1, Tuple{Gridap.Arrays.LazyArray{FillArrays.Fill{Gridap.Fields.LinearCombinationMap{Colon}, 1, Tuple{Base.OneTo{Int64}}}, TensorValue{1, 2, Float64, 2}, 1, Tuple{Gridap.Arrays.LazyArray{FillArrays.Fill{Broadcasting{Reindex{Vector{VectorValue{2, Float64}}}}, 1, Tuple{Base.OneTo{Int64}}}, Vector{VectorValue{2, Float64}}, 1, Tuple{Gridap.Arrays.LazyArray{FillArrays.Fill{Reindex{Gridap.Arrays.Table{Int32, Vector{Int32}, Vector{Int32}}}, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int32}, 1, Tuple{Vector{Int64}}}}}, Gridap.Arrays.CompressedArray{Vector{VectorValue{1, Float64}}, 1, Vector{Vector{VectorValue{1, Float64}}}, Gridap.Arrays.LazyArray{FillArrays.Fill{Reindex{Vector{Int8}}, 1, Tuple{Base.OneTo{Int64}}}, Int8, 1, Tuple{Vector{Int64}}}}}}}}, Vector{TensorValue{1, 2, Float64, 2}}, 1, Tuple{Gridap.Arrays.CompressedArray{Vector{VectorValue{1, Float64}}, 1, Vector{Vector{VectorValue{1, Float64}}}, Gridap.Arrays.LazyArray{FillArrays.Fill{Reindex{Vector{Int8}}, 1, Tuple{Base.OneTo{Int64}}}, Int8, 1, Tuple{Vector{Int64}}}}}}}}, cell_rows::Gridap.Arrays.LazyArray{FillArrays.Fill{Reindex{Gridap.Arrays.Table{Int32, Vector{Int32}, Vector{Int32}}}, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int32}, 1, Tuple{Vector{Int32}}})
@ Gridap.FESpaces C:\Users\marya\.julia\packages\Gridap\971dU\src\FESpaces\SparseMatrixAssemblers.jl:319
[10] numeric_loop_vector!(b::Vector{Float64}, a::Gridap.FESpaces.GenericSparseMatrixAssembler, vecdata::Tuple{Vector{Any}, Vector{Any}})
@ Gridap.FESpaces C:\Users\marya\.julia\packages\Gridap\971dU\src\FESpaces\SparseMatrixAssemblers.jl:306
[11] assemble_vector(a::Gridap.FESpaces.GenericSparseMatrixAssembler, vecdata::Tuple{Vector{Any}, Vector{Any}})
@ Gridap.FESpaces C:\Users\marya\.julia\packages\Gridap\971dU\src\FESpaces\SparseMatrixAssemblers.jl:79
[12] assemble_vector(f::var"#90#91", a::Gridap.FESpaces.GenericSparseMatrixAssembler, V::Gridap.FESpaces.UnconstrainedFESpace{Vector{Float64}, Gridap.FESpaces.NodeToDofGlue{VectorValue{2, Int32}}})
@ Gridap.FESpaces C:\Users\marya\.julia\packages\Gridap\971dU\src\FESpaces\Assemblers.jl:296
[13] assemble_vector(f::Function, V::Gridap.FESpaces.UnconstrainedFESpace{Vector{Float64}, Gridap.FESpaces.NodeToDofGlue{VectorValue{2, Int32}}})
@ Gridap.FESpaces C:\Users\marya\.julia\packages\Gridap\971dU\src\FESpaces\Assemblers.jl:352
[14] top-level scope
@ In[42]:11