Re: [ALPS-users] Simulation of O(N) models

Hi, from a first glance this looks good. Just type "make install" and see whether it compiles.

Matthias

On Apr 27, 2011, at 4:13 PM, Rachele Nerattini wrote:

Hi ALPS community, first of all sorry for my English and sorry if the questions that I'm going to make are trivial but I'm not a real expert in programming and using web libraries. I'm interested in doing classical MC simulations of O(N) models. I saw on the ALPS page Documentation:ClassicalMCSimulations and in particular ALPS application: spinmc that is possible to extend the application to additional O(N) models "...by editing the file mc/spins/spinmc_factory.C " I need to run a simulation for a O(4) model so I tried to edit the file spinmc_factory.C to this purpose in the following way

/*****************************************************************************

• ALPS Project Applications
• Copyright (C) 2009 by Matthias Troyer troyer@comp-phys.org,

•                   Synge Todo <wistaria@comp-phys.org> 
  

• This software is part of the ALPS Applications, published under the ALPS
• Application License; you can use, redistribute it and/or modify it under
• the terms of the license, either version 1 or (at your option) any later
• version.
• You should have received a copy of the ALPS Application License along with
• the ALPS Applications; see the file LICENSE.txt. If not, the license is also
• available from http://alps.comp-phys.org/.
• THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
• IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
• FITNESS FOR A PARTICULAR PURPOSE, TITLE AND NON-INFRINGEMENT. IN NO EVENT
• SHALL THE COPYRIGHT HOLDERS OR ANYONE DISTRIBUTING THE SOFTWARE BE LIABLE
• FOR ANY DAMAGES OR OTHER LIABILITY, WHETHER IN CONTRACT, TORT OR OTHERWISE,
• ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
• DEALINGS IN THE SOFTWARE.

*****************************************************************************/

/* $Id: spinmc_factory.C 3124 2009-02-26 09:40:13Z wistaria $ */

#include "factory.h" #include "abstractspinsim.h"

#include "spinsim.h"

#include "ising.h" #include "on.h" #include "potts.h"

void SpinFactory::print_copyright(std::ostream& out) const { out << "Generic classical Monte Carlo program using local or cluster updates\n" << " available from http://alps.comp-phys.org/%5Cn" << " copyright(c) 1999-2007 by Matthias Troyer troyer@comp-phys.org\n" << " Mathias Koerner mkoerner@comp-phys.org\n" << " for details see the publication:\n" << " A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).\n\n"; }

alps::scheduler::MCSimulation* SpinFactory::make_task( const alps::ProcessList& w, const boost::filesystem::path& fn) const { return new alps::scheduler::MCSimulation(w,fn); }

/**

• Returns the number of white-space separated elements in the string.
• \param str the string to be analyzed.

*/ int SpinFactory::countElements(const std::string& str) const { int counter = 0; int pos = 0; if (str[0] == '\0') return 0; while (true) { while (str[pos] == ' ') pos++; if (str[pos] == '\0') return counter; counter++; while ((str[pos] !='\0') && (str[pos] != ' ')) pos++; } }

/**

• Returns the number of elements in the string with the largest number of elements.

*/ int SpinFactory::findDominantMatrixString(const alps::Parameters& parms) const { int count = 0; /* should be done more elegantly ... */ if (parms.defined("J")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J"])); if (parms.defined("J'")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J'"])); if (parms.defined("J0")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J0"])); if (parms.defined("J1")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J1"])); if (parms.defined("J2")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J2"])); if (parms.defined("J3")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J3"])); if (parms.defined("D")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D"])); if (parms.defined("D'")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D'"])); if (parms.defined("D0")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D0"])); if (parms.defined("D1")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D1"])); if (parms.defined("D2")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D2"])); if (parms.defined("D3")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D3"]));

return count; }

/**

• Outputs an error message and aborts the computation if an invalid number of
• elements is provided in the input for a matrix parameter.

*/ void SpinFactory::produceError(const alps::Parameters& parms) const { std::cerr << "Invalid parameters given:\n" << "Cannot produce model " << parms["MODEL"] << " with strings " << parms["D"] << " for D and " << parms["J"] << " for J.\nAborting\n"; boost::throw_exception(std::runtime_error("invalid parameters")); std::exit(-5); }

/**

• Creates a worker, according to the parameters. The type of the produced
• worker depends on the model given in the parameter file and on the number of
• input elements for the matrix parameters. The chosen matrix type is the most
• general one (over all matrix parameters) of all most specialized ones that
• are compatible with the user input.

*/ alps::scheduler::Worker* SpinFactory::make_worker( const alps::ProcessList& where, const alps::Parameters& parms,int node) const { int maxElemCount = findDominantMatrixString(parms); if (parms["MODEL"]=="Ising") if (maxElemCount == 1) return new SpinSim<IsingMoment,MIdMatrix<double,1> >(where,parms,node); else { produceError(parms); return 0; } else if (parms["MODEL"]=="O(4)") if (maxElemCount == 1) return new SpinSim< ONMoment<4>, MIdMatrix<double,4> >(where, parms, node); else { produceError(parms); return 0; } else if (parms["MODEL"]=="Heisenberg") switch (maxElemCount) { case 1: return new SpinSim< ONMoment<3>, MIdMatrix<double,3> >(where, parms, node); break; case 3: return new SpinSim< ONMoment<3>, DiagMatrix<double,3> >(where, parms, node); break; case 6: return new SpinSim< ONMoment<3>, SquareMatrix<double,3> >(where, parms, node); break; case 9: return new SpinSim< ONMoment<3>, SquareMatrix<double,3> >(where, parms, node); break; default: produceError(parms); return 0; } else if (parms["MODEL"]=="XY") switch (maxElemCount) { case 1: return new SpinSim< XYMoment, MIdMatrix<double,2> >(where, parms, node); break; case 2: return new SpinSim< XYMoment, DiagMatrix<double,2> >(where, parms, node); break; case 3: return new SpinSim< XYMoment, SquareMatrix<double,2> >(where, parms, node); break; case 4: return new SpinSim< XYMoment, SquareMatrix<double,2> >(where,parms,node); break; default: produceError(parms); return 0; }

else if (parms["MODEL"]=="Potts") switch (int(parms["q"])) { case 3: return new SpinSim< PottsMoment<3>, MIdMatrix<double,2> >(where,parms,node); case 4: return new SpinSim< PottsMoment<4>, MIdMatrix<double,2> >(where,parms,node); case 10: return new SpinSim< PottsMoment<10>, MIdMatrix<double,2> >(where,parms,node); default: boost::throw_exception(std::runtime_error(std::string(parms["q"]) + "-state Potts model not implemented in file factory.C")); } else boost::throw_exception(std::runtime_error(std::string(parms["MODEL"]) + " model not implemented in file factory.C")); return 0; }

unfortunately I don't know this programming language so I'm not sure to have edited it in the right way.

Is it correct? And if it is so, what do I have to do after having edited the file to run the simulation?

Thank you for the help.

Best regards,

Rachele Nerattini