assert test where pytest can see it

[capitalaslash]

[nschloe]

That fails, but it's a meshio issue. It should probably warn that it can't write cell_sets yet, and that one should convert it to int_data first.

[capitalaslash]

I understand that the write out is performed in meshio, but the fact that pygmsh is not capable of printing out physical (and geometrical) tags, something that is required for any gmsh file, is a serious limitation in my opinion. #423 is a solution that does not require to touch meshio. I tried some months ago to get this to work by modifying meshio, but I gave up due to being far more involved than going the #423 way.

[nschloe]

Is sets_to_int_data not good enough then?

[capitalaslash]

I tried that but it fails. Turns out the cell_sets are wrong: that is another problem with meshio. I am looking into it. If I use gmsh.write() the output is ok.

[nschloe]

Turns out the cell_sets are wrong: that is another problem with meshio. I am looking into it.

If I understand correctly, fixing this will fix your problem.