Test on Gradient Functionality

[cindytsai]

Test on Gradient Functionality

• We expect it to be failed, since this requires nearby grids. So yt will definitly ask for non local grids.
• Related Issue: #26

TODO

• [ ] ~Write on Document that this would fail.~ (Move to #23 )

[cindytsai]

Test Run

• yt version: 4.1.0
• python version: 3.8
• gamer example: gamer/example/test_problem/Hydro/MHD_OrszagTangVortex/
• Inline Script:

  def yt_inline():
    ds = yt.frontends.libyt.libytDataset()
    ds.add_gradient_fields(("gas", "density"))
    sz = yt.ProjectionPlot(ds, 'z', "density_gradient_x", center='c')
    if yt.is_root():
        sz.save()

• Result: Failed, caused by asking non-local grids.
• Error Message:

  File "./inline_script.py", line 30, in yt_inline_inputArg
  sz = yt.ProjectionPlot(ds, 'z', "density_gradient_x", center='c')
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/visualization/plot_window.py", line 1844, in __init__
  proj = ds.proj(
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/construction_data_containers.py", line 491, in __init__
  self.get_data(field)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/construction_data_containers.py", line 245, in get_data
  self._handle_chunk(chunk, fields, tree)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/construction_data_containers.py", line 565, in _handle_chunk
  d = chunk[field] * dl
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/data_containers.py", line 258, in __getitem__
  self.get_data(f)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/selection_objects/data_selection_objects.py", line 211, in get_data
  self._generate_fields(fields_to_generate)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/selection_objects/data_selection_objects.py", line 232, in _generate_fields
  fd = self._generate_field(field)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/data_containers.py", line 299, in _generate_field
  tr = self._generate_fluid_field(field)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/data_containers.py", line 316, in _generate_fluid_field
  rv = self._generate_spatial_fluid(field, ngt_exception.ghost_zones)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/data_containers.py", line 363, in _generate_spatial_fluid
  source=gz[field][ngz:-ngz, ngz:-ngz, ngz:-ngz],
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/data_containers.py", line 258, in __getitem__
  self.get_data(f)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/construction_data_containers.py", line 833, in get_data
  self._fill_fields(fill)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/construction_data_containers.py", line 1395, in _fill_fields
  for chunk in ls.data_source.chunks(fields, "io"):
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/selection_objects/data_selection_objects.py", line 98, in chunks
  self.get_data(fields)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/data_objects/selection_objects/data_selection_objects.py", line 195, in get_data
  read_fluids, gen_fluids = self.index._read_fluid_fields(
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/geometry/geometry_handler.py", line 225, in _read_fluid_fields
  fields_to_return = self.io._read_fluid_selection(
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/frontends/libyt/io.py", line 163, in _read_fluid_selection
  data_view = self._get_data_from_libyt(g, fname)
  File "/home/cindytsai/software/python/python3.8/lib/python3.8/site-packages/yt/frontends/libyt/io.py", line 174, in _get_data_from_libyt
  assert data_convert is not None, "This MPI rank does not have grid id [%s]" % grid.id
  P000 yt : [ERROR    ] 2021-09-08 20:24:24,570 AssertionError: This MPI rank does not have grid id [153]