incorporate-qhull

fcc_lattice_test.cc
View Options


			#include <iostream>
			#include <string>
			#include <vector>


			#include "cadd_mesher.hh"


			int main(int argc, char *argv[])
			{
			Real consts[3] = {1,2,3};
			Real miller2[4] = {1, 1, -1, 1};
			Real miller[9] = {1, 1, -2, -1, 1, 0, 1, 1, 1};
			std::vector<Real> vec_miller(miller, &miller[9]);
			std::vector<Real> vec_consts(consts, &consts[2]);
			FccLattice<threeD> vec_lattice_3d(vec_consts, vec_miller);
			std::cout << "Test for vector constructor: " << vec_lattice_3d << std::endl;

			FccLattice<twoD> lattice2d(consts, miller2);
			std::cout << lattice2d << std::endl;
			FccLattice<threeD> lattice3d(consts, miller);
			std::cout << lattice3d << std::endl;

			PointContainer<twoD> atoms("atoms");
			int lattice_cords[] = { 0, 0,
			1, 0,
			-1, 0,
			0, 1,
			0,-1,
			2, 0,
			-2, 0,
			0, 2,
			0,-2,
			3, 0,
			-3, 0,
			0, 3,
			0,-3,
			4, 0,
			-4, 0,
			0, 4,
			0,-4};
			for (Uint i = 0; i < 17 ; ++i) {
			lattice2d.fillAtoms(lattice_cords + 2*i, atoms);
			}
			std::cout << atoms;

			return 0;
			}