add demo_0014_cube_varying_terial

[Aminofa70]

Dear @Kevin-Mattheus-Moerman

I implemented demo_0014_cube_varying_terial of Gibbon in Julia.
How can I share it with you.
Also I came across this problem.
In demo_0001, you have given
Fb_top = Fb[CFb_type.==1]
Fb_bottom = Fb[CFb_type.==2]

But when i use this for multi-material, the top and bottom are reversed.

With this it is correct
Fb_top = Fb[CFb_type.==2]
Fb_bottom = Fb[CFb_type.==1]

I am I wrong?
please refer to the following; bcPrescribeList_z is applied at the bottom

Node sets

bcPrescribeList_z = "bcPrescribeList_z"
bcSupportList_x = "bcSupportList_x"
bcSupportList_y = "bcSupportList_y"
bcSupportList_z = "bcSupportList_z"
aen(Mesh_node,"NodeSet",join([@sprintf("%i",x) for x ∈ elements2indices(Fb_bottom)],','); name=bcPrescribeList_z)
aen(Mesh_node,"NodeSet",join([@sprintf("%i",x) for x ∈ elements2indices(Fb_s1)],','); name=bcSupportList_x)
aen(Mesh_node,"NodeSet",join([@sprintf("%i",x) for x ∈ elements2indices(Fb_s2)],','); name=bcSupportList_y)
aen(Mesh_node,"NodeSet",join([@sprintf("%i",x) for x ∈ elements2indices(Fb_top)],','); name=bcSupportList_z)

Also I have utilized this:
E_youngs = 1
nelement = length(E)
E_youngs_elem = fill(E_youngs, nelement)
ν = 0.4
nu_elem = fill(ν, nelement)

Does Comodo.jl have: triplyPeriodicMinimal and patchCentre? to be similar to original MATLAB code
I checked the code with MATLAB and it works well for the aforementioned element data .
Thanks for your great work.

[Kevin-Mattheus-Moerman]

@Aminofa70 thanks. Sounds great. If you'd like to contribute then submitting a pull request would be the best way to do it.

[Aminofa70]

Dear @Kevin-Mattheus-Moerman
I think the code is not complete and similar to Gibbon, so I look forward to seeing your own code.

Cheers