Generowanie współrzędnych trójkątnych / sześciokątnych (xyz)

Próbuję wymyślić funkcję iteracyjną, która generuje współrzędne xyz dla siatki sześciokątnej. Z początkową pozycją szesnastkową (powiedzmy dla prostoty 0,0,0...

Konwersja Pythona między współrzędnymi