Commit 3b04239f authored by thomaskummer's avatar thomaskummer
Browse files

heart geometry with patches

parent 6f8ee5dd
......@@ -12,7 +12,7 @@ TRIBITS_ADD_EXECUTABLE_AND_TEST(
)
TRIBITS_COPY_FILES_TO_BINARY_DIR(data_heart_test_fibers_data
SOURCE_FILES dataSingleVentricle dataXVentSingle dataSingleVentConv singleVentricle.mesh xVentSingle.mesh biVentCoarse.mesh biVentMedium.mesh biVentFine.mesh singleVentConv.mesh siVent.mesh
SOURCE_FILES dataSingleVentricle dataXVentSingle dataSingleVentConv singleVentricle.mesh xVentSingle.mesh biVentCoarse.mesh biVentCoarsePatches.mesh biVentMedium.mesh biVentFine.mesh singleVentConv.mesh siVent.mesh
SOURCE_DIR ${CMAKE_CURRENT_SOURCE_DIR}
)
......
......@@ -29,67 +29,74 @@ hdf5_sheets_name = sheets
[./boundary_conditions]
list = 'EndocardiumLV Epicardium EndocardiumRV Septum Support Base SeptumEdge'
list = 'EndocardiumLV Epicardium Patch1 Patch2 EndocardiumRV Septum Base'
[./EndocardiumLV]
type = Essential
flag = 50
flag = 65
mode = Full
component = 1
function = '0.0'
[../Epicardium]
type = Essential
flag = 54
flag = 69
mode = Full
component = 1
function = '1.0'
[../EndocardiumRV]
[../Patch1]
type = Essential
flag = 51
flag = 70
mode = Full
component = 1
function = '0.75'
function = '1.0'
[../Septum]
[../Patch2]
type = Essential
flag = 52
flag = 71
mode = Full
component = 1
function = '0.75'
function = '1.0'
[../Support]
type = Natural
flag = 49
[../EndocardiumRV]
type = Essential
flag = 66
mode = Full
component = 1
function = '0.0'
[../Septum]
type = Essential
flag = 67
mode = Full
component = 1
function = '1.0'
[../Base]
type = Natural
flag = 48
flag = 64
mode = Full
component = 1
function = '0.0'
[../SeptumEdge]
type = Essential
flag = 53
type = Natural
flag = 68
mode = Full
component = 1
function = '0.75'
function = '0.0'
[../]
[../space_discretization]
mesh_dir = ./
mesh_file = biVentCoarse.mesh
elem_type = P1
mesh_file = biVentCoarsePatches.mesh
elem_type = P2
[../solver]
solver_name = AztecOO #Belos
solver = cg
solver = gmres
scaling = none
output = all # none
conv = rhs
......
......@@ -32,20 +32,14 @@ SOURCE_FILES biVentCoarse
SOURCE_DIR ${CMAKE_SOURCE_DIR}/lifev/em/examples/example_EMBiVent
)
TRIBITS_COPY_FILES_TO_BINARY_DIR(biVentFine
CREATE_SYMLINK
SOURCE_FILES biVentFine
SOURCE_DIR ${CMAKE_SOURCE_DIR}/lifev/em/examples/example_EMBiVent
)
TRIBITS_COPY_FILES_TO_BINARY_DIR(biVentCoarseP2
CREATE_SYMLINK
SOURCE_FILES biVentCoarseP2
SOURCE_DIR ${CMAKE_SOURCE_DIR}/lifev/em/examples/example_EMBiVent
)
TRIBITS_COPY_FILES_TO_BINARY_DIR(biVentMediumP2
CREATE_SYMLINK
SOURCE_FILES biVentMediumP2
SOURCE_DIR ${CMAKE_SOURCE_DIR}/lifev/em/examples/example_EMBiVent
TRIBITS_COPY_FILES_TO_BINARY_DIR(biVentCoarsePatchesP2
CREATE_SYMLINK
SOURCE_FILES biVentCoarsePatchesP2
SOURCE_DIR ${CMAKE_SOURCE_DIR}/lifev/em/examples/example_EMBiVent
)
......@@ -12,10 +12,10 @@
<Time TimeType="Single" Value="0" />
<Topology
Type="Tetrahedron"
NumberOfElements="93511"
NumberOfElements="5041"
BaseOffset="0">
<DataStructure Format="HDF"
Dimensions="93511 4"
Dimensions="5041 4"
DataType="Int"
Precision="8">
FiberDirection.h5:/Connections/Values
......@@ -23,19 +23,19 @@
</Topology>
<Geometry Type="X_Y_Z">
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
FiberDirection.h5:/PointsX.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
FiberDirection.h5:/PointsY.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
FiberDirection.h5:/PointsZ.00000/Values
......@@ -48,16 +48,16 @@
Center="Node"
Name="fibers">
<DataStructure ItemType="HyperSlab"
Dimensions="21346 3"
Dimensions="1524 3"
Type="HyperSlab">
<DataStructure Dimensions="3 2"
Format="XML">
0 0
1 1
21346 3
1524 3
</DataStructure>
<DataStructure Format="HDF"
Dimensions="147495 3"
Dimensions="9405 3"
DataType="Float"
Precision="8">
FiberDirection.h5:/fibers.00000/Values
......
......@@ -12,10 +12,10 @@
<Time TimeType="Single" Value="0" />
<Topology
Type="Tetrahedron"
NumberOfElements="93511"
NumberOfElements="5041"
BaseOffset="0">
<DataStructure Format="HDF"
Dimensions="93511 4"
Dimensions="5041 4"
DataType="Int"
Precision="8">
Potential.h5:/Connections/Values
......@@ -23,19 +23,19 @@
</Topology>
<Geometry Type="X_Y_Z">
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
Potential.h5:/PointsX.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
Potential.h5:/PointsY.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
Potential.h5:/PointsZ.00000/Values
......@@ -48,16 +48,16 @@
Center="Node"
Name="potential">
<DataStructure ItemType="HyperSlab"
Dimensions="21346 1"
Dimensions="1524 1"
Type="HyperSlab">
<DataStructure Dimensions="3 2"
Format="XML">
0 0
1 1
21346 1
1524 1
</DataStructure>
<DataStructure Format="HDF"
Dimensions="147495 1"
Dimensions="9405 1"
DataType="Float"
Precision="8">
Potential.h5:/potential.00000/Values
......
......@@ -12,10 +12,10 @@
<Time TimeType="Single" Value="0" />
<Topology
Type="Tetrahedron"
NumberOfElements="93511"
NumberOfElements="5041"
BaseOffset="0">
<DataStructure Format="HDF"
Dimensions="93511 4"
Dimensions="5041 4"
DataType="Int"
Precision="8">
Projection.h5:/Connections/Values
......@@ -23,19 +23,19 @@
</Topology>
<Geometry Type="X_Y_Z">
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
Projection.h5:/PointsX.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
Projection.h5:/PointsY.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
Projection.h5:/PointsZ.00000/Values
......@@ -48,16 +48,16 @@
Center="Node"
Name="projection">
<DataStructure ItemType="HyperSlab"
Dimensions="21346 3"
Dimensions="1524 3"
Type="HyperSlab">
<DataStructure Dimensions="3 2"
Format="XML">
0 0
1 1
21346 3
1524 3
</DataStructure>
<DataStructure Format="HDF"
Dimensions="147495 3"
Dimensions="9405 3"
DataType="Float"
Precision="8">
Projection.h5:/projection.00000/Values
......
......@@ -12,10 +12,10 @@
<Time TimeType="Single" Value="0" />
<Topology
Type="Tetrahedron"
NumberOfElements="93511"
NumberOfElements="5041"
BaseOffset="0">
<DataStructure Format="HDF"
Dimensions="93511 4"
Dimensions="5041 4"
DataType="Int"
Precision="8">
SheetsDirection.h5:/Connections/Values
......@@ -23,19 +23,19 @@
</Topology>
<Geometry Type="X_Y_Z">
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
SheetsDirection.h5:/PointsX.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
SheetsDirection.h5:/PointsY.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="21346"
Dimensions="1524"
DataType="Float"
Precision="8">
SheetsDirection.h5:/PointsZ.00000/Values
......@@ -48,16 +48,16 @@
Center="Node"
Name="sheets">
<DataStructure ItemType="HyperSlab"
Dimensions="21346 3"
Dimensions="1524 3"
Type="HyperSlab">
<DataStructure Dimensions="3 2"
Format="XML">
0 0
1 1
21346 3
1524 3
</DataStructure>
<DataStructure Format="HDF"
Dimensions="147495 3"
Dimensions="9405 3"
DataType="Float"
Precision="8">
SheetsDirection.h5:/sheets.00000/Values
......
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "FiberDirection.xdmf" [
<!ENTITY DataFile "FiberDirection.h5">
]>
<!-- FiberDirection.h5 is generated by LifeV -->
<Xdmf>
<Domain Name="FiberDirection">
<Grid Name="FiberDirectionGrid" GridType="Collection" CollectionType="Temporal">
<!-- Time 0 Iteration 00000 -->
<Grid Name="Mesh 0">
<Time TimeType="Single" Value="0" />
<Topology
Type="Tetrahedron"
NumberOfElements="317773"
BaseOffset="0">
<DataStructure Format="HDF"
Dimensions="317773 4"
DataType="Int"
Precision="8">
FiberDirection.h5:/Connections/Values
</DataStructure>
</Topology>
<Geometry Type="X_Y_Z">
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
FiberDirection.h5:/PointsX.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
FiberDirection.h5:/PointsY.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
FiberDirection.h5:/PointsZ.00000/Values
</DataStructure>
</Geometry>
<Attribute
Type="Vector"
Center="Node"
Name="fibers">
<DataStructure ItemType="HyperSlab"
Dimensions="62287 3"
Type="HyperSlab">
<DataStructure Dimensions="3 2"
Format="XML">
0 0
1 1
62287 3
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287 3"
DataType="Float"
Precision="8">
FiberDirection.h5:/fibers.00000/Values
</DataStructure>
</DataStructure>
</Attribute>
</Grid>
</Grid>
</Domain>
</Xdmf>
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Potential.xdmf" [
<!ENTITY DataFile "Potential.h5">
]>
<!-- Potential.h5 is generated by LifeV -->
<Xdmf>
<Domain Name="Potential">
<Grid Name="PotentialGrid" GridType="Collection" CollectionType="Temporal">
<!-- Time 0 Iteration 00000 -->
<Grid Name="Mesh 0">
<Time TimeType="Single" Value="0" />
<Topology
Type="Tetrahedron"
NumberOfElements="317773"
BaseOffset="0">
<DataStructure Format="HDF"
Dimensions="317773 4"
DataType="Int"
Precision="8">
Potential.h5:/Connections/Values
</DataStructure>
</Topology>
<Geometry Type="X_Y_Z">
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
Potential.h5:/PointsX.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
Potential.h5:/PointsY.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
Potential.h5:/PointsZ.00000/Values
</DataStructure>
</Geometry>
<Attribute
Type="Scalar"
Center="Node"
Name="potential">
<DataStructure ItemType="HyperSlab"
Dimensions="62287 1"
Type="HyperSlab">
<DataStructure Dimensions="3 2"
Format="XML">
0 0
1 1
62287 1
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287 1"
DataType="Float"
Precision="8">
Potential.h5:/potential.00000/Values
</DataStructure>
</DataStructure>
</Attribute>
</Grid>
</Grid>
</Domain>
</Xdmf>
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Projection.xdmf" [
<!ENTITY DataFile "Projection.h5">
]>
<!-- Projection.h5 is generated by LifeV -->
<Xdmf>
<Domain Name="Projection">
<Grid Name="ProjectionGrid" GridType="Collection" CollectionType="Temporal">
<!-- Time 0 Iteration 00000 -->
<Grid Name="Mesh 0">
<Time TimeType="Single" Value="0" />
<Topology
Type="Tetrahedron"
NumberOfElements="317773"
BaseOffset="0">
<DataStructure Format="HDF"
Dimensions="317773 4"
DataType="Int"
Precision="8">
Projection.h5:/Connections/Values
</DataStructure>
</Topology>
<Geometry Type="X_Y_Z">
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
Projection.h5:/PointsX.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
Projection.h5:/PointsY.00000/Values
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287"
DataType="Float"
Precision="8">
Projection.h5:/PointsZ.00000/Values
</DataStructure>
</Geometry>
<Attribute
Type="Vector"
Center="Node"
Name="projection">
<DataStructure ItemType="HyperSlab"
Dimensions="62287 3"
Type="HyperSlab">
<DataStructure Dimensions="3 2"
Format="XML">
0 0
1 1
62287 3
</DataStructure>
<DataStructure Format="HDF"
Dimensions="62287 3"
DataType="Float"
Precision="8">
Projection.h5:/projection.00000/Values
</DataStructure>
</DataStructure>
</Attribute>
</Grid>
</Grid>
</Domain>
</Xdmf>
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment