operators/onionpeel.py
# ----------------------------------------------------------------------------
#  CLASSES: nightly
#
#  Test Case:  onionpeel.py
#
#  Tests:      mesh      - 3D rectilinear, mulitple domain, mulitple groups
#              plots     - subset
#              operators - onionpeel
#              selection - none
#
#  Defect ID:  
#
#  Programmer: Kathleen Bonnell 
#  Date:       December 23, 2003 
#
#  Modifications:
#    Kathleen Bonnell, Thu Feb 26 15:31:31 PST 2004
#    Added TestUCD.
#
#    Kathleen Bonnell, Wed Dec 15 15:24:50 PST 2004 
#    Added TestGlobalZoneId.
#
#    Kathleen Bonnell, Tue Jan 18 19:34:41 PST 2005 
#    Added TestMatSelect.
#
#    Kathleen Bonnell, Wed Jan 19 15:45:38 PST 2005 
#    Added TestNodeId.
#
#    Kathleen Bonnell, Tue Jun 14 11:33:39 PDT 2005 
#    Added TestFilledBoundary.
#
#    Kathleen Bonnell, Thu Jul 21 09:09:26 PDT 2005 
#    Modified indices used in TestBigSil, to match corrsponding changes
#    to bigsil.silo. 
#
#    Kathleen Bonnell, Fri Aug 19 16:23:45 PDT 2005 
#    Modified indices used in TestBigSil, more changes to bigsil.silo. 
#
#    Kathleen Bonnell, Thu Sep 22 11:59:42 PDT 2005 
#    Added TestBoundary. 
#
#    Mark C. Miller, Wed Jan 20 07:37:11 PST 2010
#    Added ability to swtich between Silo's HDF5 and PDB data.
#
#    Cyrus Harrison, Thu Mar 25 09:57:34 PDT 2010
#    Added call(s) to DrawPlots() b/c of changes to the default plot state
#    behavior when an operator is added.
#
#    Kathleen Biagas, Fri Jul 25 15:10:11 MST 2014
#    Added ArbPoly test, which demonstrates use of new att: honorOriginalMesh.
#
#    Kathleen Biagas, Mon Mar 23 17:25:07 PDT 2015
#    Modified AMR test to include the default Subset plot, now that it
#    is fixed.
#
#    Kathleen Biagas, Tue Nov 29 10:40:50 PST 2022
#    Change use of obosolete text settings to new textFont settings.
#
# ----------------------------------------------------------------------------

def TestBigSil():
    OpenDatabase(silo_data_path("bigsil.silo"))


    AddPlot("Subset", "domains")
    AddPlot("Subset", "blocks")
    s = SubsetAttributes()
    s.legendFlag = 0
    s.wireframe = 1
    s.colorType = s.ColorBySingleColor
    s.singleColor = (100,100,100,255)
    SetPlotOptions(s)

    SetActivePlots(0)
    AddOperator("OnionPeel")
    op = OnionPeelAttributes()
    op.categoryName = "blocks"
    op.subsetName = "block0"
    op.index = (4, 11, 5)
    op.logical = 1
    op.adjacencyType = op.Face
    op.requestedLayer = 3
    SetOperatorOptions(op)
    DrawPlots()

    Test("ops_onionpeel_01")

    op.subsetName = "block1"
    op.index = (18, 4, 5)
    SetOperatorOptions(op)

    Test("ops_onionpeel_02")

    op.subsetName = "block2"
    op.index = (11, 18, 5)
    SetOperatorOptions(op)

    Test("ops_onionpeel_03")

    DeleteAllPlots()

def TestUCD():
    # This is to test that onion peel will work (by default)
    # on the first "On" set -- in this case "domain4".
    OpenDatabase(silo_data_path("multi_ucd3d.silo"))

    AddPlot("Pseudocolor", "d")
    domains = ("domain4", "domain5", "domain6")
    TurnDomainsOff()
    TurnDomainsOn(domains)
    DrawPlots()
    ResetView()
    SetViewExtentsType("actual")
    AddOperator("OnionPeel")
    DrawPlots()
    Test("ops_onionpeel_04")

    DeleteAllPlots()

def TestGlobalZoneId():
    # This is to test that onion peel will work when 
    # a global cell number is specified. 
    OpenDatabase(silo_data_path("global_node.silo"))

    TurnDomainsOn()
    AddPlot("Subset", "domains")
    s = SubsetAttributes()
    s.legendFlag = 0
    s.wireframe = 1
    s.colorType = s.ColorBySingleColor
    s.singleColor = (100,100,100,255)
    SetPlotOptions(s)
    AddPlot("Pseudocolor", "dist")
    AddOperator("OnionPeel")
    op = OnionPeelAttributes()
    op.categoryName = "domains"
    op.subsetName = "domain5"
    op.index = 372
    op.logical = 0
    op.adjacencyType = op.Face
    op.requestedLayer = 7
    SetOperatorOptions(op)

    DrawPlots()
    ResetView()
    SetViewExtentsType("original")

    Test("ops_onionpeel_05")

    #only change is to specify that seed cell is global
    op.useGlobalId = 1
    SetOperatorOptions(op)
    Test("ops_onionpeel_06")

    DeleteAllPlots()

def TestMatSelect():
    OpenDatabase(silo_data_path("rect3d.silo"))

    AddPlot("FilledBoundary", "mat1")
    AddOperator("OnionPeel")
    op = OnionPeelAttributes()
    op.index = 8
    op.requestedLayer = 5
    op.adjacencyType = op.Face
    SetOperatorOptions(op)
    DrawPlots()
    SetViewExtentsType("actual")
    Test("ops_onionpeel_07")

    TurnMaterialsOff("1")
    Test("ops_onionpeel_08")

    DeleteAllPlots()

    OpenDatabase(silo_data_path("globe.silo"))

    AddPlot("Pseudocolor", "u")
    TurnMaterialsOff("4")
    AddOperator("OnionPeel")
    op.requestedLayer = 2
    op.index = 27
    SetOperatorOptions(op)
    DrawPlots()
    Test("ops_onionpeel_09")
    TurnMaterialsOn()
    DeleteAllPlots()

def TestNodeId():
    OpenDatabase(silo_data_path("ucd2d.silo"))

    AddPlot("Mesh", "ucdmesh2d")
    AddPlot("Pseudocolor", "d")
    AddOperator("OnionPeel")
    op = OnionPeelAttributes()
    op.index = 0
    op.seedType = op.SeedCell
    SetOperatorOptions(op)
    DrawPlots()
    Test("ops_onionpeel_10")

    op.seedType = op.SeedNode
    SetOperatorOptions(op)
    Test("ops_onionpeel_11")

    op.requestedLayer = 2
    op.adjacencyType = op.Face
    SetOperatorOptions(op)
    Test("ops_onionpeel_12")

    DeleteAllPlots()
    OpenDatabase(silo_data_path("rect3d.silo"))

    AddPlot("Mesh", "quadmesh3d")
    AddPlot("Pseudocolor", "d")
    SetActivePlots((0, 1))
    TurnMaterialsOff("1")
    AddOperator("OnionPeel")
    op.index = 44
    op.requestedLayer = 0
    SetOperatorOptions(op)
    DrawPlots()
    SetViewExtentsType("actual")
    v = GetView3D()
    v.viewNormal = (0, 0, -1)
    SetView3D(v)
    Test("ops_onionpeel_13")
    op.requestedLayer = 12
    SetOperatorOptions(op)
    Test("ops_onionpeel_14")
    TurnMaterialsOn()
    DeleteAllPlots()

def TestFilledBoundary():
    # need to test index on high-side
    OpenDatabase(silo_data_path("rect3d.silo"))

    AddPlot("Mesh", "quadmesh3d")
    AddPlot("FilledBoundary", "mat1")
    SetActivePlots((0,1))
    AddOperator("OnionPeel")
    op = OnionPeelAttributes()
    op.index = 15735
    op.requestedLayer = 5
    op.adjacencyType = op.Face
    SetOperatorOptions(op)
    DrawPlots()
    ResetView()
    SetViewExtentsType("actual")
    Test("ops_onionpeel_15")

    op.seedType = op.SeedNode
    op.index = 32713
    SetOperatorOptions(op)
    Test("ops_onionpeel_16")
    DeleteAllPlots()

    OpenDatabase(silo_data_path("globe.silo"))

    AddPlot("Mesh", "mesh1")
    AddPlot("FilledBoundary", "mat1")
    SetActivePlots((0,1))
    AddOperator("OnionPeel")
    op.seedType = op.SeedCell
    op.requestedLayer = 2
    op.index = 742
    SetOperatorOptions(op)
    DrawPlots()
    Test("ops_onionpeel_17")
    op.index = 1117
    SetOperatorOptions(op)
    DrawPlots()
    ResetView()
    Test("ops_onionpeel_18")

    op.seedType = op.SeedNode
    op.requestedLayer = 6
    op.index = 965
    SetOperatorOptions(op)
    Test("ops_onionpeel_19")
    DeleteAllPlots()

    OpenDatabase(silo_data_path("curv3d.silo"))

    AddPlot("Mesh", "curvmesh3d")
    AddPlot("FilledBoundary", "mat1")
    SetActivePlots((0, 1))
    AddOperator("OnionPeel")
    op.seedType = op.SeedCell
    op.requestedLayer = 5
    op.index = 19412
    SetOperatorOptions(op)
    DrawPlots()
    ResetView()
    v = GetView3D()
    v.viewNormal = (0, 0, 1)
    v.focus = (-2.29977, 0.966088, 16.5)
    v.viewUp = (0, 1, 0)
    v.viewAngle = 30
    v.parallelScale = 5.61745
    v.nearPlane = -11.2349
    v.farPlane = 11.2349
    v.imagePan = (-0.00484124, 0.00104964)
    v.imageZoom = 4.08886
    v.perspective = 1
    v.eyeAngle = 2
    v.centerOfRotationSet = 0
    v.centerOfRotation = (-2.29977, 0.966088, 16.5)
    SetView3D(v)
    Test("ops_onionpeel_20")
    DeleteAllPlots()


def TestBoundary():
    # need to test index on high-side
    OpenDatabase(silo_data_path("rect3d.silo"))

    AddPlot("Mesh", "quadmesh3d")
    AddPlot("Boundary", "mat1")
    SetActivePlots((0,1))
    AddOperator("OnionPeel")
    op = OnionPeelAttributes()
    op.index = 15735
    op.requestedLayer = 5
    op.adjacencyType = op.Face
    SetOperatorOptions(op)
    DrawPlots()
    ResetView()
    SetViewExtentsType("actual")
    Test("ops_onionpeel_21")

    op.seedType = op.SeedNode
    op.index = 1370
    SetOperatorOptions(op)
    Test("ops_onionpeel_22")
    DeleteAllPlots()

    OpenDatabase(silo_data_path("globe.silo"))

    AddPlot("Mesh", "mesh1")
    AddPlot("Boundary", "mat1")
    SetActivePlots((0,1))
    AddOperator("OnionPeel")
    op.seedType = op.SeedCell
    op.requestedLayer = 2
    op.index = 17
    SetOperatorOptions(op)
    DrawPlots()
    ResetView()
    v = GetView3D()
    v.viewNormal = (1, 0, 0)
    v.viewUp = (0, 1, 0)
    SetView3D(v)
    Test("ops_onionpeel_23")
    op.index = 117
    SetOperatorOptions(op)
    DrawPlots()
    ResetView()
    Test("ops_onionpeel_24")

    op.seedType = op.SeedNode
    op.requestedLayer = 6
    op.index = 13
    SetOperatorOptions(op)
    Test("ops_onionpeel_25")
    DeleteAllPlots()

    OpenDatabase(silo_data_path("curv3d.silo"))

    AddPlot("Mesh", "curvmesh3d")
    AddPlot("Boundary", "mat1")
    SetActivePlots((0, 1))
    AddOperator("OnionPeel")
    op.seedType = op.SeedCell
    op.requestedLayer = 5
    op.index = 19472
    SetOperatorOptions(op)
    DrawPlots()
    ResetView()
    v = GetView3D()
    v.viewNormal = (0, 0, 1)
    v.focus = (-2.29977, 0.966088, 16.5)
    v.viewUp = (0, 1, 0)
    v.viewAngle = 30
    v.parallelScale = 5.61745
    v.nearPlane = -11.2349
    v.farPlane = 11.2349
    v.imagePan = (-0.00484124, 0.00104964)
    v.imageZoom = 4.08886
    v.perspective = 1
    v.eyeAngle = 2
    v.centerOfRotationSet = 0
    v.centerOfRotation = (-2.29977, 0.966088, 16.5)
    SetView3D(v)
    Test("ops_onionpeel_26")
    DeleteAllPlots()

def TestArbPoly():
    TestSection("Arbitrary Polyhedra, honor original mesh")
    OpenDatabase(silo_data_path("arbpoly-zoohybrid.silo"))
    DefineScalarExpression("gzid", "global_zoneid(<3D/mesh1>)")
    AddPlot("Mesh", "3D/mesh1")
    AddPlot("Pseudocolor", "gzid")
    pc = PseudocolorAttributes()
    pc.colorTableName="levels"
    SetPlotOptions(pc)
    AddOperator("OnionPeel")
    DrawPlots()
    ResetView()
    op = OnionPeelAttributes()

    SetQueryOutputToObject()
    # there are only 11 zones in original mesh
    # when NumZonesQuery is fixed, can use it instead
    # SetActivePlots(0)
    # numZones = Query("NumZones", use_actual_data=0)['num_zones']
    # SetActivePlots(1)
    numZones = 11
    for i in range(numZones):
        op.index = i
        SetOperatorOptions(op)
        DrawPlots()
        Test("poly_originalMesh_%02d"%i)

    TestSection("Arbitrary Polyhedra, honor actual mesh")
    op.honorOriginalMesh = 0
    SetActivePlots(0)
    numZones = Query("NumZones", use_actual_data=1)['num_zones']
    print("numZones: ", numZones)
    SetActivePlots(1)
    for i in range (numZones):
        op.index = i
        SetOperatorOptions(op)
        DrawPlots()
        Test("poly_actualMesh%02d"%i)

    DeleteAllPlots()

def TestBigSilMesh():
    # use actual spatial extents
    SetViewExtentsType(1)

    # non-logical index
    OpenDatabase(silo_data_path("bigsil.silo"))
    AddPlot("Mesh", "mesh")
    mesh = MeshAttributes()
    mesh.showInternal = 1
    SetPlotOptions(mesh)

    AddPlot("Label", "mesh")
    label = LabelAttributes()
    label.showNodes = 1
    label.showCells = 1
    label.drawLabelsFacing = label.FrontAndBack
    label.labelDisplayFormat = label.Index
    label.textFont1.useForegroundColor = 0
    label.textFont1.color = (255, 0, 0, 0)
    label.textFont2.useForegroundColor = 0
    label.textFont2.color = (0, 0, 255, 0)
    label.depthTestMode = label.LABEL_DT_NEVER
    SetPlotOptions(label)

    AddOperator("OnionPeel", 1)
    op = OnionPeelAttributes()
    op.categoryName = "domains"
    op.subsetName = "domain14"
    op.logical = 0
    op.index = (211)
    op.requestedLayer = 0
    op.seedType = op.SeedNode
    SetOperatorOptions(op,0, 1);
    DrawPlots()
    Test("bigsil_mesh_01")

    op.seedType = op.SeedCell
    op.index = (223)
    op.requestedLayer = 1
    SetOperatorOptions(op,0, 1);
    DrawPlots()
    Test("bigsil_mesh_02")

    # Blocks/Domains, logical index
    op.logical = 1
    op.index = (3, 2, 3)
    op.requestedLayer = 0
    op.seedType = op.SeedNode
    SetOperatorOptions(op,0, 1);
    DrawPlots()
    Test("bigsil_mesh_03")

    op.seedType = op.SeedCell
    op.index = (6, 3, 4)
    op.requestedLayer = 1
    SetOperatorOptions(op,0, 1);
    DrawPlots()
    Test("bigsil_mesh_04")

    DeleteAllPlots()
    CloseDatabase(silo_data_path("bigsil.silo"))
    # use original spatial extents
    SetViewExtentsType(0)

def TestAMR():
    # AMR Mesh
    OpenDatabase(data_path("samrai_test_data/sil_changes/dumps.visit"))
    # use actual spatial extents
    SetViewExtentsType(1)

    AddPlot("Mesh", "amr_mesh")
    mesh = MeshAttributes()
    mesh.showInternal = 1
    SetPlotOptions(mesh)

    AddPlot("Label", "amr_mesh")
    label = LabelAttributes()
    label.showNodes = 1
    label.showCells = 1
    label.drawLabelsFacing = label.FrontAndBack
    label.labelDisplayFormat = label.Index
    label.textFont1.useForegroundColor = 0
    label.textFont1.color = (255, 0, 0, 0)
    label.textFont2.useForegroundColor = 0
    label.textFont2.color = (0, 0, 255, 0)
    label.depthTestMode = label.LABEL_DT_NEVER
    SetPlotOptions(label)

    AddOperator("OnionPeel", 1)
    op = OnionPeelAttributes()
    op.categoryName = "levels"
    op.subsetName = "level2"
    op.logical = 1
    op.index = (11, 15, 8)
    op.seedType = op.SeedNode
    SetOperatorOptions(op,0, 1);
    DrawPlots()
    Test("amr_3d_01")

    op.seedType = op.SeedCell
    op.index = (10, 15, 8)
    op.requestedLayer = 1
    SetOperatorOptions(op,0, 1);
    Test("amr_3d_02")

    DeleteAllPlots()
    CloseDatabase(data_path("samrai_test_data/sil_changes/dumps.visit"))

    # Test Vector plot
    OpenDatabase(data_path("samrai_test_data/mats-par3/dumps.visit"))

    AddPlot("Mesh", "amr_mesh")
    AddPlot("Vector", "Velocity")
    TimeSliderNextState()
    AddOperator("OnionPeel", 1)
    op = OnionPeelAttributes()
    op.categoryName = "levels"
    op.subsetName = "level1"
    op.logical = 1
    op.index = (12, 15, 10)
    op.seedType = op.SeedCell
    op.requestedLayer = 1
    SetOperatorOptions(op,0, 1);
    DrawPlots()
    Test("amr_3d_03")

    op.seedType = op.SeedNode
    op.index = (13, 15, 10)
    op.requestedLayer = 0
    SetOperatorOptions(op,0, 1);
    Test("amr_3d_04")

    DeleteAllPlots()
    CloseDatabase(data_path("samrai_test_data/mats-par3/dumps.visit"))

    # 2D AMR Mesh
    OpenDatabase(data_path("samrai_test_data/ale2d/dumps.visit"))

    # use actual spatial extents
    SetViewExtentsType(1)

    AddPlot("Mesh", "amr_mesh")
    AddPlot("Contour", "Density")

    AddOperator("OnionPeel", 1)
    op = OnionPeelAttributes()
    op.categoryName = "levels"
    op.subsetName = "level2"
    op.logical = 1
    op.index = (384, 490)
    op.seedType = op.SeedCell
    op.requestedLayer = 1
    SetOperatorOptions(op,0, 1);
    DrawPlots()
    Test("amr_2d_01")

    DeleteAllPlots()
    CloseDatabase(data_path("samrai_test_data/ale2d/dumps.visit"))

    OpenDatabase(data_path("AMRStitchCell_test_data/AMRStitchCellTest.1.no_ghost.phi.2d.hdf5"))
    AddPlot("Subset", "levels")
    AddPlot("Mesh", "Mesh")
    AddPlot("Contour", "phi")
    AddPlot("Label", "Mesh")
    label = LabelAttributes()
    label.showNodes = 1
    label.showCells = 1
    label.drawLabelsFacing = label.FrontAndBack
    label.labelDisplayFormat = label.Index
    label.textFont1.useForegroundColor = 0
    label.textFont1.color = (255, 0, 0, 0)
    label.textFont2.useForegroundColor = 0
    label.textFont2.color = (0, 0, 255, 0)
    label.depthTestMode = label.LABEL_DT_NEVER
    SetPlotOptions(label)

    AddOperator("OnionPeel", 1)
    op = OnionPeelAttributes()
    op.categoryName = "levels"
    op.subsetName = "level2"
    op.logical = 1
    op.seedType = op.SeedCell
    op.requestedLayer = 1
    op.index = (24, 35)
    SetOperatorOptions(op, 0, 1)
    DrawPlots()
    Test("amr_2d_02")
    # use original spatial extents
    SetViewExtentsType(0)

def Main():
    TestBigSil()
    TestUCD()
    TestGlobalZoneId()
    TestMatSelect()
    TestNodeId()
    TestFilledBoundary()
    TestBoundary()
    TestArbPoly()
    TestBigSilMesh()
    TestAMR()

Main()
Exit()