Skip to content
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

elliptic cylinder #33

Open
stevengj opened this issue Nov 13, 2018 · 1 comment · May be fixed by #64
Open

elliptic cylinder #33

stevengj opened this issue Nov 13, 2018 · 1 comment · May be fixed by #64

Comments

@stevengj
Copy link
Collaborator

stevengj commented Nov 13, 2018

From NanoComp/meep#606: A recent post from the mailing list involved a request for an elliptic cylinder to be added as a new geometry object. This would have to be added first to libctl. An elliptic cylinder would be a useful compliment to the existing Cylinder and Ellipsoid objects.

cc @HomerReid

@stevengj
Copy link
Collaborator Author

Basically, this will be a subtype of cylinder that has additional properties: r2 (other semiaxis) and e1 and e2 (the two ellipse axes).

Similar to a block object, the three vectors (r1e1,r2e2,axis) define a 3x3 matrix M. You can cache inv(M) in the object (just like for a block). To check whether a point is in the elliptical cylinder, you first multiply inv(M) * (point - center) to convert the displacement to the "cylinder coordinates". In this coordinate system everything becomes easy: it is just an ordinary circular cylinder oriented along the 3rd axis. (To compute the normal vector, you will have to rotate the resulting vector back to the original coordinate system.)

@smartalecH smartalecH linked a pull request Oct 6, 2022 that will close this issue
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging a pull request may close this issue.

1 participant