operators/indexselect.py
# ----------------------------------------------------------------------------
#  CLASSES: nightly
#
#  Test Case:  indexselect.py
#
#  Tests:      mesh      - 3D curvilinear, single domain,
#                          3D rectilinear, single domain.
#                          2D rectilinear, single domain.
#                          3D rectilinear, single domain.
#              plots     - pc, mesh, subset, vector, surface
#              operators - index select (with threshold and sphere slice)
#
#  Defect ID:  none
#
#  Programmer: Hank Childs
#  Date:       October 2, 2002
#
#  Modifications:
#    Kathleen Bonnell, Thu Aug 28 14:34:57 PDT 2003
#    Remove compound var name from subset plots.
#
#    Kathleen Bonnell, Wed Sep  3 09:31:25 PDT 2003
#    Opaque mode for MeshPlot is ON by default now, so turn it off
#    for this test.
#
#    Kathleen Bonnell, Thu Sep  4 11:10:11 PDT 2003
#    Change the way MeshPlot's opaque mode is set (now an enum, not a bool).
#
#    Kathleen Bonnell, Wed Sep  8 11:29:10 PDT 2004 
#    Add two new tests, testing 'modulo' problems. ('3308).
#
#    Kathleen Bonnell, Wed Jul 20 13:32:44 PDT 2005 
#    Add two new tests, testing 'group' problems. ('6386).
#
#    Kathleen Bonnell, Wed Aug 10 10:13:56 PDT 2005 
#    Add tests 27-30 which tests extracting faces along maximum and minimum
#    boundaries of 3d data. ('6106).
#
#    Kathleen Bonnell, Fri Aug 19 15:45:46 PDT 2005 
#    Add tests 31-32 which tests domain boundaries w/i a group that coincide
#    with external boundaries of the index-selection. ('6428)  
#
#    Kathleen Bonnell, Mon Jan 30 16:40:30 PST 2006 
#    Add tests 33-38 which tests point meshes.
#
#    Kathleen Bonnell, Thu Jun  7 14:43:32 PDT 2007 
#    IndexSelect atts have changed, make corrections. 
#
#    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, Wed Jan 11 13:57:00 PST 2012
#    Added test for Index Select of point meshes with Box operator applied.
#
#    Kathleen Biagas, Mon Dec 19 15:45:38 PST 2016
#    Use FilledBoundary plot for materials instead of Subset.
#
#    Eddie Rusu, Tue Apr 13 12:08:59 PDT 2021
#    Changed Pseudocolor CurrentPlot to ActualData.
#
# ----------------------------------------------------------------------------


OpenDatabase(silo_data_path("rect3d.silo"))


atts = IndexSelectAttributes()
# Dim is an enumerated type.  0 -> 1D, 1 -> 2D, 2 -> 3D.
atts.dim = 2
atts.xMin = 5
atts.xMax = 23
atts.xIncr = 2
atts.yMin = 0
atts.yMax = 20
atts.yIncr = 1
atts.zMin = 10
atts.zMax = -1
atts.zIncr = 2

# 3D, rectilinear, zonal var.
AddPlot("Pseudocolor", "d")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
view = GetView3D()
view.viewNormal = (-0.5, 0.5, 0.7172)
SetView3D(view)
Test("ops_indexselect01")
DeleteAllPlots()

# 3D, rectilinear, nodal var.
AddPlot("Pseudocolor", "u")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect02")
DeleteAllPlots()

# 3D, rectilinear, mesh.
m = MeshAttributes()
m.opaqueMode = m.Off
SetDefaultPlotOptions(m)
AddPlot("Mesh", "quadmesh3d")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect03")
DeleteAllPlots()

# 3D, rectilinear, subset.
AddPlot("FilledBoundary", "mat1")
AddOperator("IndexSelect")
atts.xIncr = 1
atts.yIncr = 1
atts.zIncr = 1
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect04")
DeleteAllPlots()

# 3D, rectilinear, after a threshold
AddPlot("Pseudocolor", "p")
a2 = ThresholdAttributes()
a2.listedVarNames = ("u")
a2.lowerBounds = (0.25)
a2.upperBounds = (1.0)
AddOperator("Threshold")
SetOperatorOptions(a2)
atts.zMin = 5
atts.zMax = 20
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect05")
DeleteAllPlots()

OpenDatabase(silo_data_path("curv3d.silo"))


# 3D, curvilinear, mesh plot
AddPlot("Mesh", "curvmesh3d")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
ResetView()
view = GetView3D()
view.viewNormal = (-0.5, 0.5, 0.7172)
SetView3D(view)
Test("ops_indexselect06")
DeleteAllPlots()

# 3D, curvilinear, PC plot
AddPlot("Pseudocolor", "d")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect07")
DeleteAllPlots()

# 3D, curvilinear, with sphere slice
AddPlot("Pseudocolor", "d")
AddOperator("SphereSlice")
a3 = SphereSliceAttributes()
a3.origin = (2.5, -0.5, 13)
a3.radius = 7
SetOperatorOptions(a3)
AddOperator("IndexSelect")
atts.xMin = 0
atts.xMax = -1
atts.yMin = 0
atts.yMax = -1
atts.zMin = 5
atts.zMax = 18
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect08")
DeleteAllPlots()

OpenDatabase(silo_data_path("rect2d.silo"))


# 2D, rectilinear, mesh plot
AddPlot("Mesh", "quadmesh2d")
AddOperator("IndexSelect")
# Dim is an enumerated type.  0 -> 1D, 1 -> 2D, 2 -> 3D.
atts.dim = 1
atts.xMin = 5
atts.xMax = 25
atts.yMin = 5
atts.yMax = 30
SetOperatorOptions(atts)
DrawPlots()
ResetView()
Test("ops_indexselect09")
DeleteAllPlots()

# 2D, rectilinear, PC plot
AddPlot("Pseudocolor", "d")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect10")
DeleteAllPlots()

OpenDatabase(silo_data_path("curv2d.silo"))


# 2D, curvilinear, mesh plot
AddPlot("Mesh", "curvmesh2d")
AddOperator("IndexSelect")
atts.yMax = 25
SetOperatorOptions(atts)
DrawPlots()
ResetView()
Test("ops_indexselect11")
DeleteAllPlots()

# 2D, curv, PC plot
AddPlot("Pseudocolor", "d")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect12")
DeleteAllPlots()

# 2D, curv, with threshold.
AddPlot("Pseudocolor", "d")
AddOperator("Threshold")
a2.listedVarNames = ("p")
a2.lowerBounds = (0.34)
a2.upperBounds = (0.5)
SetOperatorOptions(a2)
AddOperator("IndexSelect")
atts.xMin = 10
atts.xMax = -1
atts.yMin = 0
atts.yMin = -1
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect13")
DeleteAllPlots()

# 2D, curv of vector
AddPlot("Vector", "vel")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
Test("ops_indexselect14")
DeleteAllPlots()

# 2D surface of curvilinear
AddPlot("Pseudocolor", "d")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
AddOperator("Elevate")
DrawPlots()
ResetView()
Test("ops_indexselect15")
DeleteAllPlots()

# 3D multi-rect 
# test 'modulo' problems, bug: VisIt00003308
OpenDatabase(silo_data_path("multi_rect3d.silo"))

AddPlot("Pseudocolor", "d")
AddOperator("IndexSelect")
atts.dim = 2
atts.xMin = 0
atts.xMax = 1000
atts.yMin = 0
atts.yMax = 1000
atts.zMin = 0
atts.zMax = 1000
atts.xIncr = 3
atts.yIncr = 3
atts.zIncr = 3
SetOperatorOptions(atts)
DrawPlots()
ResetView()
Test("ops_indexselect16")

atts.xMin = 1
atts.xMax = 8
atts.yMin = 1
atts.yMax = 8
atts.zMin = 1
atts.zMax = 8
atts.xIncr = 5
atts.yIncr = 5
atts.zIncr = 5
SetOperatorOptions(atts)

Test("ops_indexselect17")
DeleteAllPlots()

# bigsil 
# test 'domainIndex' problems, bug: VisIt00005559
OpenDatabase(silo_data_path("bigsil.silo"))

AddPlot("Subset", "domains")
AddPlot("Mesh", "mesh")
AddOperator("IndexSelect")
atts.dim = atts.ThreeD
atts.xMin = 0
atts.xMax = 1000
atts.yMin = 0
atts.yMax = 1000
atts.zMin = 0
atts.zMax = 1000
atts.xIncr = 1
atts.yIncr = 1
atts.zIncr = 1
atts.useWholeCollection = 0
atts.categoryName = "domains"
atts.subsetName = "domain2"
SetOperatorOptions(atts)
DrawPlots()
ResetView()

v = GetView3D()
v.viewNormal = (0.653901, -0.39603, -0.64465)
v.focus = (0.5, 0.5, 0.5)
v.viewUp = (-0.251385, -0.917378, 0.308584)
v.viewAngle = 30
v.parallelScale = 0.866025
v.nearPlane = -1.73205
v.farPlane = 1.73205
v.imagePan = (0, 0)
v.imageZoom = 1
v.perspective = 1
v.eyeAngle = 2
v.centerOfRotationSet = 0
v.centerOfRotation = (0, 0, 0)
SetView3D(v)

Test("ops_indexselect18")

# smallest domainIndex should be 1
atts.subsetName = "domain1"
SetOperatorOptions(atts)
Test("ops_indexselect19")

atts.subsetName = "domain27"
SetOperatorOptions(atts)
Test("ops_indexselect20")

DeleteAllPlots()

#test whichData = OneGroup
AddPlot("Subset", "blocks")
s = SubsetAttributes()
s.wireframe = 1
s.lineWidth = 3
SetPlotOptions(s)
AddPlot("Mesh", "mesh")
AddPlot("Pseudocolor", "dist")
AddOperator("IndexSelect")
atts.dim = atts.ThreeD
atts.xMin = 1
atts.xMax = 2
atts.yMin = 1
atts.yMax = 2
atts.zMin = 0
atts.zMax = -1
atts.categoryName = "blocks"
atts.subsetName = "block0"
SetOperatorOptions(atts)

DrawPlots()
ResetView()
v = GetView3D()
v.viewNormal = (0, 0, -1)
v.focus = (0.5, 0.5, 0.5)
v.parallelScale = 0.866025
v.viewUp = (0, 1, 0)
v.nearPlane = -1.73205
v.farPlane  = 1.73205
v.perspective = 0
SetView3D(v)
v.imageZoom = 5.39496
SetView3D(v)
v.imagePan = (-0.23029, 0.259159)
SetView3D(v)
Test("ops_indexselect21")

atts.categoryName = "domains"
atts.subsetName = "domain1"
SetOperatorOptions(atts)
Test("ops_indexselect22")

SetActivePlots(1)
DeleteActivePlots()

SetActivePlots(1)
atts.categoryName = "blocks"
atts.subsetName = "block2"
atts.xMin = 0
atts.xMax = -1
atts.yMin = 5
atts.yMax = 19
atts.zMin = 7
atts.zMax = 8
SetOperatorOptions(atts)
ResetView()
Test("ops_indexselect23")

TurnMaterialsOff(("1", "5", "8"))
Test("ops_indexselect24")

TurnMaterialsOn()
TurnDomainsOff(("domain23"))
Test("ops_indexselect25")

TurnMaterialsOff(("1", "5", "8"))
Test("ops_indexselect26")

TurnMaterialsOn()
TurnDomainsOn()

DeleteAllPlots()

OpenDatabase(silo_data_path("rect3d.silo"))

# 3D, rectilinear, zonal var, index-selected at max/min boundaries.
AddPlot("Pseudocolor", "d")
AddOperator("IndexSelect")
atts.dim = atts.ThreeD
atts.xMin = 0
atts.xMax = -1 # maximum
atts.yMin = 40 # maximum boundary
atts.yMax = 40 # maximum boundary
atts.zMin = 0
atts.zMax = -1
atts.xIncr = 1
atts.yIncr = 1
atts.zIncr = 1
atts.useWholeCollection = 1
SetOperatorOptions(atts)
pc = PseudocolorAttributes()
pc.limitsMode = pc.ActualData
SetPlotOptions(pc)
DrawPlots()
ResetView()
view = GetView3D()
view.viewNormal = (0, 1, 0)
view.viewUp = (0, 0, -1)
view.perspective = 0
SetView3D(view)
Test("ops_indexselect27")

atts.yMin = 0
atts.yMax = 0
SetOperatorOptions(atts)
Test("ops_indexselect28")

DeleteAllPlots()


OpenDatabase(silo_data_path("curv3d.silo"))

# 3D, curvilinear, zonal var, index-selected at max/min boundaries.
AddPlot("Pseudocolor", "p")
AddOperator("IndexSelect")
atts.dim = atts.ThreeD
atts.xMin = 30
atts.xMax = 30
atts.yMin = 0
atts.yMax = -1
atts.zMin = 0
atts.zMax = -1
atts.xIncr = 1
atts.yIncr = 1
atts.zIncr = 1
SetOperatorOptions(atts)
pc = PseudocolorAttributes()
pc.limitsMode = pc.ActualData
SetPlotOptions(pc)
DrawPlots()
ResetView()
view = GetView3D()
view.viewNormal = (0, 1, 0)
view.viewUp = (0, 0, -1)
view.perspective = 0
SetView3D(view)
Test("ops_indexselect29")

atts.xMin = 0
atts.xMax = 0
SetOperatorOptions(atts)
Test("ops_indexselect30")

DeleteAllPlots()

#'6428, boundary between domains disappear when on edge of
# index-selected plot 
OpenDatabase(silo_data_path("bigsil.silo"))

AddPlot("Pseudocolor", "dist")
atts.dim = atts.ThreeD
atts.xMin = 0
atts.xMax = -1 #max
atts.yMin = 8
atts.yMax = 17
atts.zMin = 0
atts.zMax = -1 #max
atts.xIncr = 1
atts.yIncr = 1
atts.zIncr = 1
atts.useWholeCollection = 0
atts.categoryName = "blocks"
# middle group for this dataset, has group boundaries and internal
# domain boundaries
atts.subsetName = "block1"

AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
ResetView()
v = GetView3D()
v.viewNormal = (0.437142, -0.719538, -0.539603)
v.viewUp = (-0.309012, -0.683599, 0.661214)
v.perspective = 0
SetView3D(v)
Test("ops_indexselect31")

DeleteAllPlots()

# Filled boundary takes different path through IndexSelect than
# does PC, so test it too. 
AddPlot("FilledBoundary", "mat")
AddOperator("IndexSelect")
SetOperatorOptions(atts)
DrawPlots()
SetView3D(v)
Test("ops_indexselect32")
DeleteAllPlots()


# Test point meshes
TestSection("IndexSelect on Point meshes.")
OpenDatabase(data_path("vtk_test_data/ugrid_points.vtk"))

AddPlot("Pseudocolor", "u")
p = PseudocolorAttributes()
p.pointType = p.Box
p.pointSize = 0.2
SetPlotOptions(p)
DrawPlots()

v = GetView2D()
v.windowCoords = (-0.4, 9.4, -0.4, 9.4)
SetView2D(v)

Test("ops_indexselect33")

AddOperator("IndexSelect")
isa = IndexSelectAttributes()
isa.dim = isa.OneD
isa.xMin = 60
isa.xMax = 70
isa.xIncr = 1
SetOperatorOptions(isa)
DrawPlots()
Test("ops_indexselect34")

isa.xMin = 0
isa.xMax = -1
isa.xIncr = 2
SetOperatorOptions(isa)

Test("ops_indexselect35")

isa.xIncr = 3
SetOperatorOptions(isa)
Test("ops_indexselect36")

isa.xMin = 4
isa.xMax = -1
isa.xIncr = 2
SetOperatorOptions(isa)

Test("ops_indexselect37")

isa.xMin = 10
isa.xMax = 89
isa.xIncr = 7
SetOperatorOptions(isa)

Test("ops_indexselect38")

DeleteAllPlots()

OpenDatabase(silo_data_path("noise.silo"))

AddPlot("Pseudocolor", "PointVar")
AddOperator("IndexSelect")
# make sure we are using default atts
isa = IndexSelectAttributes()
SetOperatorOptions(isa)
box = BoxAttributes()
box.minx = -5
box.maxx = 5
box.miny = -5
box.maxy = 5
box.minz = -5
box.maxz = 5
AddOperator("Box")
SetOperatorOptions(box)
DrawPlots()

Test("ops_indexselect39")

DeleteAllPlots()


Exit()