<?xml version="1.0" encoding="UTF-8"?>
<ChasteParameters xmlns="https://chaste.comlab.ox.ac.uk/nss/parameters/2017_1"
    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
    xsi:schemaLocation="https://chaste.comlab.ox.ac.uk/nss/parameters/2017_1 heart/src/io/ChasteParameters_2017_1.xsd">

    <!--
        Note: all unit attributes must be included, but their content cannot be changed.
        They exist for documentary purposes only.
        See https://chaste.cs.ox.ac.uk/trac/wiki/ChasteUnits for details
        on the units used by Chaste.
    -->

    <Simulation>
        <!-- 
            Problem definition 
        -->
        <SpaceDimension>3</SpaceDimension>
        <SimulationDuration unit="ms">50.0</SimulationDuration>
        <Domain>Mono</Domain>
        <IonicModels>
            <Default>
                <Hardcoded>FaberRudy2000</Hardcoded>
            </Default>
        </IonicModels>

        <!-- 
            Mesh definition 
        -->
        <Mesh unit="cm">
            <Slab x="4.0" y="0.1" z="2.0" inter_node_space="0.1"/>
        </Mesh>

        <!-- 
            Stimuli (as many <Stimulus> definitions as needed) 
           -->
        <Stimuli>
            <Stimulus>
                <!-- #1 -->
                <Strength unit="uA/cm^3">-25500.0</Strength>
                <Duration unit="ms">2.0</Duration>
                <Delay unit="ms">0</Delay>
                <Location unit="cm">
                    <Cuboid>
                        <LowerCoordinates x="0" y="0" z="0"/>
                        <UpperCoordinates x="0.1" y="0.1" z="2.0"/>
                    </Cuboid>
                </Location>
            </Stimulus>
        </Stimuli>

        <!-- 
            Ionic cell model heterogeneities (as many <CellHeterogeneity> definitions as needed) 
           -->
        <CellHeterogeneities>
            <CellHeterogeneity>
                <!-- #1 -->
                <ScaleFactorGks unit="dimensionless">0.462</ScaleFactorGks>
                <ScaleFactorIto unit="dimensionless">0.0</ScaleFactorIto>
                <ScaleFactorGkr unit="dimensionless">1.0</ScaleFactorGkr>
                <Location unit="cm">
                    <Cuboid>
                        <LowerCoordinates x="0" y="0" z="0"/>
                        <UpperCoordinates x="1.5" y="0.1" z="2.0"/>
                    </Cuboid>
                </Location>
            </CellHeterogeneity>

            <CellHeterogeneity>
                <!-- #2 -->
                <ScaleFactorGks unit="dimensionless">1.154</ScaleFactorGks>
                <ScaleFactorIto unit="dimensionless">0.85</ScaleFactorIto>
                <ScaleFactorGkr unit="dimensionless">1.0</ScaleFactorGkr>
                <Location unit="cm">
                    <Cuboid>
                        <LowerCoordinates x="1.5" y="0" z="0"/>
                        <UpperCoordinates x="2.5" y="0.1" z="2.0"/>
                    </Cuboid>
                </Location>
            </CellHeterogeneity>

            <CellHeterogeneity>
                <!-- #3 -->
                <ScaleFactorGks unit="dimensionless">1.154</ScaleFactorGks>
                <ScaleFactorIto unit="dimensionless">1.0</ScaleFactorIto>
                <ScaleFactorGkr unit="dimensionless">1.0</ScaleFactorGkr>
                <Location unit="cm">
                    <Cuboid>
                        <LowerCoordinates x="1.5" y="0" z="0"/>
                        <UpperCoordinates x="4.0" y="0.1" z="2.0"/>
                    </Cuboid>
                </Location>
            </CellHeterogeneity>
        </CellHeterogeneities>

        <!-- 
            Output 
        -->
        <OutputDirectory>ChasteResults</OutputDirectory>
        <OutputFilenamePrefix>SimulationResults</OutputFilenamePrefix>
        <OutputVisualizer meshalyzer="yes"/>
    </Simulation>

    <Physiological>
        <!--
            Conductivities
        -->
        <IntracellularConductivities longi="1.75" trans="1.75" normal="1.75" unit="mS/cm"/>
        <ExtracellularConductivities longi="7.0" trans="7.0" normal="7.0" unit="mS/cm"/>
        <SurfaceAreaToVolumeRatio unit="1/cm"> 1400 </SurfaceAreaToVolumeRatio>
        <Capacitance unit="uF/cm^2"> 1.0 </Capacitance>
    </Physiological>

    <Numerical>
        <TimeSteps ode="0.005" pde="0.01" printing="0.1" unit="ms"/>
        <KSPTolerances>
            <KSPRelative>1e-6</KSPRelative>
        </KSPTolerances>
        <KSPSolver>gmres</KSPSolver>
        <KSPPreconditioner>jacobi</KSPPreconditioner>
    </Numerical>

    <PostProcessing> </PostProcessing>
</ChasteParameters>