Skip to content
Snippets Groups Projects
Commit 556d2d9b authored by Michael Wedel's avatar Michael Wedel
Browse files

Refs #11006. Removing comments, performance improvements.

parent d03774c7
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Code/Mantid/Framework/Geometry/inc/MantidGeometry/Crystal/PointGroup.h
+ 46
− 199
View file @ 556d2d9b
......@@ -16,6 +16,46 @@
namespace Mantid {
namespace Geometry {
class MANTID_GEOMETRY_DLL IndexPermutation {
public:
IndexPermutation(const Kernel::IntMatrix &matrix) : m_factors() {
for (size_t i = 0; i < matrix.numRows(); ++i) {
m_factors.push_back(IndexSum(i, matrix[i]));
}
}
Kernel::V3D getPermutation(const Kernel::V3D &hkl) const {
Kernel::V3D hklPrime(0, 0, 0);
for (auto it = m_factors.begin(); it != m_factors.end(); ++it) {
hklPrime[(*it).getIndex()] = (*it).transformed(hkl);
}
return hklPrime;
}
private:
class IndexSum {
public:
IndexSum(size_t index, const int *matrixRow)
: m_index(index), m_f1(static_cast<double>(matrixRow[0])),
m_f2(static_cast<double>(matrixRow[1])),
m_f3(static_cast<double>(matrixRow[2])) {}
size_t getIndex() const { return m_index; }
double transformed(const Kernel::V3D &hkl) const {
return m_f1 * hkl[0] + m_f2 * hkl[1] + m_f3 * hkl[2];
}
private:
size_t m_index;
double m_f1, m_f2, m_f3;
};
std::vector<IndexSum> m_factors;
};
/** A class containing the Point Groups for a crystal.
*
* @author Vickie Lynch
......@@ -48,220 +88,27 @@ public:
CrystalSystem crystalSystem() const { return Cubic; }
/// Return true if the hkls are in same group
bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
bool isEquivalent(const Kernel::V3D &hkl, const Kernel::V3D &hkl2) const;
/// Returns a vector with all equivalent hkls
std::vector<Kernel::V3D> getEquivalents(const Kernel::V3D &hkl) const;
/// Returns the same hkl for all equivalent hkls
Kernel::V3D getReflectionFamily(const Kernel::V3D &hkl) const;
void init() { }
protected:
bool groupHasNoTranslations(const Group &group) const;
std::set<Kernel::V3D> getEquivalentSet(const Kernel::V3D &hkl) const;
std::string m_symbolHM;
std::string m_name;
};
//------------------------------------------------------------------------
/** -1 (Triclinic) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue1 : public PointGroup {
public:
PointGroupLaue1();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** 1 2/m 1 (Monoclinic, unique axis b) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue2 : public PointGroup {
public:
PointGroupLaue2();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** 1 1 2/m (Monoclinic, unique axis c) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue3 : public PointGroup {
public:
PointGroupLaue3();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** mmm (Orthorombic) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue4 : public PointGroup {
public:
PointGroupLaue4();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** 4/m (Tetragonal) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue5 : public PointGroup {
public:
PointGroupLaue5();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** 4/mmm (Tetragonal) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue6 : public PointGroup {
public:
PointGroupLaue6();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** -3 (Trigonal - Hexagonal) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue7 : public PointGroup {
public:
PointGroupLaue7();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** -3m1 (Trigonal - Rhombohedral) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue8 : public PointGroup {
public:
PointGroupLaue8();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** -31m (Trigonal - Rhombohedral) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue9 : public PointGroup {
public:
PointGroupLaue9();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** 6/m (Hexagonal) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue10 : public PointGroup {
public:
PointGroupLaue10();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
//------------------------------------------------------------------------
/** 6/mmm (Hexagonal) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue11 : public PointGroup {
public:
PointGroupLaue11();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
std::vector<Kernel::V3D> getEquivalentSet(const Kernel::V3D &hkl) const;
virtual void init();
};
void createPermutationsFromOperations();
//------------------------------------------------------------------------
/** m-3 (Cubic) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue12 : public PointGroup {
public:
PointGroupLaue12();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
std::vector<IndexPermutation> m_permutations;
virtual void init();
std::string m_symbolHM;
std::string m_name;
};
//------------------------------------------------------------------------
/** m-3m (Cubic) PointGroup
class MANTID_GEOMETRY_DLL PointGroupLaue13 : public PointGroup {
public:
PointGroupLaue13();
/// Name of the point group
virtual std::string getName() const;
/// Return true if the hkls are equivalent.
virtual bool isEquivalent(const Kernel::V3D &hkl,
const Kernel::V3D &hkl2) const;
virtual PointGroup::CrystalSystem crystalSystem() const;
virtual void init();
};
*/
/// Shared pointer to a PointGroup
typedef boost::shared_ptr<PointGroup> PointGroup_sptr;
......
Code/Mantid/Framework/Geometry/src/Crystal/PointGroup.cpp
+ 29
− 467
View file @ 556d2d9b
This diff is collapsed.
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment