dc1394’s gists

dc1394 / my_primes_2.cpp

Last active July 7, 2024 02:37

C++20で素数を求めるコードその2

	// コンパイルオプション：g++ -std=c++20 -Wall -O3 -march=native -mtune=native my_primes_2.cpp
	#include <chrono> // for std::chrono
	#include <cmath> // for std::log
	#include <cstdint> // for std::int64_t
	#include <format> // for std::format
	#include <iostream> // for std::iostream
	#include <vector> // for std::vector

	namespace {
	static auto constexpr MAX = 10000000000LL;

dc1394 / my_primes.cpp

Last active July 6, 2024 09:16

C++20とBoost C++ライブラリで素数を求めるコード

	// コンパイルオプション：g++ -std=c++20 -Wall -O3 -march=native -mtune=native my_primes.cpp
	#include <chrono> // for std::chrono
	#include <cmath> // for std::log
	#include <cstdint> // for std::int32_t
	#include <format> // for std::format
	#include <iostream> // for std::iostream
	#include <vector> // for std::vector
	#include <boost/math/quadrature/gauss.hpp> // boost::math::quadrature::gauss

	namespace {

dc1394 / vmc_dsfmt.cpp

Created January 4, 2024 09:21

変分モンテカルロ法のプログラム（dSFMT-AVX512を使ったC++版）

	// 元のコード→ https://miyantarumi.hatenablog.com/entry/2022/04/08/080000
	#include <array> // for std::array
	#include <cmath> // for std::asin, std::cos, std::exp, std::fabs, std::hypot, std::log, std::sqrt, std::sin
	#include <cstdint> // for std::int32_t
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#include <utility> // for std::make_pair, std::pair
	#include <vector> // for std::vector
	#define HAVE_SSE2

dc1394 / vmc_drand48.cpp

Created January 4, 2024 09:13

変分モンテカルロ法のプログラム（C非標準のdrand48()を使ったC++版）

	// 元のコード→ https://miyantarumi.hatenablog.com/entry/2022/04/08/080000
	#include <array> // for std::array
	#include <cmath> // for std::asin, std::cos, std::exp, std::fabs, std::hypot, std::log, std::sqrt, std::sin
	#include <cstdint> // for std::int32_t
	#include <cstdlib> // for drand48, srand48
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#include <utility> // for std::make_pair, std::pair
	#include <vector> // for std::vector

dc1394 / vmc_xorshift.cpp

Last active January 4, 2024 09:17

変分モンテカルロ法のプログラム（Xoshiro256PlusAVX2を使ったC++版）

	// 元のコード→ https://miyantarumi.hatenablog.com/entry/2022/04/08/080000
	#include <array> // for std::array
	#include <cmath> // for std::asin, std::cos, std::exp, std::fabs, std::hypot, std::log, std::sqrt, std::sin
	#include <cstdint> // for std::int32_t
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#include <utility> // for std::make_pair, std::pair
	#include <vector> // for std::vector
	#define __AVX2_AVAILABLE__

dc1394 / vmc_mt19937.cpp

Last active January 4, 2024 09:15

変分モンテカルロ法のプログラム（std::mt19937を使ったC++版）

	// 元のコード→ https://miyantarumi.hatenablog.com/entry/2022/04/08/080000
	#include <algorithm> // for std::fill
	#include <cmath> // for std::asin, std::cos, std::exp, std::fabs, std::hypot, std::log, std::sqrt
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#include <random> // for std::mt19937, std::uniform_real_distribution
	#include <utility> // for std::make_pair, std::pair
	#include <vector> // for std::vector

dc1394 / vmc_crand.cpp

Last active January 4, 2024 09:09

変分モンテカルロ法のプログラム（C備え付けのrand()を使ったC++版）

	// 元のコード→ https://miyantarumi.hatenablog.com/entry/2022/04/08/080000
	#include <array> // for std::array
	#include <cmath> // for std::asin, std::cos, std::exp, std::fabs, std::hypot, std::log, std::sqrt, std::sin
	#include <cstdint> // for std::int32_t
	#include <cstdlib> // for std::rand, std::srand
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#include <utility> // for std::make_pair, std::pair
	#include <vector> // for std::vector

dc1394 / mc_dsfmt_vec.cpp

Last active January 4, 2024 09:19

Twitterのモンテカルロ法のC++版の速度比較コード（dSFMT-AVX512使用＋手動でSIMD (AVX-512)ベクトル化）

	#include <array> // for std::array
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#define HAVE_SSE2
	#define HAVE_AVX2
	#include "dSFMT.h"
	#include "dSFMT-2203-avx512.h"

	namespace {

dc1394 / mc_dsfmt.cpp

Created December 31, 2023 13:32

Twitterのモンテカルロ法のC++版の速度比較コード（dSFMT-AVX512使用）

	#include <array> // for std::array
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#define HAVE_SSE2
	#define HAVE_AVX2
	#include "dSFMT.h"
	#include "dSFMT-2203-avx512.h"

	namespace {

dc1394 / mc_xorshift_vec.cpp

Created December 31, 2023 13:25

Twitterのモンテカルロ法のC++版の速度比較コード（Xoshiro256PlusSIMD使用＋手動でSIMD (AVX2)ベクトル化）

	#include <cstdint> // for std::uint32_t
	#include <iomanip> // for std::setprecision
	#include <ios> // for std::ios::fixed, std::ios::floatfield
	#include <iostream> // for std::cout, std::endl
	#define __AVX2_AVAILABLE__
	#include "SIMDInstructionSet.h"
	#include "Xoshiro256Plus.h"

	namespace {
	inline double mcpi();

dc1394 dc1394