Stencils with 2nd Neighborhood
von @kuron:
für ein bestimmtes Problem benötige ich einen Stencil mit mehr als einer Schicht Nachbarn, z.B. einen D3Q125. LBM will ich damit zum Glück nicht machen, nur den PDE-Löser benutzen, und ich brauche auch nur eine Untermenge der 125 Nachbarn. Da der Stencil-Generator von dir ist: wie lässt sich das am besten und möglichst minimal-invasiv umsetzen? Ich müsste wohl in Directions.h die 98 fehlenden Richtungen in alle Arrays schreiben und in generate.py die Richtungen und die Vektoren der elementaren Richtungen eintragen. Die wesentliche Frage: wie soll ich die neuen Richtungen benennen? (2,0,0) könnte man z.B. EE nennen — aber das erfordert z.B. Änderungen an isSubDirection() in generate.py.