-
Notifications
You must be signed in to change notification settings - Fork 155
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Feature/topcharge #265
Feature/topcharge #265
Conversation
Provides "ext_topologicalchargelattice" extension that computes the topological charge using a lattice-based formalism, rather than the discretisation of the continuum form.
Re-introduced cell volumes to obtain true density for the topological charge. Comparable to output from ext_topologicalcharge
Convex corners are now treaded correctly in topologicalchargelattice. topologicalchargefinitelattice is now redundant, and hence removed.
@JeroenMulkers How does it score accuracy-wise in tests? |
It gives exactly the same results as Joo-Von's implementation in the unit tests. These unit tests have a trivial geometry, in which case both implementations are equivalent, and hence we expect that it yields the same results The only thing which is different here, is that the contribution of a triangle is counted twice in a convex corner of the geometry. |
@JeroenMulkers Good. Please give credit to the Joo-Von in the sources and make sure we ask to cite his paper if feature is used. |
This is an alternative for the
topologicalchargelattice
implementation proposed in pull request #262.The local topological charge is computed by taking a weighted sum of the topological charge contributions of the four elementary triangles defined by the four nearest neighbours.
In general, the elementary triangles of neighbouring cells overlap, in which case the charge contributions are counted twice. This is solved by using a weight of 1/2. (Exactly the same is done in the topologicalchargelattice implementation in #262)
However, the triangle at a convex corner of the geometry does not overlap any other triangles. Hence, at convex corners we use a weight of 1.