Nico Schlömer

I really don't get it. I think the data is read alright, but I don't understand why the cell data is converted into the wrong endianness now. Can I get...

Okay, did apply.

This has to go: ``` cells = np.asarray(cells,dtype='object') ``` optimesh only operates on triangular cells, so you need to extract those first.

> optimesh only operates on triangular cells, so you need to extract those first.

You need to provide a numpy array with the triangle connectivity.

Thanks for the report. I've tried to reproduce this, but the script just ran forever, filling the memory until it eventually choked. Somethings not right with `Union`. I'll check it...

I've improved some things in domain specifications and released a new version, but the above mesh will still fail. I found that all cuboid unions fail where the two cuboids...

Yup, feature edges aren't preserved by default. Right now, _all_ primitive domains are implemented in terms of level-sets, even polygonal domains like cuboids. As a consequence, feature edges are _not_...

This really happens? Seeing is believing, please attach a surface mesh so I can try it out for myself. :smile_cat: I can try and fix it then.

You can zip the file and attach it.