Skip to content
Snippets Groups Projects
Commit ff925638 authored by Stephan Hilb's avatar Stephan Hilb
Browse files

support function names

parent 252e7332
Branches
Tags
No related merge requests found
...@@ -7,5 +7,6 @@ version = "0.1.0" ...@@ -7,5 +7,6 @@ version = "0.1.0"
Ferrite = "c061ca5d-56c9-439f-9c0e-210fe06d3992" Ferrite = "c061ca5d-56c9-439f-9c0e-210fe06d3992"
Fewaku = "eaeb555d-1c71-5d84-9937-0fdf4ae0c3db" Fewaku = "eaeb555d-1c71-5d84-9937-0fdf4ae0c3db"
SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf" SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
UnicodePlots = "b8865327-cd53-5732-bb35-84acbb429228" UnicodePlots = "b8865327-cd53-5732-bb35-84acbb429228"
WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192" WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
...@@ -14,11 +14,12 @@ end ...@@ -14,11 +14,12 @@ end
struct FeFunction struct FeFunction
mapper::Mapper mapper::Mapper
data::Vector{Float64} data::Vector{Float64}
name::String
end end
function FeFunction(mapper) function FeFunction(mapper, name=string(gensym("f")))
data = Vector{Float64}(undef, mapper.ndofs) data = Vector{Float64}(undef, mapper.ndofs)
return FeFunction(mapper, data) return FeFunction(mapper, data, name)
end end
# P1 Elements (1st order Lagrange) # P1 Elements (1st order Lagrange)
...@@ -76,11 +77,11 @@ append_data!(vtkfile, f::FeFunction) = append_data!(vtkfile, f, f.mapper) ...@@ -76,11 +77,11 @@ append_data!(vtkfile, f::FeFunction) = append_data!(vtkfile, f, f.mapper)
function append_data!(vtkfile, f::FeFunction, ::Mapper{P1}) function append_data!(vtkfile, f::FeFunction, ::Mapper{P1})
# FIXME: vector-valued data # FIXME: vector-valued data
fdata = reshape(f.data, 1, :) fdata = reshape(f.data, 1, :)
vtk_point_data(vtkfile, fdata, "f") vtk_point_data(vtkfile, fdata, f.name)
end end
function append_data!(vtkfile, f::FeFunction, ::Mapper{DP0}) function append_data!(vtkfile, f::FeFunction, ::Mapper{DP0})
# FIXME: vector-valued data # FIXME: vector-valued data
fdata = reshape(f.data, 1, :) fdata = reshape(f.data, 1, :)
vtk_cell_data(vtkfile, fdata, "f") vtk_cell_data(vtkfile, fdata, f.name)
end end
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment