# Create an orthorombic (default) cell
uc <- unit_cell(10,30,15)
# Create the related reciprocal cell
ruc <- create_rec_unit_cell(uc)
# Create symmetry (P n c 2)
csym <- cryst_symm(30)
# Create a few records (these include syst. absences)
records <- expand.grid(H=-2:2,K=-2:2,L=-2:2)
print(length(records[,1]))
# dtypes are all H
dtypes <- c("H","H","H")
# Create merged_reflections object with H, K, L
# Systematic absences have been eliminated
mrefs <- merged_reflections(ruc,csym,records,dtypes)
print(length(mrefs$records[,1]))
Run the code above in your browser using DataLab