Hello,
I defined a region as the RegionIntersection of a thin and large cylinder and a shell and I want to know how many nodes of a lattice are contained inside.
My lattice is defined by 3 vectors of different lengths and one angle is 99 deg (monoclinic lattice).
I have no idea where to start. Any idea ?
EDIT 1 :
I created a list of the node coordinates :
d = 5;
noeuds = Table[
h*b1 + k*b2 + l*b3, {h, -d, d, 1}, {k, -d, d, 1}, {l, -d, d, 1}];
I will find a way to test if they are in the region.
EDIT 2 :
I finally used the fonction Select, to extract the points that are contained into the region of interest :
noeuds2 = Flatten[noeuds, 2];
Length[noeuds2]
bonsnoeuds = Select[noeuds2, 0.99*lk <= Norm[#] <= lk*1.01 &];
Length[bonsnoeuds]
Now I want to know if those selected points are contained inbetween 2 plans. I used a cylinder to represent the volume in-between those 2 plans but it seems it is not the fastest method as my results don't arrive instantly. I will have to repeat those calculations fore about 4k different plans hehe.