Utility routines for manipulating color matrices
[Utility routines]


Files
file	reunit.cc
	Reunitarize in place a color matrix to SU(N).
file	reunit.h
	Reunitarize in place a color matrix to SU(N).
file	su3proj.cc
	Project a GL(3,C) color matrix onto SU(3).
file	su3proj.h
	Project a GL(3,C) color matrix onto SU(3).
file	sun_proj.cc
	Project a complex Nc x Nc matrix W onto SU(Nc) by maximizing Tr(VW).
file	sun_proj.h
	Project a complex Nc x Nc matrix W onto SU(Nc) by maximizing Tr(VW).
Namespaces
namespace	Chroma::GaugeInitEnv
	Registration aggregator.
namespace	Chroma::StoutLinkTimings
	Timings.
namespace	Chroma::Stouting
	Utilities.
namespace	Chroma::StoutLinkTimings
	Timings.
namespace	Chroma::StoutLinkTimings
	Timings.
namespace	Chroma::Stouting
	Utilities.
Classes
struct	Chroma::CPPACSGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::CPPACSGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
struct	Chroma::DisorderedGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::DisorderedGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
class	Chroma::GaugeInit
	Base class for gauge initialization. More...
struct	Chroma::KYUGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::KYUGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
struct	Chroma::MILCGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::MILCGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
struct	Chroma::NERSCGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::NERSCGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
struct	Chroma::SFGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::SFGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
class	Chroma::StagPhases::alphaClass
class	Chroma::StagPhases::betaClass
struct	Chroma::SZINGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::SZINGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
struct	Chroma::SZINQIOGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::SZINQIOGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
struct	Chroma::UnitGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::UnitGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
struct	Chroma::WeakGaugeInitEnv::Params
	Params for initializing config. More...
class	Chroma::WeakGaugeInitEnv::GaugeIniter
	Gauge initialization. More...
Typedefs
typedef SingletonHolder< ObjectFactory< GaugeInit, std::string, TYPELIST_2(XMLReader &, const std::string &), GaugeInit ()(XMLReader &, const std::string &), StringFactoryError > >	Chroma::TheGaugeInitFactory
	Gauge initialization factory (foundry).
Functions
void	Chroma::conjgauge (multi1d< LatticeColorMatrix > &u)
	Take complex conjugate of gauge field.
void	Chroma::constgauge (multi1d< LatticeColorMatrix > &u, const multi2d< Real > theta)
	Const diagonal gauge field.
void	Chroma::read (XMLReader &xml, const string &path, CPPACSGaugeInitEnv::Params &param)
	Reader.
void	Chroma::write (XMLWriter &xml, const string &path, const CPPACSGaugeInitEnv::Params &param)
	Writer.
void	Chroma::eesu2 (LatticeColorMatrix &m)
	Exactly exponentiate a SU(2) lie algebra element.
void	Chroma::eesu3 (LatticeColorMatrix &iQ)
	Exactly exponentiate a SU(3) lie algebra element.
void	Chroma::eeu1 (LatticeColorMatrix &lambda)
	Exactly exponentiate a U(1) lie algebra element.
void	Chroma::expm12 (LatticeColorMatrix &a)
	12-th order exponentiation of a lattice color matrix
void	Chroma::expmat (LatticeColorMatrix &a, enum ExpMat_t opt)
	Exponentiate a SU(n) lie algebra element by some method.
void	Chroma::gaugeStartup (XMLReader &gauge_file_xml, XMLReader &gauge_xml, multi1d< LatticeColorMatrix > &u, Cfg_t &cfg)
	Initialize the gauge fields.
void	Chroma::HotSt (multi1d< LatticeColorMatrix > &u)
	Set a gauge field from a sample of (almost) Haar measure.
void	Chroma::read (XMLReader &xml, const string &path, KYUGaugeInitEnv::Params &param)
	Reader.
void	Chroma::write (XMLWriter &xml, const string &path, const KYUGaugeInitEnv::Params &param)
	Writer.
void	Chroma::read (XMLReader &xml, const string &path, MILCGaugeInitEnv::Params &param)
	Reader.
void	Chroma::write (XMLWriter &xml, const string &path, const MILCGaugeInitEnv::Params &param)
	Writer.
void	Chroma::read (XMLReader &xml, const string &path, NERSCGaugeInitEnv::Params &param)
	Reader.
void	Chroma::write (XMLWriter &xml, const string &path, const NERSCGaugeInitEnv::Params &param)
	Writer.
void	Chroma::rgauge (multi1d< LatticeColorMatrix > &u, LatticeColorMatrix &g)
	Do a random gauge transformation on the u fields.
void	Chroma::rgauge (multi1d< LatticeColorMatrix > &u)
	Do a random gauge transformation on the u fields.
void	Chroma::read (XMLReader &xml, const string &path, SFGaugeInitEnv::Params &param)
	Reader.
void	Chroma::write (XMLWriter &xml, const string &path, const SFGaugeInitEnv::Params &param)
	Writer.
LatticeColorMatrix	Chroma::shift2 (const LatticeColorMatrix &s1, int isign, int dir, int level)
	A simple not-fancy power of 2 shift.
void	Chroma::Stouting::getQs (const multi1d< LatticeColorMatrix > &u, LatticeColorMatrix &Q, LatticeColorMatrix &QQ, int mu, const multi1d< bool > &smear_in_this_dirP, const multi2d< Real > &rho)
	Given field U, form Q and Q^2.
void	Chroma::Stouting::getQsandCs (const multi1d< LatticeColorMatrix > &u, LatticeColorMatrix &Q, LatticeColorMatrix &QQ, LatticeColorMatrix &C, int mu, const multi1d< bool > &smear_in_this_dirP, const multi2d< Real > &rho)
	Given field U, construct the staples into C, form Q and Q^2 and compute c0 and c1.
void	Chroma::Stouting::deriv_recurse (multi1d< LatticeColorMatrix > &F, const multi1d< bool > &smear_in_this_dirP, const multi2d< Real > &rho, const multi1d< LatticeColorMatrix > &u)
	Do the force recursion from level i+1, to level i.
void	Chroma::Stouting::getFs (const LatticeColorMatrix &Q, const LatticeColorMatrix &QQ, multi1d< LatticeComplex > &f)
	Given c0 and c1 compute the f-s and b-s.
void	Chroma::Stouting::getFsAndBs (const LatticeColorMatrix &Q, const LatticeColorMatrix &QQ, multi1d< LatticeComplex > &f, multi1d< LatticeComplex > &b1, multi1d< LatticeComplex > &b2, bool do_bs=true)
	Given c0 and c1 compute the f-s and b-s.
void	Chroma::Stouting::smear_links (const multi1d< LatticeColorMatrix > &current, multi1d< LatticeColorMatrix > &next, const multi1d< bool > &smear_in_this_dirP, const multi2d< Real > &rho)
	Do the smearing from level i to level i+1.
void	Chroma::Stouting::stout_smear (LatticeColorMatrix &next, const multi1d< LatticeColorMatrix > &current, int mu, const multi1d< bool > &smear_in_this_dirP, const multi2d< Real > &rho)
	Stout smear in a specific link direction.
template<typename S>
void	Chroma::su2Extract_t (multi1d< LatticeReal > &r, const LatticeColorMatrix &source, int su2_index, const S &s)
	Su2_extract: r_0,r_1,r_2,r_3 <- source(su2_index) [SU(N) field] under a subset.
void	Chroma::su2Extract (multi1d< LatticeReal > &r, const LatticeColorMatrix &source, int su2_index, const Subset &s)
	Su2_extract: r_0,r_1,r_2,r_3 <- source(su2_index) [SU(N) field] under a subset.
void	Chroma::sunFill (LatticeColorMatrix &dest, const multi1d< LatticeReal > &r, int su2_index, const Subset &s)
	Fill a dest(su2_index) <- r_0,r_1,r_2,r_3 under a subset.
void	Chroma::read (XMLReader &xml, const string &path, SZINGaugeInitEnv::Params &param)
	Reader.
void	Chroma::write (XMLWriter &xml, const string &path, const SZINGaugeInitEnv::Params &param)
	Writer.
void	Chroma::read (XMLReader &xml, const string &path, SZINQIOGaugeInitEnv::Params &param)
	Reader.
void	Chroma::write (XMLWriter &xml, const string &path, const SZINQIOGaugeInitEnv::Params &param)
	Writer.
void	Chroma::taproj (LatticeColorMatrix &a)
	Take the traceless antihermitian projection of a color matrix.
void	Chroma::unitarityCheck (const multi1d< LatticeColorMatrix > &u)
	Check the unitarity of color matrix in SU(N).
void	Chroma::weakField (multi1d< LatticeColorMatrix > &u)
	Construct a weak field.

Detailed Description

Utility routines for simple manipulation of gauge fields.

Typedef Documentation

typedef SingletonHolder< ObjectFactory<GaugeInit, std::string, TYPELIST_2(XMLReader&, const std::string&), GaugeInit* (*)(XMLReader&, const std::string&), StringFactoryError> > Chroma::TheGaugeInitFactory

Gauge initialization factory (foundry).

Definition at line 25 of file gauge_init_factory.h.

Function Documentation

void Chroma::conjgauge ( multi1d< LatticeColorMatrix > & u )

Take complex conjugate of gauge field.

Arguments:

Parameters:

	u	Gauge field (Modify)
	g	Gauge transformation matrices (Write)

Definition at line 21 of file conjgauge.cc.

References END_CODE, and START_CODE.

Referenced by main().

void Chroma::constgauge	(	multi1d< LatticeColorMatrix > &	u,
		const multi2d< Real >	theta
	)

Const diagonal gauge field.

Arguments:

Parameters:

	u	Gauge field (Modify)
	theta	Angles (Read)

Definition at line 24 of file constgauge.cc.

References END_CODE, and START_CODE.

Referenced by main().

void Chroma::Stouting::deriv_recurse	(	multi1d< LatticeColorMatrix > &	F,
		const multi1d< bool > &	smear_in_this_dirP,
		const multi2d< Real > &	rho,
		const multi1d< LatticeColorMatrix > &	u
	)

Do the force recursion from level i+1, to level i.

Definition at line 135 of file stout_utils.cc.

References BACKWARD, Chroma::BaryonSpinMats::C(), END_CODE, FORWARD, Chroma::Stouting::getFsAndBs(), Chroma::Stouting::getQsandCs(), and START_CODE.

Referenced by Chroma::StoutGaugeState< P, Q >::fatForceToThin(), and Chroma::StoutFermState< LatticeStaggeredPropagator, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >::fatForceToThin().

void Chroma::eesu2 ( LatticeColorMatrix & m )

Exactly exponentiate a SU(2) lie algebra element.

Parameters:

LatticeColorMatrix (Modify)

Definition at line 19 of file eesu2.cc.

References END_CODE, START_CODE, Chroma::su2Extract(), and Chroma::sunFill().

Referenced by Chroma::expmat().

void Chroma::eesu3 ( LatticeColorMatrix & m )

Exactly exponentiate a SU(3) lie algebra element.

Input: 3x3 anti-Hermitian, traceless matrix iQ Output: exp(iQ)

Formula for exp(iQ) = f0 + f1*Q + f2*Q*Q is found in section III of hep-lat/0311018.

Definition at line 21 of file eesu3.cc.

References END_CODE, START_CODE, and u.

Referenced by Chroma::expmat().

void Chroma::eeu1 ( LatticeColorMatrix & lambda )

Exactly exponentiate a U(1) lie algebra element.

Parameters:

lambda

LatticeColorMatrix (Modify)

Definition at line 19 of file eeu1.cc.

References END_CODE, and START_CODE.

Referenced by Chroma::expmat().

void Chroma::expm12 ( LatticeColorMatrix & a )

12-th order exponentiation of a lattice color matrix

In place a_ = 1 + a_ + (1/2)*a_^2 + ...+ (1/n!)*(a_)^n n = power

Arguments:

Parameters:

LatticeColorMatrix (Modify)

Definition at line 22 of file expm12.cc.

References END_CODE, and START_CODE.

Referenced by Chroma::expmat(), Chroma::HotSt(), and Chroma::rgauge().

void Chroma::expmat	(	LatticeColorMatrix &	a,
		enum ExpMat_t	opt
	)

Exponentiate a SU(n) lie algebra element by some method.

This routine is a driver for other routines. For example, expsu3 will exponentiate EFFICIENTLY an SU(3) matrix and, if desired, check for errors in the truncation of the Taylor series. The routine eesu3 EXACTLY exponentiates a SU(3) matrix.

Parameters:

	a	LatticeColorMatrix (Modify)
	opt	Method of exponentiation (Read)

Definition at line 40 of file expmat.cc.

References Chroma::eesu2(), Chroma::eesu3(), Chroma::eeu1(), END_CODE, Chroma::EXP_EXACT, Chroma::EXP_TWELFTH_ORDER, Chroma::expm12(), Chroma::expsu3(), Chroma::REUNITARIZE_LABEL, START_CODE, and Chroma::ExpMatEnv::time_spent.

Referenced by Chroma::LCMMDIntegratorSteps::leapQ(), and main().

void Chroma::gaugeStartup	(	XMLReader &	gauge_file_xml,
		XMLReader &	gauge_xml,
		multi1d< LatticeColorMatrix > &	u,
		Cfg_t &	cfg
	)

Initialize the gauge fields.

Parameters:

	gauge_file_xml	File xml
	gauge_xml	Record xml
	u	Gauge fields
	cfg	Configuration structure

Definition at line 35 of file gauge_startup.cc.

References Chroma::Cfg_t::cfg_file, Chroma::Cfg_t::cfg_type, Chroma::CFG_TYPE_CLASSICAL_SF, Chroma::CFG_TYPE_CPPACS, Chroma::CFG_TYPE_DISORDERED, Chroma::CFG_TYPE_KYU, Chroma::CFG_TYPE_MILC, Chroma::CFG_TYPE_NERSC, Chroma::CFG_TYPE_SCIDAC, Chroma::CFG_TYPE_SZIN, Chroma::CFG_TYPE_SZINQIO, Chroma::CFG_TYPE_UNIT, Chroma::CFG_TYPE_WEAK_FIELD, Chroma::SchrGaugeBCParams::decay_dir, END_CODE, Chroma::HotSt(), Chroma::SchrGaugeBCParams::loop_extent, Chroma::readCPPACS(), Chroma::readGauge(), Chroma::readKYU(), Chroma::readMILC(), Chroma::readSzin(), Chroma::reunit(), Chroma::SchrGaugeBCParams::SchrPhiMult, Chroma::SchrNonPertGaugeBC::SFBndFld(), START_CODE, Chroma::weakField(), and Chroma::write().

Referenced by main().

void Chroma::Stouting::getFs	(	const LatticeColorMatrix &	Q,
		const LatticeColorMatrix &	QQ,
		multi1d< LatticeComplex > &	f
	)

Given c0 and c1 compute the f-s and b-s.

Definition at line 332 of file stout_utils.cc.

References END_CODE, Chroma::Stouting::getFsAndBs(), and START_CODE.

Referenced by Chroma::Stouting::smear_links(), and Chroma::Stouting::stout_smear().

void Chroma::Stouting::getFsAndBs	(	const LatticeColorMatrix &	Q,
		const LatticeColorMatrix &	QQ,
		multi1d< LatticeComplex > &	f,
		multi1d< LatticeComplex > &	b1,
		multi1d< LatticeComplex > &	b2,
		bool	do_bs = `true`
	)

Given c0 and c1 compute the f-s and b-s.

Only compute b-s if do_bs is set to true (default)

Definition at line 818 of file stout_utils.cc.

References END_CODE, Chroma::StoutLinkTimings::functions_secs, Chroma::Stouting::StoutUtils::getFsAndBsSiteLoop(), and START_CODE.

Referenced by Chroma::Stouting::deriv_recurse(), Chroma::Stouting::getFs(), and main().

void Chroma::Stouting::getQs	(	const multi1d< LatticeColorMatrix > &	u,
		LatticeColorMatrix &	Q,
		LatticeColorMatrix &	QQ,
		int	mu,
		const multi1d< bool > &	smear_in_this_dirP,
		const multi2d< Real > &	rho
	)

Given field U, form Q and Q^2.

Definition at line 38 of file stout_utils.cc.

References Chroma::BaryonSpinMats::C(), END_CODE, Chroma::Stouting::getQsandCs(), and START_CODE.

Referenced by Chroma::Stouting::smear_links(), and Chroma::Stouting::stout_smear().

void Chroma::Stouting::getQsandCs	(	const multi1d< LatticeColorMatrix > &	u,
		LatticeColorMatrix &	Q,
		LatticeColorMatrix &	QQ,
		LatticeColorMatrix &	C,
		int	mu,
		const multi1d< bool > &	smear_in_this_dirP,
		const multi2d< Real > &	rho
	)

Given field U, construct the staples into C, form Q and Q^2 and compute c0 and c1.

Definition at line 54 of file stout_utils.cc.

References BACKWARD, END_CODE, FORWARD, and START_CODE.

Referenced by Chroma::Stouting::deriv_recurse(), Chroma::Stouting::getQs(), and main().

void Chroma::HotSt ( multi1d< LatticeColorMatrix > & u )

Set a gauge field from a sample of (almost) Haar measure.

u = exp(A), where A = random traceless antihermitian matrix.

Arguments:

Parameters:

Gauge field (Write)

NOTE!!!!: we are using a hack. The twelth_order exponentiation followed by a reunitarization should be replace by an exact exponentiation. So, the gauge field distriubution is NOT from Haar measure, but if someone can think of a really good reason why it should be Haar measure, then something can be done about it.

Definition at line 35 of file hotst.cc.

References END_CODE, Chroma::expm12(), Chroma::reunit(), START_CODE, and Chroma::taproj().

Referenced by Chroma::gaugeStartup(), main(), and Chroma::DisorderedGaugeInitEnv::GaugeIniter::operator()().

void Chroma::read	(	XMLReader &	xml,
		const string &	path,
		SZINQIOGaugeInitEnv::Params &	param
	)

Reader.

Definition at line 16 of file szinqio_gauge_init.cc.

void Chroma::read	(	XMLReader &	xml,
		const string &	path,
		SZINGaugeInitEnv::Params &	param
	)

Reader.

Definition at line 16 of file szin_gauge_init.cc.

void Chroma::read	(	XMLReader &	xml,
		const string &	path,
		SFGaugeInitEnv::Params &	param
	)

Reader.

Definition at line 19 of file sf_gauge_init.cc.

void Chroma::read	(	XMLReader &	xml,
		const string &	path,
		NERSCGaugeInitEnv::Params &	param
	)

Reader.

Definition at line 16 of file nersc_gauge_init.cc.

void Chroma::read	(	XMLReader &	xml,
		const string &	path,
		MILCGaugeInitEnv::Params &	param
	)

Reader.

Definition at line 16 of file milc_gauge_init.cc.

void Chroma::read	(	XMLReader &	xml,
		const string &	path,
		KYUGaugeInitEnv::Params &	param
	)

Reader.

Definition at line 16 of file kyu_gauge_init.cc.

void Chroma::read	(	XMLReader &	xml,
		const string &	path,
		CPPACSGaugeInitEnv::Params &	param
	)

Reader.

Definition at line 16 of file cppacs_gauge_init.cc.

void Chroma::rgauge ( multi1d< LatticeColorMatrix > & u )

Do a random gauge transformation on the u fields.

Convenience function: does not return gauge transformation matrices

Arguments:

Parameters:

Gauge field (Modify)

Definition at line 64 of file rgauge.cc.

References END_CODE, Chroma::rgauge(), and START_CODE.

void Chroma::rgauge	(	multi1d< LatticeColorMatrix > &	u,
		LatticeColorMatrix &	g
	)

Do a random gauge transformation on the u fields.

Arguments:

Parameters:

	u	Gauge field (Modify)
	g	Gauge transformation matrices (Write)

Definition at line 25 of file rgauge.cc.

References END_CODE, Chroma::expm12(), FORWARD, Chroma::reunit(), START_CODE, and Chroma::taproj().

Referenced by main(), Chroma::InlineRandomTransfGaugeEnv::InlineMeas::operator()(), and Chroma::rgauge().

LatticeColorMatrix Chroma::shift2	(	const LatticeColorMatrix &	s1,
		int	isign,
		int	dir,
		int	level
	)

A simple not-fancy power of 2 shift.

Definition at line 14 of file shift2.cc.

References END_CODE, and START_CODE.

Referenced by Chroma::APE_Smear(), Chroma::block(), Chroma::gluecor(), and Chroma::polycor().

void Chroma::Stouting::smear_links	(	const multi1d< LatticeColorMatrix > &	current,
		multi1d< LatticeColorMatrix > &	next,
		const multi1d< bool > &	smear_in_this_dirP,
		const multi2d< Real > &	rho
	)

Do the smearing from level i to level i+1.

Definition at line 848 of file stout_utils.cc.

References END_CODE, Chroma::Stouting::getFs(), Chroma::Stouting::getQs(), and START_CODE.

Referenced by Chroma::StoutGaugeState< P, Q >::create(), Chroma::StoutFermState< LatticeStaggeredPropagator, multi1d< LatticeColorMatrix >, multi1d< LatticeColorMatrix > >::create(), and Chroma::StoutLinkSmearingEnv::LinkSmear::operator()().

void Chroma::Stouting::stout_smear	(	LatticeColorMatrix &	next,
		const multi1d< LatticeColorMatrix > &	current,
		int	mu,
		const multi1d< bool > &	smear_in_this_dirP,
		const multi2d< Real > &	rho
	)

Stout smear in a specific link direction.

Definition at line 882 of file stout_utils.cc.

References END_CODE, Chroma::Stouting::getFs(), Chroma::Stouting::getQs(), and START_CODE.

Referenced by main().

void Chroma::su2Extract	(	multi1d< LatticeReal > &	r,
		const LatticeColorMatrix &	source,
		int	su2_index,
		const Subset &	s
	)

Su2_extract: r_0,r_1,r_2,r_3 <- source(su2_index) [SU(N) field] under a subset.

Extract components r_k proportional to SU(2) submatrix su2_index from the "SU(Nc)" matrix V. The SU(2) matrix is parametrized in the sigma matrix basis.

There are Nc*(Nc-1)/2 unique SU(2) submatrices in an SU(Nc) matrix. The user does not need to know exactly which one is which, just that they are unique.

Arguments:

Parameters:

	r	su(2) matrix represented in the O(4) rep. - an array of LatticeReal
	source	extract the su2 matrix from this su(n) gauge field
	su2_index	int lying in [0, Nc*(Nc-1)/2)
	s	subset for operations (Read)

Definition at line 87 of file su2extract.cc.

References Chroma::su2Extract_t().

Referenced by Chroma::eesu2(), Chroma::grelax(), Chroma::su2_hb_update(), Chroma::su3hb(), Chroma::su3over(), and Chroma::su3proj_t().

template<typename S>

void Chroma::su2Extract_t	(	multi1d< LatticeReal > &	r,
		const LatticeColorMatrix &	source,
		int	su2_index,
		const S &	s
	)			`[inline]`

Su2_extract: r_0,r_1,r_2,r_3 <- source(su2_index) [SU(N) field] under a subset.

Extract components r_k proportional to SU(2) submatrix su2_index from the "SU(Nc)" matrix V. The SU(2) matrix is parametrized in the sigma matrix basis.

There are Nc*(Nc-1)/2 unique SU(2) submatrices in an SU(Nc) matrix. The user does not need to know exactly which one is which, just that they are unique.

Arguments:

Parameters:

	r	su(2) matrix represented in the O(4) rep. - an array of LatticeReal
	source	extract the su2 matrix from this su(n) gauge field
	su2_index	int lying in [0, Nc*(Nc-1)/2)
	s	subset for operations (Read)

Definition at line 35 of file su2extract.cc.

References END_CODE, and START_CODE.

Referenced by Chroma::su2Extract().

void Chroma::sunFill	(	LatticeColorMatrix &	dest,
		const multi1d< LatticeReal > &	r,
		int	su2_index,
		const Subset &	s
	)

Fill a dest(su2_index) <- r_0,r_1,r_2,r_3 under a subset.

Fill an SU(Nc) matrix V with the SU(2) submatrix su2_index paramtrized by b_k in the sigma matrix basis.

Fill in B from B_SU(2) = b0 + i sum_k bk sigma_k

There are Nc*(Nc-1)/2 unique SU(2) submatrices in an SU(Nc) matrix. The user does not need to know exactly which one is which, just that they are unique.

Arguments:

Parameters:

	dest	su(n) matrix
	r	su2 matrix represented in the O(4) rep. - an array of LatticeReal
	su2_index	int lying in [0, Nc*(Nc-1)/2)
	s	subset for operations (Read)