# Instrucciones:

# 0) Cargar las funciones en bravais.m
bravais

# 1) Ingresar un conjunto de vectores primitivos para la celda 
#    unidad en forma de matriz, utilizando vectores fila. Por
#    ejemplo:
Rsc  = [1.0, 0.0, 0.0; 0.0, 1.0, 0.0; 0.0, 0.0, 1.0]; # cubica simple
Rbcc = [0.5, 0.5,-0.5; 0.5,-0.5, 0.5;-0.5, 0.5, 0.5]; # cubica centrada en el cuerpo
Rfcc = [0.5, 0.5, 0.0; 0.0, 0.5, 0.5; 0.5, 0.0, 0.5]; # cubica centrada en las caras 
Rhex = [sqrt(3)/2, 0.5, 0.0; sqrt(3)/2, -0.5, 0.0; 0.0, 0.0, sqrt(8/3)]; # hexagonal c/a = 1.633 (ideal)

# 2) Ingresar un conjunto de vectores para base en forma de matriz,
#    utilizando vectores columna (si no requiere base, poner solo un 0) \n");

Sari = [0.0, 0.0, 0.0; 0.0, 0.5, 0.0; ...
        0.0, 0.0, 0.5; 0.5, 0.0, 0.0]; # + SC = cubica centrada en las aristas

Snacl = [0.0, 0.0, 0.0;  0.5, 0.0, 0.0]; # + FCC = NaCl

Shcp = [0.0, 0.0, 0.0; sqrt(3)/3, 0.0, sqrt(2/3)];

# 3) Ver la red directa (elegir la herramienta "rotar" para cambiar el angulo de vision)

verRed3d(Rbcc, 0, [3,3,3], 1); # ver red BCC en la figura 1

Sbcc = [0.0, 0.0, 0.0; 0.5, 0.5, 0.5];
verRed3d(Rsc, Sbcc, [3,3,3], 2); # compara la misma red con base en la figura 2

# 4) Comparar la recíproca de la BCC con la FCC directa

verRed3d(reciproco3d(Rbcc)/(2*pi), 0, [2,2,2], 3);
verRed3d(Rfcc, 0, [2,2,2], 4);

# 5) Calcular distancias a vecinos hasta 3 unidades de la red diamante (Fcc + base)

vecinos(Rfcc, [0.0, 0.0, 0.0; 0.25,0.25,0.25], [3,3,3])

# 6) Simular un difractograma de polvo de NaCl (Na+ = 10 e-, Cl- = 18 e-, a = 5.64 amstrongs)

xrdPolvo (Rfcc , [0,0,0;0.5,0.5,0.5], [18,10], [5,5,5], 4); # adimensionalizado
xrdPolvo (5.64 * Rfcc , 5.64 * [0,0,0;0.5,0.5,0.5], [18,10], [1,1,1], 5); # en Amstrongs