<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">my preference goes towards using numpy for this type to query<br>
<br>
use the where() function.<br>
<br>
# example 1:<br>
<br>
import vtk<br>
from vtk.numpy_interface import dataset_adapter as dsa<br>
<br>
# assuming your mesh data is in the object 'vtkgrid'<br>
np_grid = dsa.WrapDataObject(vtkgrid)<br>
<br>
# assuming your scalar variable is called 'scalars'<br>
np_scalars = np_grid.PointData['scalars']<br>
<br>
# min and max are:<br>
print np.min(np_scalars), np.max(np_scalars)<br>
<br>
#Their position(s) in the array is(are):<br>
np.where(np_scalars == np.max(np_scalars))<br>
np.where(np_scalars == np.min(np_scalars))<br>
<br>
# example 2:<br>
<br>
a = np.array([0, 0, 1, 2, 8, 9, 10, 9, 0, 0, 0, 9, 10, 9, 8, 6, 4, 3, 0, 0])<br>
<br>
<div>In [51]: np.where(a == a.min())<br>
Out[51]: (array([ 0, 1, 8, 9, 10, 18, 19]),)<br>
<br>
In [52]: np.where(a == a.max())<br>
Out[52]: (array([ 6, 12]),)<br>
<br>
References:<br>
<br>
<a href="http://kitware.com/blog/home/post/714" target="_blank">http://kitware.com/blog/home/post/714</a><br>
<br>
<div class="BodyFragment"><font size="2"><span style="font-size:10pt;">
<div class="PlainText">-----------------<br>
Jean<br>
CSCS<br>
</div>
</span></font></div>
</div>
</div>
</body>
</html>