mirror of
https://github.com/triqs/dft_tools
synced 2024-12-26 14:23:38 +01:00
70 lines
1.6 KiB
C++
70 lines
1.6 KiB
C++
//#define TRIQS_ARRAYS_ENFORCE_BOUNDCHECK
|
|
|
|
#include <triqs/gf/imfreq.hpp>
|
|
#include <triqs/gf/imtime.hpp>
|
|
#include <triqs/gf/local/functions.hpp>
|
|
|
|
namespace tql= triqs::clef;
|
|
namespace tqa= triqs::arrays;
|
|
using tqa::range;
|
|
using triqs::arrays::make_shape;
|
|
using triqs::gf::Fermion;
|
|
using triqs::gf::imfreq;
|
|
using triqs::gf::imtime;
|
|
using triqs::gf::make_gf;
|
|
|
|
#define TEST(X) std::cout << BOOST_PP_STRINGIZE((X)) << " ---> "<< (X) <<std::endl<<std::endl;
|
|
|
|
const int nl_interne = 1000;
|
|
const int N = 1000;
|
|
|
|
typedef double VALUE_TYPE;
|
|
//typedef int VALUE_TYPE;
|
|
inline VALUE_TYPE fnt(size_t i) { return i*(i+2.0)*(i-8.0);}
|
|
//inline VALUE_TYPE fnt(size_t i) { return i;} //*(i+2.0)*(i-8);}
|
|
//inline VALUE_TYPE fnt(size_t i) { return i*(i+2.0)*(i-8);}
|
|
|
|
struct with_sliding_view {
|
|
void operator()() {
|
|
|
|
double beta =1;
|
|
auto G = make_gf<imfreq> (beta, Fermion, make_shape(2,2),N);
|
|
G() =0;
|
|
|
|
//auto slv = G.data_getter.slv;
|
|
|
|
for (int u =0; u<nl_interne; ++u)
|
|
//for (int i =0; i<N-1; ++i) { G.data_getter.slv[i](0,0) = fnt(i);}
|
|
//for (int i =0; i<N-1; ++i) { slv.set(i); slv(0,0) = fnt(i);}
|
|
//for (int i =0; i<N-1; ++i) G.data_getter(i)(0,0) = fnt(i);
|
|
//for (int i =0; i<N-1; ++i) G.on_mesh(i)(0,0) = fnt(i);
|
|
for (int i =0; i<N-1; ++i) G[i](0,0) = fnt(i);
|
|
|
|
}
|
|
};
|
|
|
|
struct array_code {
|
|
void operator()() {
|
|
|
|
double beta =1;
|
|
auto G = make_gf<imfreq> (beta, Fermion, make_shape(2,2),N);
|
|
G() =0;
|
|
auto V = G.data();
|
|
|
|
for (int u =0; u<nl_interne; ++u)
|
|
for (int i =0; i<N-1; ++i) V(0,0,i) = fnt(i);
|
|
|
|
}
|
|
};
|
|
|
|
|
|
|
|
#include "./speed_tester.hpp"
|
|
int main() {
|
|
speed_tester<array_code> (5000);
|
|
speed_tester<with_sliding_view> (5000);
|
|
//speed_tester<with_slices> (5000);
|
|
return 0;
|
|
}
|
|
|