alpsize-src-20130306-r2942/0000775000127400003120000000000012117341400014116 5ustar wistarialabalpsize-src-20130306-r2942/07_lattice/0000775000127400003120000000000012117341400016051 5ustar wistarialabalpsize-src-20130306-r2942/07_lattice/lattice.C0000664000127400003120000001027611364020401017606 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include int main() { alps::Parameters params(std::cin); alps::graph_helper<> graph(params); std::cout << "lattice name = " << params["LATTICE"] << std::endl; // // number of sites // std::cout << "number of sites = " << graph.num_sites() << std::endl; // // number of bonds // std::cout << "number of bonds = " << graph.num_bonds() << std::endl; // // iterate over all the sites and print out indices // // (1) for (int i = 0; i < graph.num_sites(); ++i) std::cout << "site index = " << i << std::endl; // (2) alps::graph_helper<>::site_iterator itr0, itr0_end; for (boost::tie(itr0, itr0_end) = graph.sites(); itr0 != itr0_end; ++itr0) { std::cout << "site index = " << graph.index(*itr0) << std::endl; // this is equivalent to alps::graph_helper<>::site_descriptor desc = *itr0; std::cout << "site index = " << graph.index(desc) << std::endl; // or simply std::cout << "site index = " << desc << std::endl; // or just std::cout << "site index = " << *itr0 << std::endl; } // (3) (using BOOST_FOREACH) BOOST_FOREACH(const alps::graph_helper<>::site_descriptor& s, graph.sites()) { std::cout << "site index = " << graph.index(s) << std::endl; // or just std::cout << "site index = " << s << std::endl; } // (4) since site_descriptor is equivalent to int the following code also works BOOST_FOREACH(int s, graph.sites()) std::cout << "site index = " << s << std::endl; // // iterate over all the bonds and print out indices // // (1) for (int i = 0; i < graph.num_bonds(); ++i) std::cout << "bond index = " << i << std::endl; // (2) alps::graph_helper<>::bond_iterator itr1, itr1_end; for (boost::tie(itr1, itr1_end) = graph.bonds(); itr1 != itr1_end; ++itr1) { std::cout << "bond index = " << graph.index(*itr1) << std::endl; // this is equivalent to alps::graph_helper<>::bond_descriptor desc = *itr1; std::cout << "bond index = " << graph.index(desc) << std::endl; // Warning: the following codes do not output the bond index, but pair of its source and target std::cout << "bond index != " << desc << std::endl; std::cout << "bond index != " << *itr1 << std::endl; } // (3) (using BOOST_FOREACH) BOOST_FOREACH(const alps::graph_helper<>::bond_descriptor& b, graph.bonds()) { std::cout << "bond index = " << graph.index(b) << std::endl; // Warning: the following code does not output the bond index, but pair of its source and target std::cout << "bond index != " << b << std::endl; } return 0; } alpsize-src-20130306-r2942/07_lattice/wolff.C0000664000127400003120000000772011364020401017276 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include #include #include #include #include #include #include int main() { alps::Parameters params(std::cin); const double T = params.value_or_default("T", 2.2); const int MCSTEP = params.value_or_default("SWEEPS", 1 << 15); const int MCTHRM = params.value_or_default("THERMALIZATION", MCSTEP >> 3); const unsigned int SEED = params.value_or_default("SEED", 93812); // setting up square lattice alps::graph_helper<> graph(params); const int N = graph.num_sites(); // random number generator boost::mt19937 eng(SEED); boost::variate_generator > random_01(eng, boost::uniform_real<>()); // spin configuration std::vector spin(N, 1); int sz = N; // stack for uninspected sites std::stack stck; // connecting probability double pc = 1 - std::exp(-2./T); // measurement alps::ObservableSet measurements; measurements << alps::RealObservable("Magnetization"); measurements << alps::RealObservable("Magnetization^2"); measurements << alps::RealObservable("Magnetization^4"); // timer boost::timer tm; for (int mcs = 0; mcs < MCSTEP + MCTHRM; ++mcs) { if (mcs == MCTHRM) measurements.reset(true); int s = static_cast(random_01() * N); int so = spin[s]; spin[s] = -so; stck.push(s); int cs = 0; while (!stck.empty()) { ++cs; int sc = stck.top(); stck.pop(); BOOST_FOREACH(alps::graph_helper<>::site_descriptor const& sn, graph.neighbors(sc)) { if (spin[sn] == so && random_01() < pc) { stck.push(sn); spin[sn] = -so; } } } sz -= 2 * so * cs; double dsz = sz / static_cast(N); measurements["Magnetization"] << dsz; measurements["Magnetization^2"] << dsz * dsz; measurements["Magnetization^4"] << dsz * dsz * dsz * dsz; } // calculate Binder parameter alps::RealObsevaluator m2 = measurements["Magnetization^2"]; alps::RealObsevaluator m4 = measurements["Magnetization^4"]; alps::RealObsevaluator binder("Binder Ratio of Magnetization"); binder = m2 * m2 / m4; measurements.addObservable(binder); // output results std::cout << measurements; std::cerr << "Elapsed time = " << tm.elapsed() << " sec\n"; return 0; } alpsize-src-20130306-r2942/07_lattice/wolff.op0000664000127400003120000000311211071061125017524 0ustar wistarialabBinder Ratio of Magnetization: 0.958369 +/- 0.0011 Magnetization: -0.00363082 +/- 0.00231; tau = -0.36 bin #1 : 32768 entries: error = 0.00436488 bin #2 : 16384 entries: error = 0.00276793 bin #3 : 8192 entries: error = 0.00261629 bin #4 : 4096 entries: error = 0.00247348 bin #5 : 2048 entries: error = 0.00234902 bin #6 : 1024 entries: error = 0.00230857 bin #7 : 512 entries: error = 0.00221888 bin #8 : 256 entries: error = 0.00226437 bin #9 : 128 entries: error = 0.00231019 Magnetization^2: 0.624276 +/- 0.00147; tau = 1.59 bin #1 : 32768 entries: error = 0.000718779 bin #2 : 16384 entries: error = 0.000894127 bin #3 : 8192 entries: error = 0.00107617 bin #4 : 4096 entries: error = 0.00122415 bin #5 : 2048 entries: error = 0.00131216 bin #6 : 1024 entries: error = 0.00140491 bin #7 : 512 entries: error = 0.00141392 bin #8 : 256 entries: error = 0.00144545 bin #9 : 128 entries: error = 0.00146852 Magnetization^4: 0.406649 +/- 0.00153; tau = 1.43 bin #1 : 32768 entries: error = 0.000777249 bin #2 : 16384 entries: error = 0.000963571 bin #3 : 8192 entries: error = 0.00115421 bin #4 : 4096 entries: error = 0.0013069 bin #5 : 2048 entries: error = 0.00138956 bin #6 : 1024 entries: error = 0.00148461 bin #7 : 512 entries: error = 0.00150367 bin #8 : 256 entries: error = 0.00152819 bin #9 : 128 entries: error = 0.00152892 alpsize-src-20130306-r2942/07_lattice/wolff.ip0000664000127400003120000000005211071061125017516 0ustar wistarialabLATTICE = "square lattice" L = 32 T = 2.2 alpsize-src-20130306-r2942/07_lattice/lattice.op0000664000127400003120000001363011071061125020042 0ustar wistarialablattice name = square lattice number of sites = 16 number of bonds = 32 site index = 0 site index = 1 site index = 2 site index = 3 site index = 4 site index = 5 site index = 6 site index = 7 site index = 8 site index = 9 site index = 10 site index = 11 site index = 12 site index = 13 site index = 14 site index = 15 site index = 0 site index = 0 site index = 0 site index = 0 site index = 1 site index = 1 site index = 1 site index = 1 site index = 2 site index = 2 site index = 2 site index = 2 site index = 3 site index = 3 site index = 3 site index = 3 site index = 4 site index = 4 site index = 4 site index = 4 site index = 5 site index = 5 site index = 5 site index = 5 site index = 6 site index = 6 site index = 6 site index = 6 site index = 7 site index = 7 site index = 7 site index = 7 site index = 8 site index = 8 site index = 8 site index = 8 site index = 9 site index = 9 site index = 9 site index = 9 site index = 10 site index = 10 site index = 10 site index = 10 site index = 11 site index = 11 site index = 11 site index = 11 site index = 12 site index = 12 site index = 12 site index = 12 site index = 13 site index = 13 site index = 13 site index = 13 site index = 14 site index = 14 site index = 14 site index = 14 site index = 15 site index = 15 site index = 15 site index = 15 site index = 0 site index = 0 site index = 1 site index = 1 site index = 2 site index = 2 site index = 3 site index = 3 site index = 4 site index = 4 site index = 5 site index = 5 site index = 6 site index = 6 site index = 7 site index = 7 site index = 8 site index = 8 site index = 9 site index = 9 site index = 10 site index = 10 site index = 11 site index = 11 site index = 12 site index = 12 site index = 13 site index = 13 site index = 14 site index = 14 site index = 15 site index = 15 site index = 0 site index = 1 site index = 2 site index = 3 site index = 4 site index = 5 site index = 6 site index = 7 site index = 8 site index = 9 site index = 10 site index = 11 site index = 12 site index = 13 site index = 14 site index = 15 bond index = 0 bond index = 1 bond index = 2 bond index = 3 bond index = 4 bond index = 5 bond index = 6 bond index = 7 bond index = 8 bond index = 9 bond index = 10 bond index = 11 bond index = 12 bond index = 13 bond index = 14 bond index = 15 bond index = 16 bond index = 17 bond index = 18 bond index = 19 bond index = 20 bond index = 21 bond index = 22 bond index = 23 bond index = 24 bond index = 25 bond index = 26 bond index = 27 bond index = 28 bond index = 29 bond index = 30 bond index = 31 bond index = 0 bond index = 0 bond index != (0,1) bond index != (0,1) bond index = 1 bond index = 1 bond index != (0,4) bond index != (0,4) bond index = 2 bond index = 2 bond index != (1,2) bond index != (1,2) bond index = 3 bond index = 3 bond index != (1,5) bond index != (1,5) bond index = 4 bond index = 4 bond index != (2,3) bond index != (2,3) bond index = 5 bond index = 5 bond index != (2,6) bond index != (2,6) bond index = 6 bond index = 6 bond index != (3,0) bond index != (3,0) bond index = 7 bond index = 7 bond index != (3,7) bond index != (3,7) bond index = 8 bond index = 8 bond index != (4,5) bond index != (4,5) bond index = 9 bond index = 9 bond index != (4,8) bond index != (4,8) bond index = 10 bond index = 10 bond index != (5,6) bond index != (5,6) bond index = 11 bond index = 11 bond index != (5,9) bond index != (5,9) bond index = 12 bond index = 12 bond index != (6,7) bond index != (6,7) bond index = 13 bond index = 13 bond index != (6,10) bond index != (6,10) bond index = 14 bond index = 14 bond index != (7,4) bond index != (7,4) bond index = 15 bond index = 15 bond index != (7,11) bond index != (7,11) bond index = 16 bond index = 16 bond index != (8,9) bond index != (8,9) bond index = 17 bond index = 17 bond index != (8,12) bond index != (8,12) bond index = 18 bond index = 18 bond index != (9,10) bond index != (9,10) bond index = 19 bond index = 19 bond index != (9,13) bond index != (9,13) bond index = 20 bond index = 20 bond index != (10,11) bond index != (10,11) bond index = 21 bond index = 21 bond index != (10,14) bond index != (10,14) bond index = 22 bond index = 22 bond index != (11,8) bond index != (11,8) bond index = 23 bond index = 23 bond index != (11,15) bond index != (11,15) bond index = 24 bond index = 24 bond index != (12,13) bond index != (12,13) bond index = 25 bond index = 25 bond index != (12,0) bond index != (12,0) bond index = 26 bond index = 26 bond index != (13,14) bond index != (13,14) bond index = 27 bond index = 27 bond index != (13,1) bond index != (13,1) bond index = 28 bond index = 28 bond index != (14,15) bond index != (14,15) bond index = 29 bond index = 29 bond index != (14,2) bond index != (14,2) bond index = 30 bond index = 30 bond index != (15,12) bond index != (15,12) bond index = 31 bond index = 31 bond index != (15,3) bond index != (15,3) bond index = 0 bond index != (0,1) bond index = 1 bond index != (0,4) bond index = 2 bond index != (1,2) bond index = 3 bond index != (1,5) bond index = 4 bond index != (2,3) bond index = 5 bond index != (2,6) bond index = 6 bond index != (3,0) bond index = 7 bond index != (3,7) bond index = 8 bond index != (4,5) bond index = 9 bond index != (4,8) bond index = 10 bond index != (5,6) bond index = 11 bond index != (5,9) bond index = 12 bond index != (6,7) bond index = 13 bond index != (6,10) bond index = 14 bond index != (7,4) bond index = 15 bond index != (7,11) bond index = 16 bond index != (8,9) bond index = 17 bond index != (8,12) bond index = 18 bond index != (9,10) bond index = 19 bond index != (9,13) bond index = 20 bond index != (10,11) bond index = 21 bond index != (10,14) bond index = 22 bond index != (11,8) bond index = 23 bond index != (11,15) bond index = 24 bond index != (12,13) bond index = 25 bond index != (12,0) bond index = 26 bond index != (13,14) bond index = 27 bond index != (13,1) bond index = 28 bond index != (14,15) bond index = 29 bond index != (14,2) bond index = 30 bond index != (15,12) bond index = 31 bond index != (15,3) alpsize-src-20130306-r2942/07_lattice/CMakeLists.txt0000664000127400003120000000105612115500225020613 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(lattice NONE) # find ALPS Library find_package(ALPS REQUIRED PATHS ${ALPS_ROOT_DIR} $ENV{ALPS_HOME} NO_SYSTEM_ENVIRONMENT_PATH) message(STATUS "Found ALPS: ${ALPS_ROOT_DIR} (revision: ${ALPS_VERSION})") include(${ALPS_USE_FILE}) # enable C++ compilers enable_language(CXX) # lattice add_executable(lattice lattice.C) target_link_libraries(lattice ${ALPS_LIBRARIES}) add_alps_test(lattice) # wolff add_executable(wolff wolff.C) target_link_libraries(wolff ${ALPS_LIBRARIES}) add_alps_test(wolff) alpsize-src-20130306-r2942/07_lattice/lattice.ip0000664000127400003120000000004111071061125020024 0ustar wistarialabLATTICE = "square lattice" L = 4 alpsize-src-20130306-r2942/LICENSE0000664000127400003120000000247211071061125015131 0ustar wistarialabBoost Software License - Version 1.0 - August 17th, 2003 Permission is hereby granted, free of charge, to any person or organization obtaining a copy of the software and accompanying documentation covered by this license (the "Software") to use, reproduce, display, distribute, execute, and transmit the Software, and to prepare derivative works of the Software, and to permit third-parties to whom the Software is furnished to do so, all subject to the following: The copyright notices in the Software and this entire statement, including the above license grant, this restriction and the following disclaimer, must be included in all copies of the Software, in whole or in part, and all derivative works of the Software, unless such copies or derivative works are solely in the form of machine-executable object code generated by a source language processor. 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. alpsize-src-20130306-r2942/ChangeLog0000664000127400003120000000147711401642216015705 0ustar wistarialab2010-06-03 Synge Todo - Version 2.0b1 for ALPS 2.0 tutorial 2010-04-22 Synge Todo - Updated for configuring with cmake 2008-11-09 Synge Todo - ALPS/parapack has been merged into the ALPS Library 2008-10-14 Synge Todo - Updated for new alps/parapack worker interface 2008-08-21 Synge Todo - Added plot XML sample: scheduler/plot-binder - Fixed lattice/lattice.C 2008-08-19 Synge Todo * version 1.0 released for ALPS Programming Workshop - Added measurement of Magnetization^4 and Binder ratio - Added scheduler/hello program - Added missing files to distribution tarball 2008-08-19 Synge Todo * version 1.0b1 released alpsize-src-20130306-r2942/08_scheduler/0000775000127400003120000000000012117341400016403 5ustar wistarialabalpsize-src-20130306-r2942/08_scheduler/hello.C0000664000127400003120000000326611364020401017617 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include int main(int argc, char** argv) { return alps::parapack::start(argc, argv); } alpsize-src-20130306-r2942/08_scheduler/wolff.C0000664000127400003120000000326611364020401017631 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include int main(int argc, char** argv) { return alps::parapack::start(argc, argv); } alpsize-src-20130306-r2942/08_scheduler/wolff.op0000664000127400003120000001262611413306713020075 0ustar wistarialab[input parameters] ALGORITHM = wolff; SEED = 93812; LATTICE = "square lattice"; L = 32; T = 2.2; DIR_NAME = .; BASE_NAME = task1; WORKER_SEED = 93812; DISORDER_SEED = 93812; [results] Binder Ratio of Magnetization: 0.95905 +/- 0.001 Magnetization: 0.00188583 +/- 0.002; tau = -0.396 bin #1 : 32768 entries: error = 0.00437008 bin #2 : 16384 entries: error = 0.00275661 bin #3 : 8192 entries: error = 0.00259806 bin #4 : 4096 entries: error = 0.00240655 bin #5 : 2048 entries: error = 0.00229005 bin #6 : 1024 entries: error = 0.00215442 bin #7 : 512 entries: error = 0.00214279 bin #8 : 256 entries: error = 0.00219167 bin #9 : 128 entries: error = 0.00199722 Magnetization^2: 0.625756 +/- 0.00152; tau = 1.75 bin #1 : 32768 entries: error = 0.000714308 bin #2 : 16384 entries: error = 0.000887705 bin #3 : 8192 entries: error = 0.00107367 bin #4 : 4096 entries: error = 0.00122594 bin #5 : 2048 entries: error = 0.00133008 bin #6 : 1024 entries: error = 0.00139734 bin #7 : 512 entries: error = 0.00146874 bin #8 : 256 entries: error = 0.00155281 bin #9 : 128 entries: error = 0.00151519 Magnetization^4: 0.408288 +/- 0.00164; tau = 1.73 bin #1 : 32768 entries: error = 0.000775868 bin #2 : 16384 entries: error = 0.000960415 bin #3 : 8192 entries: error = 0.00115747 bin #4 : 4096 entries: error = 0.00131761 bin #5 : 2048 entries: error = 0.00141932 bin #6 : 1024 entries: error = 0.00148684 bin #7 : 512 entries: error = 0.00155568 bin #8 : 256 entries: error = 0.00166161 bin #9 : 128 entries: error = 0.0016398 [input parameters] ALGORITHM = wolff; SEED = 93812; LATTICE = "square lattice"; L = 32; T = 2.4; DIR_NAME = .; BASE_NAME = task2; WORKER_SEED = 93812; DISORDER_SEED = 93812; [results] Binder Ratio of Magnetization: 0.496104 +/- 0.00425 WARNING: ERRORS NOT CONVERGED!!! Magnetization: -0.00432694 +/- 0.00341; tau = 0.898 WARNING: check error convergence bin #1 : 32768 entries: error = 0.00203631 bin #2 : 16384 entries: error = 0.00228139 bin #3 : 8192 entries: error = 0.00252055 bin #4 : 4096 entries: error = 0.00267857 bin #5 : 2048 entries: error = 0.00280677 bin #6 : 1024 entries: error = 0.00295545 bin #7 : 512 entries: error = 0.00302941 bin #8 : 256 entries: error = 0.003109 bin #9 : 128 entries: error = 0.00340507 Magnetization^2: 0.135885 +/- 0.00221; tau = 3.75 WARNING: ERRORS NOT CONVERGED!!! bin #1 : 32768 entries: error = 0.000756465 bin #2 : 16384 entries: error = 0.00096098 bin #3 : 8192 entries: error = 0.00119274 bin #4 : 4096 entries: error = 0.00143126 bin #5 : 2048 entries: error = 0.00164607 bin #6 : 1024 entries: error = 0.00180887 bin #7 : 512 entries: error = 0.00195939 bin #8 : 256 entries: error = 0.0020358 bin #9 : 128 entries: error = 0.00220578 Magnetization^4: 0.0372146 +/- 0.000986; tau = 3.62 WARNING: ERRORS NOT CONVERGED!!! bin #1 : 32768 entries: error = 0.000343404 bin #2 : 16384 entries: error = 0.000433857 bin #3 : 8192 entries: error = 0.000536235 bin #4 : 4096 entries: error = 0.00064157 bin #5 : 2048 entries: error = 0.000730072 bin #6 : 1024 entries: error = 0.000800745 bin #7 : 512 entries: error = 0.00086277 bin #8 : 256 entries: error = 0.000890464 bin #9 : 128 entries: error = 0.000986325 [input parameters] ALGORITHM = wolff; SEED = 93812; LATTICE = "square lattice"; L = 32; T = 2.6; DIR_NAME = .; BASE_NAME = task3; WORKER_SEED = 93812; DISORDER_SEED = 93812; [results] Binder Ratio of Magnetization: 0.362371 +/- 0.00604 Magnetization: -0.00364673 +/- 0.00301; tau = 3.71 bin #1 : 32768 entries: error = 0.00103616 bin #2 : 16384 entries: error = 0.00138227 bin #3 : 8192 entries: error = 0.00178963 bin #4 : 4096 entries: error = 0.00221601 bin #5 : 2048 entries: error = 0.0025987 bin #6 : 1024 entries: error = 0.00282834 bin #7 : 512 entries: error = 0.0029904 bin #8 : 256 entries: error = 0.00297498 bin #9 : 128 entries: error = 0.0030068 Magnetization^2: 0.0351918 +/- 0.000662; tau = 2.8 bin #1 : 32768 entries: error = 0.000257812 bin #2 : 16384 entries: error = 0.000336445 bin #3 : 8192 entries: error = 0.000422699 bin #4 : 4096 entries: error = 0.000505001 bin #5 : 2048 entries: error = 0.00057723 bin #6 : 1024 entries: error = 0.00062865 bin #7 : 512 entries: error = 0.00067105 bin #8 : 256 entries: error = 0.00068286 bin #9 : 128 entries: error = 0.000662371 Magnetization^4: 0.00341633 +/- 0.000119; tau = 2.09 bin #1 : 32768 entries: error = 5.23258e-05 bin #2 : 16384 entries: error = 6.76359e-05 bin #3 : 8192 entries: error = 8.28588e-05 bin #4 : 4096 entries: error = 9.6803e-05 bin #5 : 2048 entries: error = 0.000108733 bin #6 : 1024 entries: error = 0.000118842 bin #7 : 512 entries: error = 0.000124315 bin #8 : 256 entries: error = 0.000126642 bin #9 : 128 entries: error = 0.000119006 alpsize-src-20130306-r2942/08_scheduler/hello.ip0000664000127400003120000000020411413306713020042 0ustar wistarialabALGORITHM = "hello" SEED = 0 { WORLD = "world" } { WORLD = "alps" } { WORLD = "looper" } { WORLD = "japan" } { WORLD = "wistaria" } alpsize-src-20130306-r2942/08_scheduler/hello.op0000664000127400003120000000141511413306713020055 0ustar wistarialab[input parameters] ALGORITHM = hello; SEED = 0; WORLD = world; DIR_NAME = .; BASE_NAME = task1; WORKER_SEED = 0; DISORDER_SEED = 0; hello, world [results] [input parameters] ALGORITHM = hello; SEED = 0; WORLD = alps; DIR_NAME = .; BASE_NAME = task2; WORKER_SEED = 0; DISORDER_SEED = 0; hello, alps [results] [input parameters] ALGORITHM = hello; SEED = 0; WORLD = looper; DIR_NAME = .; BASE_NAME = task3; WORKER_SEED = 0; DISORDER_SEED = 0; hello, looper [results] [input parameters] ALGORITHM = hello; SEED = 0; WORLD = japan; DIR_NAME = .; BASE_NAME = task4; WORKER_SEED = 0; DISORDER_SEED = 0; hello, japan [results] [input parameters] ALGORITHM = hello; SEED = 0; WORLD = wistaria; DIR_NAME = .; BASE_NAME = task5; WORKER_SEED = 0; DISORDER_SEED = 0; hello, wistaria [results] alpsize-src-20130306-r2942/08_scheduler/hello_params0000664000127400003120000000017311364017211021000 0ustar wistarialabALGORITHM = "hello" { WORLD = "world" } { WORLD = "alps" } { WORLD = "looper" } { WORLD = "japan" } { WORLD = "wistaria" } alpsize-src-20130306-r2942/08_scheduler/wolff_worker.h0000664000127400003120000000775711364020401021300 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include class wolff_worker : public alps::parapack::lattice_mc_worker<> { private: typedef alps::parapack::lattice_mc_worker<> super_type; public: wolff_worker(alps::Parameters const& params) : super_type(params) { T = params.value_or_default("T", 2.2); mcs = 0; MCSTEP = params.value_or_default("SWEEPS", 1 << 15); MCTHRM = params.value_or_default("THERMALIZATION", MCSTEP >> 3); N = num_sites(); // number of lattice sites spin.resize(N, 1); // spin configuration sz = N; pc = 1 - std::exp(-2./T); // connecting probability } virtual ~wolff_worker() {} void init_observables(alps::Parameters const&, alps::ObservableSet& obs) { obs << alps::RealObservable("Magnetization") << alps::RealObservable("Magnetization^2") << alps::RealObservable("Magnetization^4"); } bool is_thermalized() const { return mcs >= MCTHRM; } double progress() const { return 1.0 * mcs / (MCTHRM + MCSTEP); } void run(alps::ObservableSet& obs) { ++mcs; int s = static_cast(random_01() * N); int so = spin[s]; spin[s] = -so; stck.push(s); int cs = 0; while (!stck.empty()) { ++cs; int sc = stck.top(); stck.pop(); BOOST_FOREACH(alps::graph_helper<>::site_descriptor const& sn, neighbors(sc)) { if (spin[sn] == so && random_01() < pc) { stck.push(sn); spin[sn] = -so; } } } sz -= 2 * so * cs; double dsz = sz / static_cast(N); obs["Magnetization"] << dsz; obs["Magnetization^2"] << dsz * dsz; obs["Magnetization^4"] << dsz * dsz * dsz * dsz; } void save(alps::ODump& dp) const { dp << mcs << spin << sz; } void load(alps::IDump& dp) { dp >> mcs >> spin >> sz; } private: double T; int mcs; int MCSTEP; int MCTHRM; int N; // number of lattice sites std::vector spin; // spin configuration int sz; std::stack stck; // stack for uninspected sites double pc; // connecting probability }; class wolff_evaluator : public alps::parapack::simple_evaluator { public: wolff_evaluator(alps::Parameters const&) {} void evaluate(alps::ObservableSet& obs) const { alps::RealObsevaluator m2 = obs["Magnetization^2"]; alps::RealObsevaluator m4 = obs["Magnetization^4"]; alps::RealObsevaluator binder("Binder Ratio of Magnetization"); binder = m2 * m2 / m4; obs.addObservable(binder); } }; alpsize-src-20130306-r2942/08_scheduler/wolff.ip0000664000127400003120000000014711413306713020062 0ustar wistarialabALGORITHM = "wolff" SEED = 93812 LATTICE = "square lattice" L = 32 { T = 2.2 } { T = 2.4 } { T = 2.6 } alpsize-src-20130306-r2942/08_scheduler/wolff_worker.C0000664000127400003120000000365411364020401021223 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include "wolff_worker.h" PARAPACK_SET_VERSION("Tutorial: How to ALPSize your applications"); PARAPACK_SET_COPYRIGHT("Tutorial: How to ALPSize your applications\n" " Copyright (C) 2005-2010 by Synge Todo "); PARAPACK_REGISTER_WORKER(wolff_worker, "wolff"); PARAPACK_REGISTER_EVALUATOR(wolff_evaluator, "wolff"); alpsize-src-20130306-r2942/08_scheduler/wolff_params0000664000127400003120000000015211364017211021007 0ustar wistarialabALGORITHM = "wolff" LATTICE = "square lattice" L = 32 NUM_CLONES = 16 { T = 2.2 } { T = 2.4 } { T = 2.6 } alpsize-src-20130306-r2942/08_scheduler/hello_worker.h0000664000127400003120000000417711364020401021257 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include class hello_worker : public alps::parapack::abstract_worker { public: hello_worker(alps::Parameters const& params) { std::cout << "hello, " << params["WORLD"] << std::endl; } virtual ~hello_worker() {} void init_observables(alps::Parameters const&, alps::ObservableSet&) { // initialize observable set } bool is_thermalized() const { return true; } double progress() const { return 1.0; } void run(alps::ObservableSet&) { // do some work } void save(alps::ODump&) const {} void load(alps::IDump&) {} }; alpsize-src-20130306-r2942/08_scheduler/hello_worker.C0000664000127400003120000000344311364020401021205 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include "hello_worker.h" PARAPACK_SET_VERSION("my version"); PARAPACK_SET_COPYRIGHT("my copyright"); PARAPACK_REGISTER_WORKER(hello_worker, "hello"); PARAPACK_REGISTER_EVALUATOR(alps::parapack::simple_evaluator, "hello"); alpsize-src-20130306-r2942/08_scheduler/CMakeLists.txt0000664000127400003120000000110312017313126021141 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(scheduler NONE) # find ALPS Library find_package(ALPS REQUIRED PATHS ${ALPS_ROOT_DIR} $ENV{ALPS_HOME} NO_SYSTEM_ENVIRONMENT_PATH) message(STATUS "Found ALPS: ${ALPS_ROOT_DIR} (revision: ${ALPS_VERSION})") include(${ALPS_USE_FILE}) # enable C++ compilers enable_language(CXX) # hello add_executable(hello hello.C hello_worker.C) target_link_libraries(hello ${ALPS_LIBRARIES}) add_alps_test(hello) # wolff add_executable(wolff wolff.C wolff_worker.C) target_link_libraries(wolff ${ALPS_LIBRARIES}) add_alps_test(wolff) alpsize-src-20130306-r2942/06_alea/0000775000127400003120000000000012117341400015325 5ustar wistarialabalpsize-src-20130306-r2942/06_alea/wolff.C0000664000127400003120000001024611364017211016554 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include #include #include #include #include #include int main() { alps::Parameters params(std::cin); const int L = params.value_or_default("L", 32); const int N = L * L; const double T = params.value_or_default("T", 2.2); const int MCSTEP = params.value_or_default("SWEEPS", 1 << 15); const int MCTHRM = params.value_or_default("THERMALIZATION", MCSTEP >> 3); const unsigned int SEED = params.value_or_default("SEED", 93812); // setting up square lattice std::vector > nn(N); for (int y = 0; y < L; ++y) for (int x = 0; x < L; ++x) { nn[x+L*y][0] = ((x+L-1)%L) + L*y; nn[x+L*y][1] = ((x+1)%L) + L*y; nn[x+L*y][2] = x + L*((y+L-1)%L); nn[x+L*y][3] = x + L*((y+1)%L); } // random number generator boost::mt19937 eng(SEED); boost::variate_generator > random_01(eng, boost::uniform_real<>()); // spin configuration std::vector spin(N, 1); int sz = N; // stack for uninspected sites std::stack stck; // connecting probability double pc = 1 - std::exp(-2./T); // measurement alps::ObservableSet measurements; measurements << alps::RealObservable("Magnetization"); measurements << alps::RealObservable("Magnetization^2"); measurements << alps::RealObservable("Magnetization^4"); // timer boost::timer tm; for (int mcs = 0; mcs < MCSTEP + MCTHRM; ++mcs) { if (mcs == MCTHRM) measurements.reset(true); int s = static_cast(random_01() * N); int so = spin[s]; spin[s] = -so; stck.push(s); int cs = 0; while (!stck.empty()) { ++cs; int sc = stck.top(); stck.pop(); for (int k = 0; k < 4; ++k) { int sn = nn[sc][k]; if (spin[sn] == so && random_01() < pc) { stck.push(sn); spin[sn] = -so; } } } sz -= 2 * so * cs; double dsz = sz / static_cast(N); measurements["Magnetization"] << dsz; measurements["Magnetization^2"] << dsz * dsz; measurements["Magnetization^4"] << dsz * dsz * dsz * dsz; } // calculate Binder parameter alps::RealObsevaluator m2 = measurements["Magnetization^2"]; alps::RealObsevaluator m4 = measurements["Magnetization^4"]; alps::RealObsevaluator binder("Binder Ratio of Magnetization"); binder = m2 * m2 / m4; measurements.addObservable(binder); // output results std::cout << measurements; std::cerr << "Elapsed time = " << tm.elapsed() << " sec\n"; return 0; } alpsize-src-20130306-r2942/06_alea/wolff.op0000664000127400003120000000311011071061125016776 0ustar wistarialabBinder Ratio of Magnetization: 0.959454 +/- 0.000919 Magnetization: -0.000368834 +/- 0.00213; tau = -0.381 bin #1 : 32768 entries: error = 0.004371 bin #2 : 16384 entries: error = 0.0027974 bin #3 : 8192 entries: error = 0.00264547 bin #4 : 4096 entries: error = 0.00252408 bin #5 : 2048 entries: error = 0.00239812 bin #6 : 1024 entries: error = 0.00234122 bin #7 : 512 entries: error = 0.00225817 bin #8 : 256 entries: error = 0.00221755 bin #9 : 128 entries: error = 0.00213235 Magnetization^2: 0.626016 +/- 0.00135; tau = 1.29 bin #1 : 32768 entries: error = 0.000710919 bin #2 : 16384 entries: error = 0.000885743 bin #3 : 8192 entries: error = 0.00107123 bin #4 : 4096 entries: error = 0.00123753 bin #5 : 2048 entries: error = 0.0013472 bin #6 : 1024 entries: error = 0.0014023 bin #7 : 512 entries: error = 0.0014018 bin #8 : 256 entries: error = 0.00138314 bin #9 : 128 entries: error = 0.00134558 Magnetization^4: 0.408456 +/- 0.00147; tau = 1.31 bin #1 : 32768 entries: error = 0.000774323 bin #2 : 16384 entries: error = 0.000961241 bin #3 : 8192 entries: error = 0.00116115 bin #4 : 4096 entries: error = 0.00133756 bin #5 : 2048 entries: error = 0.00145559 bin #6 : 1024 entries: error = 0.00152099 bin #7 : 512 entries: error = 0.0015283 bin #8 : 256 entries: error = 0.00150523 bin #9 : 128 entries: error = 0.00147241 alpsize-src-20130306-r2942/06_alea/wolff.ip0000664000127400003120000000001711071061125016773 0ustar wistarialabL = 32 T = 2.2 alpsize-src-20130306-r2942/06_alea/CMakeLists.txt0000664000127400003120000000072112017313126020070 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(alea NONE) # find ALPS Library find_package(ALPS REQUIRED PATHS ${ALPS_ROOT_DIR} $ENV{ALPS_HOME} NO_SYSTEM_ENVIRONMENT_PATH) message(STATUS "Found ALPS: ${ALPS_ROOT_DIR} (revision: ${ALPS_VERSION})") include(${ALPS_USE_FILE}) # enable C++ compilers enable_language(CXX) # rule for generating wolff program add_executable(wolff wolff.C) target_link_libraries(wolff ${ALPS_LIBRARIES}) add_alps_test(wolff) alpsize-src-20130306-r2942/01_original-c/0000775000127400003120000000000012117341400016442 5ustar wistarialabalpsize-src-20130306-r2942/01_original-c/wolff.c0000664000127400003120000000660211364020401017725 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include #define L 32 #define N (L*L) #define T 2.2 #define MCSTEP (1 << 15) #define MCTHRM (MCSTEP >> 3) #define SEED 93812 double random_01() { return (double)(rand()) / RAND_MAX; } int main() { int x, y, k, s, mcs; /* setting up square lattice */ int nn[N][4]; for (y = 0; y < L; ++y) for (x = 0; x < L; ++x) { nn[x+L*y][0] = ((x+L-1)%L) + L*y; nn[x+L*y][1] = ((x+1)%L) + L*y; nn[x+L*y][2] = x + L*((y+L-1)%L); nn[x+L*y][3] = x + L*((y+1)%L); } /* random number generator */ srand(SEED); /* spin configuration */ int spin[N]; for (s = 0; s < N; ++s) spin[s] = 1; int sz = N; /* stack for uninspected sites */ int stck[N]; int is = 0; /* connecting probability */ double pc = 1 - exp(-2./T); /* measurement */ double m = 0; double m2 = 0; double m4 = 0; /* timer */ clock_t tm = clock(); for (mcs = 0; mcs < MCSTEP + MCTHRM; ++mcs) { s = random_01() * N; int so = spin[s]; spin[s] = -so; stck[0] = s; is = 1; int cs = 0; while (is) { ++cs; int sc = stck[--is]; for (k = 0; k < 4; ++k) { int sn = nn[sc][k]; if (spin[sn] == so && random_01() < pc) { stck[is++] = sn; spin[sn] = -so; } } } sz -= 2 * so * cs; if (mcs >= MCTHRM) { double dsz = sz / (double)N; m += dsz; m2 += dsz * dsz; m4 += dsz * dsz * dsz * dsz; } } /* output results */ printf("Magnetization = %f\n", m / MCSTEP); printf("Magnetization^2 = %f\n", m2 / MCSTEP); printf("Magnetization^4 = %f\n", m4 / MCSTEP); printf("Binder Ratio of Magnetization = %f\n", m2 * m2 / m4 / MCSTEP); fprintf(stderr, "Elapsed time = %f sec\n", difftime(clock(), tm) / CLOCKS_PER_SEC); return 0; } alpsize-src-20130306-r2942/01_original-c/CMakeLists.txt0000664000127400003120000000025011364017211021202 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(original-c) # rule for generating wolff program add_executable(wolff wolff.c) target_link_libraries(wolff m) alpsize-src-20130306-r2942/README0000664000127400003120000000101311401430753014776 0ustar wistarialabTutorial: How to ALPSize your applications Copyright (C) 2005-2010 by Synge Todo Folder Contents --------------- - cmake: setting up a source tree using ALPS build framework - original-c: Wolff algorithm written in C - basic-cpp: naive C++ version - stl: standard tempelate library - boost: Boost C++ library - parameters: using ALPS Parameters class - alea: using alea for Monte Carlo averages - lattice: using ALPS Lattice library - scheduler: integrating with ALPS/parapack scheduler alpsize-src-20130306-r2942/00_cmake/0000775000127400003120000000000012117341400015475 5ustar wistarialabalpsize-src-20130306-r2942/00_cmake/hello.C0000664000127400003120000000322311364020401016702 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include int main() { std::cout << "hello, world" << std::endl; } alpsize-src-20130306-r2942/00_cmake/hello.op0000664000127400003120000000001511364013704017142 0ustar wistarialabhello, world alpsize-src-20130306-r2942/00_cmake/CMakeLists.txt0000664000127400003120000000065712017313126020250 0ustar wistarialabcmake_minimum_required(VERSION 2.8 FATAL_ERROR) project(hello NONE) # find ALPS Library find_package(ALPS REQUIRED PATHS ${ALPS_ROOT_DIR} $ENV{ALPS_HOME} NO_SYSTEM_ENVIRONMENT_PATH) message(STATUS "Found ALPS: ${ALPS_ROOT_DIR} (revision: ${ALPS_VERSION})") include(${ALPS_USE_FILE}) # enable C and C++ compilers enable_language(CXX) # rule for generating 'hello world' program add_executable(hello hello.C) add_alps_test(hello) alpsize-src-20130306-r2942/03_stl/0000775000127400003120000000000012117341400015222 5ustar wistarialabalpsize-src-20130306-r2942/03_stl/wolff.C0000664000127400003120000000720211364020401016442 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include #include #include #define L 32 #define N (L*L) #define T 2.2 #define MCSTEP (1 << 15) #define MCTHRM (MCSTEP >> 3) #define SEED 93812 double random_01() { return (double)(std::rand()) / RAND_MAX; } struct vec4 { int& operator[](int i) { return nn_[i]; } int operator[](int i) const { return nn_[i]; } int nn_[4]; }; int main() { // setting up square lattice std::vector nn(N); for (int y = 0; y < L; ++y) for (int x = 0; x < L; ++x) { nn[x+L*y][0] = ((x+L-1)%L) + L*y; nn[x+L*y][1] = ((x+1)%L) + L*y; nn[x+L*y][2] = x + L*((y+L-1)%L); nn[x+L*y][3] = x + L*((y+1)%L); } // random number generator srand(SEED); // spin configuration std::vector spin(N, 1); int sz = N; // stack for uninspected sites std::stack stck; // connecting probability double pc = 1 - std::exp(-2./T); // measurement double m = 0; double m2 = 0; double m4 = 0; // timer std::clock_t tm = std::clock(); for (int mcs = 0; mcs < MCSTEP + MCTHRM; ++mcs) { int s = static_cast(random_01() * N); int so = spin[s]; spin[s] = -so; stck.push(s); int cs = 0; while (!stck.empty()) { ++cs; int sc = stck.top(); stck.pop(); for (int k = 0; k < 4; ++k) { int sn = nn[sc][k]; if (spin[sn] == so && random_01() < pc) { stck.push(sn); spin[sn] = -so; } } } sz -= 2 * so * cs; if (mcs >= MCTHRM) { double dsz = sz / static_cast(N); m += dsz; m2 += dsz * dsz; m4 += dsz * dsz * dsz * dsz; } } // output results std::cout << "Magnetization = " << m / MCSTEP << std::endl; std::cout << "Magnetization^2 = " << m2 / MCSTEP << std::endl; std::cout << "Magnetization^4 = " << m4 / MCSTEP << std::endl; std::cout << "Binder Ratio of Magnetization = " << m2 * m2 / m4 / MCSTEP << std::endl; std::cerr << "Elapsed time = " << std::difftime(std::clock(), tm) / CLOCKS_PER_SEC << " sec\n"; return 0; } alpsize-src-20130306-r2942/03_stl/CMakeLists.txt0000664000127400003120000000020211364017211017757 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(stl) # rule for generating wolff program add_executable(wolff wolff.C) alpsize-src-20130306-r2942/05_parameters/0000775000127400003120000000000012117341400016565 5ustar wistarialabalpsize-src-20130306-r2942/05_parameters/wolff.C0000664000127400003120000000755211364020401020015 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include #include #include #include #include int main() { alps::Parameters params(std::cin); const int L = params.value_or_default("L", 32); const int N = L * L; const double T = params.value_or_default("T", 2.2); const int MCSTEP = params.value_or_default("SWEEPS", 1 << 15); const int MCTHRM = params.value_or_default("THERMALIZATION", MCSTEP >> 3); const unsigned int SEED = params.value_or_default("SEED", 93812); // setting up square lattice std::vector > nn(N); for (int y = 0; y < L; ++y) for (int x = 0; x < L; ++x) { nn[x+L*y][0] = ((x+L-1)%L) + L*y; nn[x+L*y][1] = ((x+1)%L) + L*y; nn[x+L*y][2] = x + L*((y+L-1)%L); nn[x+L*y][3] = x + L*((y+1)%L); } // random number generator boost::mt19937 eng(SEED); boost::variate_generator > random_01(eng, boost::uniform_real<>()); // spin configuration std::vector spin(N, 1); int sz = N; // stack for uninspected sites std::stack stck; // connecting probability double pc = 1 - std::exp(-2./T); // measurement double m = 0; double m2 = 0; double m4 = 0; // timer boost::timer tm; for (int mcs = 0; mcs < MCSTEP + MCTHRM; ++mcs) { int s = static_cast(random_01() * N); int so = spin[s]; spin[s] = -so; stck.push(s); int cs = 0; while (!stck.empty()) { ++cs; int sc = stck.top(); stck.pop(); for (int k = 0; k < 4; ++k) { int sn = nn[sc][k]; if (spin[sn] == so && random_01() < pc) { stck.push(sn); spin[sn] = -so; } } } sz -= 2 * so * cs; if (mcs >= MCTHRM) { double dsz = sz / static_cast(N); m += dsz; m2 += dsz * dsz; m4 += dsz * dsz * dsz * dsz; } } // output results std::cout << "Magnetization = " << m / MCSTEP << std::endl; std::cout << "Magnetization^2 = " << m2 / MCSTEP << std::endl; std::cout << "Magnetization^4 = " << m4 / MCSTEP << std::endl; std::cout << "Binder Ratio of Magnetization = " << m2 * m2 / m4 / MCSTEP << std::endl; std::cerr << "Elapsed time = " << tm.elapsed() << " sec\n"; return 0; } alpsize-src-20130306-r2942/05_parameters/wolff.op0000664000127400003120000000017411071061125020245 0ustar wistarialabMagnetization = -0.000368834 Magnetization^2 = 0.626016 Magnetization^4 = 0.408456 Binder Ratio of Magnetization = 0.959457 alpsize-src-20130306-r2942/05_parameters/wolff.ip0000664000127400003120000000001711071061125020233 0ustar wistarialabL = 32 T = 2.2 alpsize-src-20130306-r2942/05_parameters/CMakeLists.txt0000664000127400003120000000072712017313126021336 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(parameters NONE) # find ALPS Library find_package(ALPS REQUIRED PATHS ${ALPS_ROOT_DIR} $ENV{ALPS_HOME} NO_SYSTEM_ENVIRONMENT_PATH) message(STATUS "Found ALPS: ${ALPS_ROOT_DIR} (revision: ${ALPS_VERSION})") include(${ALPS_USE_FILE}) # enable C++ compilers enable_language(CXX) # rule for generating wolff program add_executable(wolff wolff.C) target_link_libraries(wolff ${ALPS_LIBRARIES}) add_alps_test(wolff) alpsize-src-20130306-r2942/02_basic-cpp/0000775000127400003120000000000012117341400016260 5ustar wistarialabalpsize-src-20130306-r2942/02_basic-cpp/wolff.C0000664000127400003120000000674711364017211017522 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include #define L 32 #define N (L*L) #define T 2.2 #define MCSTEP (1 << 15) #define MCTHRM (MCSTEP >> 3) #define SEED 93812 double random_01() { return (double)(std::rand()) / RAND_MAX; } int main() { // setting up square lattice int nn[N][4]; for (int y = 0; y < L; ++y) for (int x = 0; x < L; ++x) { nn[x+L*y][0] = ((x+L-1)%L) + L*y; nn[x+L*y][1] = ((x+1)%L) + L*y; nn[x+L*y][2] = x + L*((y+L-1)%L); nn[x+L*y][3] = x + L*((y+1)%L); } // random number generator srand(SEED); // spin configuration int spin[N]; for (int s = 0; s < N; ++s) spin[s] = 1; int sz = N; // stack for uninspected sites int stck[N]; int is = 0; // connecting probability double pc = 1 - std::exp(-2./T); // measurement double m = 0; double m2 = 0; double m4 = 0; // timer std::clock_t tm = std::clock(); for (int mcs = 0; mcs < MCSTEP + MCTHRM; ++mcs) { int s = static_cast(random_01() * N); int so = spin[s]; spin[s] = -so; stck[0] = s; is = 1; int cs = 0; while (is) { ++cs; int sc = stck[--is]; for (int k = 0; k < 4; ++k) { int sn = nn[sc][k]; if (spin[sn] == so && random_01() < pc) { stck[is++] = sn; spin[sn] = -so; } } } sz -= 2 * so * cs; if (mcs >= MCTHRM) { double dsz = sz / static_cast(N); m += dsz; m2 += dsz * dsz; m4 += dsz * dsz * dsz * dsz; } } // output results std::cout << "Magnetization = " << m / MCSTEP << std::endl; std::cout << "Magnetization^2 = " << m2 / MCSTEP << std::endl; std::cout << "Magnetization^4 = " << m4 / MCSTEP << std::endl; std::cout << "Binder Ratio of Magnetization = " << m2 * m2 / m4 / MCSTEP << std::endl; std::cerr << "Elapsed time = " << std::difftime(std::clock(), tm) / CLOCKS_PER_SEC << " sec\n"; return 0; } alpsize-src-20130306-r2942/02_basic-cpp/CMakeLists.txt0000664000127400003120000000021011364017211021014 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(basic-cpp) # rule for generating wolff program add_executable(wolff wolff.C) alpsize-src-20130306-r2942/04_boost/0000775000127400003120000000000012117341400015547 5ustar wistarialabalpsize-src-20130306-r2942/04_boost/wolff.C0000664000127400003120000000711411364017211016776 0ustar wistarialab/***************************************************************************** * * Tutorial: How to ALPSize your applications * * Copyright (C) 2005-2010 by Synge Todo * * Permission is hereby granted, free of charge, to any person or organization * obtaining a copy of the software and accompanying documentation covered by * this license (the "Software") to use, reproduce, display, distribute, * execute, and transmit the Software, and to prepare derivative works of the * Software, and to permit third-parties to whom the Software is furnished to * do so, all subject to the following: * * The copyright notices in the Software and this entire statement, including * the above license grant, this restriction and the following disclaimer, * must be included in all copies of the Software, in whole or in part, and * all derivative works of the Software, unless such copies or derivative * works are solely in the form of machine-executable object code generated by * a source language processor. * * 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. * *****************************************************************************/ #include #include #include #include #include #include #include #define L 32 #define N (L*L) #define T 2.2 #define MCSTEP (1 << 15) #define MCTHRM (MCSTEP >> 3) #define SEED 93812 int main() { // setting up square lattice std::vector > nn(N); for (int y = 0; y < L; ++y) for (int x = 0; x < L; ++x) { nn[x+L*y][0] = ((x+L-1)%L) + L*y; nn[x+L*y][1] = ((x+1)%L) + L*y; nn[x+L*y][2] = x + L*((y+L-1)%L); nn[x+L*y][3] = x + L*((y+1)%L); } // random number generator boost::mt19937 eng(SEED); boost::variate_generator > random_01(eng, boost::uniform_real<>()); // spin configuration std::vector spin(N, 1); int sz = N; // stack for uninspected sites std::stack stck; // connecting probability double pc = 1 - std::exp(-2./T); // measurement double m = 0; double m2 = 0; double m4 = 0; // timer boost::timer tm; for (int mcs = 0; mcs < MCSTEP + MCTHRM; ++mcs) { int s = static_cast(random_01() * N); int so = spin[s]; spin[s] = -so; stck.push(s); int cs = 0; while (!stck.empty()) { ++cs; int sc = stck.top(); stck.pop(); for (int k = 0; k < 4; ++k) { int sn = nn[sc][k]; if (spin[sn] == so && random_01() < pc) { stck.push(sn); spin[sn] = -so; } } } sz -= 2 * so * cs; if (mcs >= MCTHRM) { double dsz = sz / static_cast(N); m += dsz; m2 += dsz * dsz; m4 += dsz * dsz * dsz * dsz; } } // output results std::cout << "Magnetization = " << m / MCSTEP << std::endl; std::cout << "Magnetization^2 = " << m2 / MCSTEP << std::endl; std::cout << "Magnetization^4 = " << m4 / MCSTEP << std::endl; std::cout << "Binder Ratio of Magnetization = " << m2 * m2 / m4 / MCSTEP << std::endl; std::cerr << "Elapsed time = " << tm.elapsed() << " sec\n"; return 0; } alpsize-src-20130306-r2942/04_boost/wolff.op0000664000127400003120000000017411071061125017227 0ustar wistarialabMagnetization = -0.000368834 Magnetization^2 = 0.626016 Magnetization^4 = 0.408456 Binder Ratio of Magnetization = 0.959457 alpsize-src-20130306-r2942/04_boost/CMakeLists.txt0000664000127400003120000000064312017313126020315 0ustar wistarialabcmake_minimum_required(VERSION 2.8.0 FATAL_ERROR) project(boost NONE) # find ALPS Library find_package(ALPS REQUIRED PATHS ${ALPS_ROOT_DIR} $ENV{ALPS_HOME} NO_SYSTEM_ENVIRONMENT_PATH) message(STATUS "Found ALPS: ${ALPS_ROOT_DIR} (revision: ${ALPS_VERSION})") include(${ALPS_USE_FILE}) # enable C++ compilers enable_language(CXX) # rule for generating wolff program add_executable(wolff wolff.C) add_alps_test(wolff)