#include #include #include #include namespace kel{ namespace lbm { namespace sch { using namespace saw::schema; } } saw::error_or lbm_main(int argc, char** argv){ (void) argc; (void) argv; using namespace lbm; using namespace acpp; saw::data>> particles{1024u}; for(saw::data i{0u}; i < saw::data{32u}; ++i){ for(saw::data j{0u}; j < saw::data{32u}; ++j){ auto& part = particles.at(i*32ul+j); auto& body = part.template get<"rigid_body">(); auto& pos = body.template get<"position">(); auto& old_pos = body.template get<"position_old">(); auto& acceleration = body.template get<"acceleration">(); auto& p_size = part.template get<"size">(); p_size = {0.5f}; if(j.get() % 2u == 0) acceleration.at({{1u}}) = {-9.81}; pos.at({{0u}}) = {i.template cast_to()}; pos.at({{1u}}) = {j.template cast_to()}; old_pos = pos; } } for(saw::data dt{0u}; dt < saw::data{1024ul*32ul}; ++dt){ for(saw::data i{0u}; i < particles.size(); ++i){ auto& part_i = particles.at(i); verlet_step_lambda(part_i,{1.0f}); } for(saw::data i{0u}; i < particles.size(); ++i){ auto& part_i = particles.at(i); for(saw::data j{i+1ul}; j < particles.size(); ++j){ auto& part_j = particles.at(j); auto res = broadphase_collision_distance(part_i, part_j); if(res.first){ std::cout<<"Collision"<>, saw::encode::Json> j_codec; } return saw::make_void(); } } int main(int argc, char** argv){ auto eov = kel::lbm_main(argc, argv); if(eov.is_error()){ auto& err = eov.get_error(); std::cerr<<"[Error] "< 0u){ std::cerr<<" - "<