std 27 data/sphere/2014paperALE_sphere/improve_sphere.cpp using namespace std; std 6 data/sphere/first_revision/test.cpp using namespace std; std 26 data/sphere/level_set_test.cpp using Vertex = std::valarray<double>; std 27 data/sphere/level_set_test.cpp using Vertices = std::vector<Vertex>; std 28 data/sphere/level_set_test.cpp const Vertices get_vertices(const std::string& filename); std 41 data/sphere/level_set_test.cpp std::valarray<double> lsf_values (v_list.size()); std 49 data/sphere/level_set_test.cpp std::cout << std::scientific; std 50 data/sphere/level_set_test.cpp std::cout << "Mean value: " << lsf_values.sum() / v_list.size() << std::endl; std 52 data/sphere/level_set_test.cpp std::cout << "Median: "; std 53 data/sphere/level_set_test.cpp std::sort(std::begin(lsf_values), std::end(lsf_values)); std 55 data/sphere/level_set_test.cpp std::cout << lsf_values[ v_list.size()/2 ]; std 57 data/sphere/level_set_test.cpp std::cout << (lsf_values[ v_list.size()/2 - 1 ] + lsf_values[ v_list.size()/2])/2; std 58 data/sphere/level_set_test.cpp std::cout << std::endl; std 60 data/sphere/level_set_test.cpp catch(std::exception& e){ std 61 data/sphere/level_set_test.cpp std::cerr << e.what() << '\n'; std 69 data/sphere/level_set_test.cpp std::ostringstream err_msg; std 72 data/sphere/level_set_test.cpp throw std::runtime_error{err_msg.str()}; std 79 data/sphere/level_set_test.cpp const Vertex convert_to_vertex(const std::string& str_vertex); std 81 data/sphere/level_set_test.cpp const Vertices get_vertices(const std::string& filename){ std 82 data/sphere/level_set_test.cpp std::ifstream dgf_file {filename}; std 83 data/sphere/level_set_test.cpp if(!dgf_file) throw std::runtime_error {"File \"" + filename + "\" does not exist."}; std 85 data/sphere/level_set_test.cpp std::vector<std::string> lines; std 86 data/sphere/level_set_test.cpp for(std::string line; std::getline(dgf_file,line); ) std 89 data/sphere/level_set_test.cpp const std::string begin_vertices = "VERTEX"; std 90 data/sphere/level_set_test.cpp auto first_it = std::find(std::cbegin(lines), std::cend(lines), begin_vertices); std 91 data/sphere/level_set_test.cpp if(first_it == std::cend(lines)) std 92 data/sphere/level_set_test.cpp throw std::runtime_error std 96 data/sphere/level_set_test.cpp const std::string end_vertices = "#"; std 97 data/sphere/level_set_test.cpp const auto second_it = std::find(first_it, std::cend(lines), end_vertices); std 99 data/sphere/level_set_test.cpp throw std::runtime_error std 109 data/sphere/level_set_test.cpp catch(const std::exception& e){ std 110 data/sphere/level_set_test.cpp std::cerr << e.what() << '\n'; std 111 data/sphere/level_set_test.cpp const std::string err_msg {"line" + std::to_string(line_no)}; std 112 data/sphere/level_set_test.cpp throw std::runtime_error {err_msg}; std 120 data/sphere/level_set_test.cpp const Vertex convert_to_vertex(const std::string& str_vertex){ std 121 data/sphere/level_set_test.cpp std::istringstream iss {str_vertex}; std 123 data/sphere/level_set_test.cpp if(!(iss >> d1 >> d2 >> d3) || !(iss >> std::ws).eof() ) std 124 data/sphere/level_set_test.cpp throw std::runtime_error {"convert_to_vertex() failed."}; std 134 data/sphere/level_set_test.cpp using std::begin; std 135 data/sphere/level_set_test.cpp using std::end; std 136 data/sphere/level_set_test.cpp return std::inner_product(begin(vertex), end(vertex), begin(vertex), 0.); std 32 include/esfem_error.h struct Parameter_error : std::runtime_error{ std 33 include/esfem_error.h explicit Parameter_error(const std::string& msg) std 34 include/esfem_error.h : std::runtime_error {"Esfem::Io::Parameter.\n" + msg} std 38 include/esfem_error.h struct SolutionDriven_error : std::runtime_error{ std 39 include/esfem_error.h explicit SolutionDriven_error(const std::string& msg) std 40 include/esfem_error.h : std::runtime_error std 45 include/esfem_error.h struct BrusselatorScheme_error : std::runtime_error{ std 46 include/esfem_error.h explicit BrusselatorScheme_error(const std::string& msg) std 47 include/esfem_error.h : std::runtime_error std 52 include/esfem_error.h struct Grid_error : std::runtime_error{ std 53 include/esfem_error.h explicit Grid_error(const std::string& msg) std 54 include/esfem_error.h : std::runtime_error std 60 include/esfem_error.h struct Rhs_error : std::runtime_error{ std 62 include/esfem_error.h explicit Rhs_error(const std::string& msg) std 63 include/esfem_error.h :std::runtime_error {"Esfem::Rhs.\n" + msg} std 67 include/esfem_error.h struct InitData_error : std::runtime_error{ std 69 include/esfem_error.h explicit InitData_error(const std::string& msg) std 70 include/esfem_error.h :std::runtime_error {"Esfem::InitData.\n" + msg} std 80 include/esfem_error.h (const bool, const std::string&){}; std 83 include/esfem_error.h (const bool, const std::string&){}; std 86 include/esfem_error.h (const bool, const std::string&){}; std 89 include/esfem_error.h (const bool, const std::string&){}; std 92 include/esfem_error.h (const bool, const std::string&){}; std 95 include/esfem_error.h (const bool, const std::string&){}; std 66 include/grid.h explicit Deformation(const std::string& fname); std 79 include/grid.h std::unique_ptr<Data> d_ptr; std 107 include/grid.h explicit Grid_and_time(const Io::Parameter&, const std::string& dgf_file, std 135 include/grid.h std::unique_ptr<Data> d_ptr; std 143 include/grid.h explicit Scal_FEfun(const std::string& name, const Grid_and_time& gt); std 171 include/grid.h explicit Vec_FEfun(const std::string& name, const Grid_and_time&); std 40 include/io_dgf.h Handler(const std::string& dgf_filename); std 43 include/io_dgf.h void write(const std::string& out_filename, const Grid::Scal_FEfun&) const; std 44 include/io_dgf.h void write(const std::string& out_filename, const Grid::Vec_FEfun&) const; std 46 include/io_dgf.h void read(const std::string& in_filename, Grid::Vec_FEfun&) const; std 47 include/io_dgf.h void read(const std::string& in_filename, Grid::Scal_FEfun&) const; std 48 include/io_dgf.h std::size_t no_of_nodes() const noexcept; std 50 include/io_dgf.h using Vertex = std::vector<double>; std 51 include/io_dgf.h using Vertices = std::vector<Vertex>; std 52 include/io_dgf.h using Simplex = std::vector<short>; std 53 include/io_dgf.h using Simplices = std::vector<Simplex>; std 54 include/io_dgf.h using Lines = std::vector<std::string>; std 57 include/io_dgf.h std::unique_ptr<Data> d_ptr; std 63 include/io_dgf.h Handler::Simplices create_simplices(const std::string& filename); std 64 include/io_dgf.h Handler::Lines create_stringSimplices(const std::string& filename); std 66 include/io_dgf.h std::size_t get_precision(const std::string& filename); std 67 include/io_dgf.h std::size_t calculate_dof(const std::string& dgf_file); std 69 include/io_dgf.h Handler::Vertices create_vertices(const std::string& filename); std 70 include/io_dgf.h Handler::Lines create_stringVertices(const std::string& filename); std 79 include/io_dgf.h It find(It first, It last, const std::string& keyword); std 81 include/io_dgf.h It close_list(It first, It last, const std::string& keyword); std 90 include/io_dgf.h void to_dgfFile(const CC&, std::ostream&, const std::size_t precision); std 93 include/io_dgf.h C to(const std::string&); std 95 include/io_dgf.h std::string from(const C&, const std::size_t precision); std 101 include/io_dgf.h to_vertices(It first, It last, const int dim, const std::size_t size); std 107 include/io_dgf.h It find(It first, It last, const std::string& keyword){ std 108 include/io_dgf.h It rv = std::find(first, last, keyword); std 109 include/io_dgf.h if(rv == last) throw std::runtime_error std 115 include/io_dgf.h It close_list(It first, It last, const std::string& keyword){ std 116 include/io_dgf.h It rv = std::find(first, last, keyword); std 117 include/io_dgf.h if(rv == last) throw std::runtime_error std 122 include/io_dgf.h void to_dgfFile(const CC& ccontainer, std::ostream& os, std 123 include/io_dgf.h const std::size_t precision){ std 125 include/io_dgf.h os << from(container, precision) << std::endl; std 128 include/io_dgf.h C to(const std::string& s){ std 130 include/io_dgf.h std::istringstream iss {s}; std 136 include/io_dgf.h std::string from(const C& container, const std::size_t precision){ std 137 include/io_dgf.h std::ostringstream oss {}; std 138 include/io_dgf.h oss << std::scientific; std 140 include/io_dgf.h if( !(oss << std::setprecision(precision) << it << ' ') ) std 141 include/io_dgf.h throw std::runtime_error {"Esfem::Io::Dgf::from()."}; std 146 include/io_dgf.h using std::begin; std 147 include/io_dgf.h using std::end; std 149 include/io_dgf.h if(size % dim) throw std::logic_error std 155 include/io_dgf.h const int dim, const std::size_t size) try{ std 165 include/io_dgf.h rv.emplace_back(std::move(v)); std 169 include/io_dgf.h catch(const std::exception&){ std 170 include/io_dgf.h std::throw_with_nested(std::runtime_error {"Esfem::Io::Dgf::to_vertices()."}); std 175 include/io_dgf.h using std::begin; std 176 include/io_dgf.h using std::end; std 188 include/io_dgf.h if(o_it != end(fef)) throw std::logic_error std 27 include/io_dof.h const std::string& filename); std 30 include/io_dof.h void file_to_dof(It first, const std::string& filename, std::size_t dof_no); std 33 include/io_dof.h const std::string& filename); std 41 include/io_dof.h inline void err_open_file(const std::string& fname){ std 42 include/io_dof.h throw std::logic_error {"Could not open file " + fname}; std 44 include/io_dof.h inline void err_file_to_small(const std::string& fname){ std 45 include/io_dof.h throw std::runtime_error {"File \"" + fname + "\" not big enough."}; std 50 include/io_dof.h const std::string& filename) try{ std 51 include/io_dof.h std::ifstream ifs {filename}; std 60 include/io_dof.h catch(const std::exception&){ std 61 include/io_dof.h std::throw_with_nested(std::runtime_error std 65 include/io_dof.h void file_to_dof(It first, const std::string& filename, std 66 include/io_dof.h std::size_t dof_no) try{ std 67 include/io_dof.h std::ifstream ifs {filename}; std 75 include/io_dof.h if(!(ifs >> std::ws).eof()) throw std::runtime_error std 78 include/io_dof.h catch(const std::exception&){ std 79 include/io_dof.h std::throw_with_nested(std::runtime_error std 84 include/io_dof.h const std::string& filename) try{ std 85 include/io_dof.h std::ofstream ofs {filename}; std 88 include/io_dof.h if( !(ofs << *first << std::endl) ) std 93 include/io_dof.h catch(const std::exception&){ std 94 include/io_dof.h std::throw_with_nested(std::runtime_error std 36 include/io_errorStream.h explicit Error_stream(const std::string& suffix, const Parameter&); std 41 include/io_errorStream.h using Basic_ostream = std::basic_ostream<char, std::char_traits<char> >; std 54 include/io_errorStream.h std::ofstream ofs; std 57 include/io_errorStream.h explicit Error_stream(const std::string& fname); std 53 include/io_l2h1Calculator.h std::unique_ptr<Data> d_ptr; std 31 include/io_parameter.h explicit Parameter(int argc, char** argv, const std::string& parameter_file_name); std 38 include/io_parameter.h const std::string& grid() const noexcept; std 39 include/io_parameter.h const std::string& error_log() const noexcept; std 40 include/io_parameter.h const std::string& paraview() const noexcept; std 65 include/io_parameter.h const std::vector<double>& bdf_alphas() const noexcept; std 66 include/io_parameter.h const std::vector<double>& bdf_gammas() const noexcept; std 86 include/io_parameter.h const std::string& u_init_dof() const noexcept; std 88 include/io_parameter.h const std::string& w_init_dof() const noexcept; std 92 include/io_parameter.h friend std::ostream& operator<<(std::ostream&, const Parameter&); std 96 include/io_parameter.h std::unique_ptr<Data> d_ptr; std 38 include/io_paraview.h std::unique_ptr<Data> d_ptr; std 67 include/secOrd_op_brusselator.h std::unique_ptr<Data> d_ptr; std 46 include/secOrd_op_identity.h std::unique_ptr<Data> d_ptr; std 46 include/secOrd_op_initData.h std::unique_ptr<Data> d_ptr; std 114 include/secOrd_op_initData.h std::unique_ptr<Data> d_ptr; std 84 include/secOrd_op_linearHeat.h std::unique_ptr<Data> d_ptr; std 56 include/secOrd_op_rhs.h std::unique_ptr<Data> d_ptr; std 73 include/secOrd_op_rhs.h std::unique_ptr<Data> d_ptr; std 78 include/secOrd_op_solutionDriven.h std::unique_ptr<Data> d_ptr; std 17 old_code/2015linfty/2015linfty.cpp std::cerr << "Error: " << e << std::endl; std 40 old_code/2015linfty/backup_dune_heat_algo.hpp DataOutputParameters(const std::string name, const int step) std 46 old_code/2015linfty/backup_dune_heat_algo.hpp std::string prefix() const { std 47 old_code/2015linfty/backup_dune_heat_algo.hpp std::stringstream s; std 53 old_code/2015linfty/backup_dune_heat_algo.hpp std::string name_; std 68 old_code/2015linfty/backup_dune_heat_algo.hpp const std::string gridkey = std 70 old_code/2015linfty/backup_dune_heat_algo.hpp const std::string gridfile = std 71 old_code/2015linfty/backup_dune_heat_algo.hpp Dune::Fem::Parameter::getValue< std::string >( gridkey ); std 73 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "Loading macro grid: " << gridfile << std::endl; std 85 old_code/2015linfty/backup_dune_heat_algo.hpp std::ofstream l2h1error_ofs std 86 old_code/2015linfty/backup_dune_heat_algo.hpp {Dune::Fem::Parameter::getValue<std::string>("fem.io.errorFile", std 117 old_code/2015linfty/backup_dune_heat_algo.hpp DataOutputParameters(Dune::Fem::Parameter::getValue<std::string> std 123 old_code/2015linfty/backup_dune_heat_algo.hpp auto write_error = [&](std::ostream& os){ std 124 old_code/2015linfty/backup_dune_heat_algo.hpp os << std::defaultfloat << timeProvider.time() << ' ' std 125 old_code/2015linfty/backup_dune_heat_algo.hpp << std::scientific std 127 old_code/2015linfty/backup_dune_heat_algo.hpp << h1norm.distance(exactSolution, U_n) << std::endl; std 147 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::scientific << "\n==========\n" std 148 old_code/2015linfty/backup_dune_heat_algo.hpp << "U_n.dbegin() = " << *(U_n.dbegin()) << std::endl; std 149 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "rhs.dbegin() = " << *(rhs.dbegin()) << std::endl; std 155 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::defaultfloat << "==========\n" std 157 old_code/2015linfty/backup_dune_heat_algo.hpp << "time = " << timeProvider.time() << std::endl; std 166 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "assembleRHS(f, load_vector) = " << *(load_vector.dbegin()) << std::endl; std 175 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "rhs.axpy = " << *(rhs.dbegin()) << std::endl; std 182 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "solver(rhs, U_n) = " << *(U_n.dbegin()) << std::endl; std 190 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "mass_matrix(U_n, rhs) = " << *(rhs.dbegin()) << std::endl; std 196 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::scientific std 197 old_code/2015linfty/backup_dune_heat_algo.hpp << "U_n.dbegin() = " << *(U_n.dbegin()) << std::endl; std 198 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "rhs.dbegin() = " << *(rhs.dbegin()) << std::endl; std 199 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "load_vector.dbegin() = " << *(load_vector.dbegin()) std 200 old_code/2015linfty/backup_dune_heat_algo.hpp << std::endl; std 207 old_code/2015linfty/backup_dune_heat_algo.hpp throw std::runtime_error("ERROR in BDF::bdf_algorithm():" std 208 old_code/2015linfty/backup_dune_heat_algo.hpp " Bad bdf_no = " + std::to_string(bdf_no) + "\n" std 220 old_code/2015linfty/backup_dune_heat_algo.hpp const std::string gridkey = std 222 old_code/2015linfty/backup_dune_heat_algo.hpp const std::string gridfile = std 223 old_code/2015linfty/backup_dune_heat_algo.hpp Dune::Fem::Parameter::getValue<std::string>(gridkey); std 225 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "Loading macro grid: " << gridfile << std::endl; std 237 old_code/2015linfty/backup_dune_heat_algo.hpp std::ofstream l2h1error_ofs std 238 old_code/2015linfty/backup_dune_heat_algo.hpp {Dune::Fem::Parameter::getValue<std::string>("fem.io.errorFile", std 252 old_code/2015linfty/backup_dune_heat_algo.hpp std::vector<DiscreteFunctionType> prev_M_n_U_n; std 255 old_code/2015linfty/backup_dune_heat_algo.hpp DiscreteFunctionType {"U_n-" + std::to_string(i+1), dfSpace}); std 262 old_code/2015linfty/backup_dune_heat_algo.hpp DataOutputParameters(Dune::Fem::Parameter::getValue<std::string> std 268 old_code/2015linfty/backup_dune_heat_algo.hpp auto write_error = [&](std::ostream& os){ std 269 old_code/2015linfty/backup_dune_heat_algo.hpp os << std::defaultfloat << timeProvider.time() << ' ' std 270 old_code/2015linfty/backup_dune_heat_algo.hpp << std::scientific std 272 old_code/2015linfty/backup_dune_heat_algo.hpp << h1norm.distance(exactSolution, U_n) << std::endl; std 300 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::scientific << "\n==========\n" std 304 old_code/2015linfty/backup_dune_heat_algo.hpp << "prev_M_n_U_n.at(0) = " << *(prev_M_n_U_n.at(0).dbegin()) << std::endl; std 310 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::defaultfloat << "==========\n" std 312 old_code/2015linfty/backup_dune_heat_algo.hpp << "time = " << timeProvider.time() << std::endl; std 321 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "assembleRHS(f, load_vector) = " << *(load_vector.dbegin()) << std::endl; std 330 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "rhs.axpy = " << *(rhs.dbegin()) << std::endl; std 337 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "solver(rhs, U_n) = " << *(U_n.dbegin()) << std::endl; std 345 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "mass_matrix(U_n, rhs) = " << *(rhs.dbegin()) << std::endl; std 351 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::scientific std 352 old_code/2015linfty/backup_dune_heat_algo.hpp << "U_n.dbegin() = " << *(U_n.dbegin()) << std::endl; std 353 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "rhs.dbegin() = " << *(rhs.dbegin()) << std::endl; std 354 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "load_vector.dbegin() = " << *(load_vector.dbegin()) std 355 old_code/2015linfty/backup_dune_heat_algo.hpp << std::endl; std 374 old_code/2015linfty/backup_dune_heat_algo.hpp const std::vector<double> bdf_alpha_coeff { NUMERIK::bdf_alpha_coeff(bdf_no) }; std 375 old_code/2015linfty/backup_dune_heat_algo.hpp const std::vector<double> bdf_gamma_coeff { NUMERIK::bdf_gamma_coeff(bdf_no) }; std 378 old_code/2015linfty/backup_dune_heat_algo.hpp const std::string gridkey = std 380 old_code/2015linfty/backup_dune_heat_algo.hpp const std::string gridfile = std 381 old_code/2015linfty/backup_dune_heat_algo.hpp Dune::Fem::Parameter::getValue< std::string >( gridkey ); std 383 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "Loading macro grid: " << gridfile << std::endl; std 395 old_code/2015linfty/backup_dune_heat_algo.hpp std::ofstream l2h1error_ofs std 396 old_code/2015linfty/backup_dune_heat_algo.hpp {Dune::Fem::Parameter::getValue<std::string>("fem.io.errorFile", std 424 old_code/2015linfty/backup_dune_heat_algo.hpp std::vector<DiscreteFunctionType> prev_steps_vec; std 427 old_code/2015linfty/backup_dune_heat_algo.hpp DiscreteFunctionType {"U_n" + std::to_string(i+1), dfSpace} ); std 450 old_code/2015linfty/backup_dune_heat_algo.hpp DataOutputParameters(Dune::Fem::Parameter::getValue<std::string> std 455 old_code/2015linfty/backup_dune_heat_algo.hpp auto write_error = [&](std::ostream& os){ std 456 old_code/2015linfty/backup_dune_heat_algo.hpp os << std::defaultfloat << timeProvider.time() << ' ' std 457 old_code/2015linfty/backup_dune_heat_algo.hpp << std::scientific std 459 old_code/2015linfty/backup_dune_heat_algo.hpp << h1norm.distance(exact_solution, U_n) << std::endl; std 464 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << '\n' << std::endl; std 470 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::defaultfloat << "==========\n" std 472 old_code/2015linfty/backup_dune_heat_algo.hpp << "time = " << timeProvider.time() << std::endl; std 497 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::scientific std 500 old_code/2015linfty/backup_dune_heat_algo.hpp << "ie_rhs.dbegin() = " << *(ie_rhs.dbegin()) << std::endl; std 541 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::defaultfloat << "==========\n" std 543 old_code/2015linfty/backup_dune_heat_algo.hpp << "time = " << timeProvider.time() << std::endl; std 570 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "assembleRHS(f, ie_load_vector) = " << *(ie_load_vector.dbegin()) << std::endl; std 572 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "ie_rhs.axpy = " << *(ie_rhs.dbegin()) << std::endl; std 574 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "solver(ie_rhs, ie_U_n) = " << *(ie_U_n.dbegin()) << std::endl; std 576 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "mass_matrix(ie_U_n, ie_rhs) = " << *(ie_rhs.dbegin()) << std::endl; std 580 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << std::scientific std 582 old_code/2015linfty/backup_dune_heat_algo.hpp << "ie_U_n.dbegin() = " << *(ie_U_n.dbegin()) << std::endl; std 583 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "rhs.dbegin() = " << *(rhs.dbegin()) << "\t\t" std 584 old_code/2015linfty/backup_dune_heat_algo.hpp << "ie_rhs.dbegin() = " << *(ie_rhs.dbegin()) << std::endl; std 585 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "load_vector.dbegin() = " << *(load_vector.dbegin()) << '\t' std 587 old_code/2015linfty/backup_dune_heat_algo.hpp << std::endl; std 588 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "bdf_tmp.dbegin() = " << *(bdf_tmp.dbegin()) << std::endl; std 591 old_code/2015linfty/backup_dune_heat_algo.hpp std::cout << "prev_steps_vec.at(" << i << ").dbegin() = " std 592 old_code/2015linfty/backup_dune_heat_algo.hpp << *(prev_steps_vec.at(i).dbegin()) << std::endl; std 36 old_code/2015linfty/dune_bdf.hpp EvoMapType(const std::string& input_filename); std 42 old_code/2015linfty/dune_bdf.hpp const std::array<double,3>& operator[] (const std::array<double,3>& vec) const; std 44 old_code/2015linfty/dune_bdf.hpp std::unordered_map<std::string, std::array<double, 3> > evoMap; std 46 old_code/2015linfty/dune_bdf.hpp std::string filename; std 60 old_code/2015linfty/dune_bdf.hpp std::vector<std::string> getVertexList(const std::string& filename) std 65 old_code/2015linfty/dune_bdf.hpp auto error = [](const std::string& msg){ std 67 old_code/2015linfty/dune_bdf.hpp std::runtime_error("ERROR: BDF::getVertexList()\n" std 70 old_code/2015linfty/dune_bdf.hpp auto clean_whitespace = [](const std::string& line) std 73 old_code/2015linfty/dune_bdf.hpp std::string clean_line; std 74 old_code/2015linfty/dune_bdf.hpp std::istringstream iss {line}; std 75 old_code/2015linfty/dune_bdf.hpp for(std::string s; iss >> s; ) std 82 old_code/2015linfty/dune_bdf.hpp std::ifstream ifs {filename}; std 85 old_code/2015linfty/dune_bdf.hpp std::vector<std::string> lines; std 86 old_code/2015linfty/dune_bdf.hpp for(std::string line; getline(ifs, line); ) std 91 old_code/2015linfty/dune_bdf.hpp std::vector<std::string>::iterator first = std 92 old_code/2015linfty/dune_bdf.hpp std::find(lines.begin(), lines.end(), "VERTEX"); std 96 old_code/2015linfty/dune_bdf.hpp std::vector<std::string>::iterator last = std 97 old_code/2015linfty/dune_bdf.hpp std::find(first, lines.end(), "#"); std 101 old_code/2015linfty/dune_bdf.hpp return std::vector<std::string> {first, last}; std 105 old_code/2015linfty/dune_bdf.hpp size_t get_precision(const std::vector<std::string>& vertex_list){ std 108 old_code/2015linfty/dune_bdf.hpp std::vector<size_t> precision_vec; std 109 old_code/2015linfty/dune_bdf.hpp for(const std::string& line : vertex_list){ std 110 old_code/2015linfty/dune_bdf.hpp std::istringstream interpreter {line}; std 111 old_code/2015linfty/dune_bdf.hpp for(std::string s; interpreter >> s; ){ std 116 old_code/2015linfty/dune_bdf.hpp std::vector<size_t>::iterator p_precision = std 117 old_code/2015linfty/dune_bdf.hpp std::max_element(precision_vec.begin(), precision_vec.end()); std 122 old_code/2015linfty/dune_bdf.hpp std::string array_to_string(const std::array<T, N> ad, std 124 old_code/2015linfty/dune_bdf.hpp std::ostringstream oss; std 125 old_code/2015linfty/dune_bdf.hpp oss << std::scientific; std 128 old_code/2015linfty/dune_bdf.hpp for(std::size_t s = 0; s < ad.size() -1 ; ++s) std 129 old_code/2015linfty/dune_bdf.hpp oss << std::setprecision(precision) << ad[s] << ' '; std 135 old_code/2015linfty/dune_bdf.hpp std::array<T, N> string_to_array(const std::string s){ std 136 old_code/2015linfty/dune_bdf.hpp auto error = [] (const std::string& msg) { std 137 old_code/2015linfty/dune_bdf.hpp throw std::runtime_error("ERROR: BDF::string_to_array()\n" + msg); std 139 old_code/2015linfty/dune_bdf.hpp std::istringstream interpreter {s}; std 140 old_code/2015linfty/dune_bdf.hpp std::array<T, N> return_array; std 143 old_code/2015linfty/dune_bdf.hpp if(!(interpreter >> std::ws).eof()) // stuff left in stream? std 150 old_code/2015linfty/dune_bdf.hpp BDF::EvoMapType::EvoMapType( const std::string& input_filename ) : std 153 old_code/2015linfty/dune_bdf.hpp std::vector<std::string> lines = getVertexList(filename); std 161 old_code/2015linfty/dune_bdf.hpp const std::array<double, 3>& std 162 old_code/2015linfty/dune_bdf.hpp BDF::EvoMapType::operator[] (const std::array<double, 3>& node) const { std 174 old_code/2015linfty/dune_bdf.hpp std::array<double, 3> it_new; // translate it back std 175 old_code/2015linfty/dune_bdf.hpp std::copy(x_0.data(), x_0.data() + x_0.size(), it_new.begin()); std 40 old_code/2015linfty/dune_heat_algorithm.hpp DataOutputParameters(const std::string name, const int step) std 46 old_code/2015linfty/dune_heat_algorithm.hpp std::string prefix() const { std 47 old_code/2015linfty/dune_heat_algorithm.hpp std::stringstream s; std 53 old_code/2015linfty/dune_heat_algorithm.hpp std::string name_; std 71 old_code/2015linfty/dune_heat_algorithm.hpp const std::vector<double> bdf_alpha_coeff { NUMERIK::bdf_alpha_coeff(bdf_no) }; std 72 old_code/2015linfty/dune_heat_algorithm.hpp const std::vector<double> bdf_gamma_coeff { NUMERIK::bdf_gamma_coeff(bdf_no) }; std 84 old_code/2015linfty/dune_heat_algorithm.hpp const std::string gridkey = std 86 old_code/2015linfty/dune_heat_algorithm.hpp const std::string gridfile = std 87 old_code/2015linfty/dune_heat_algorithm.hpp Dune::Fem::Parameter::getValue< std::string >( gridkey ); std 89 old_code/2015linfty/dune_heat_algorithm.hpp std::cout << "Loading macro grid: " << gridfile << std::endl; std 111 old_code/2015linfty/dune_heat_algorithm.hpp std::vector<DiscreteFunctionType> prev_steps_U_nmk; // All previous U_{n-k} std 114 old_code/2015linfty/dune_heat_algorithm.hpp {"U_nm" + std::to_string(bdf_no - (i+1)), dfSpace} ); std 115 old_code/2015linfty/dune_heat_algorithm.hpp std::vector<DiscreteFunctionType> prev_steps_M_U_nmk; // All previous M_U_{n-k} std 118 old_code/2015linfty/dune_heat_algorithm.hpp {"M_U_nm" + std::to_string(bdf_no - (i+1)), dfSpace}); std 122 old_code/2015linfty/dune_heat_algorithm.hpp throw std::runtime_error("ERROR in bdf_cycle()."); std 127 old_code/2015linfty/dune_heat_algorithm.hpp std::ofstream l2h1error_ofs std 128 old_code/2015linfty/dune_heat_algorithm.hpp {Dune::Fem::Parameter::getValue<std::string>("fem.io.errorFile", std 130 old_code/2015linfty/dune_heat_algorithm.hpp std::ios_base::app}; std 136 old_code/2015linfty/dune_heat_algorithm.hpp DataOutputParameters(Dune::Fem::Parameter::getValue<std::string> std 154 old_code/2015linfty/dune_heat_algorithm.hpp auto write_error = [&](std::ostream& os) std 158 old_code/2015linfty/dune_heat_algorithm.hpp os << std::defaultfloat << timeProvider.deltaT() << ' ' std 159 old_code/2015linfty/dune_heat_algorithm.hpp << std::scientific std 162 old_code/2015linfty/dune_heat_algorithm.hpp << std::endl; std 187 old_code/2015linfty/dune_heat_algorithm.hpp time_step_no < std::min(prev_steps_U_nmk.size(), size_t(time_step_no_max)); std 28 old_code/2015linfty/iodof.h IO_dune_fem(const std::string& filename_series_prefix, std 34 old_code/2015linfty/iodof.h std::string filename; std 44 old_code/2015linfty/iodof.h IO_dune_fem::IO_dune_fem(const std::string& filename_series_prefix, std 45 old_code/2015linfty/iodof.h const std::string& input_dir, std 46 old_code/2015linfty/iodof.h const std::string& output_dir, std 54 old_code/2015linfty/iodof.h throw std::runtime_error std 62 old_code/2015linfty/iodof.h auto err_handling = [](const std::string& msg){ std 63 old_code/2015linfty/iodof.h std::string full_msg {"Error in IO_dune_fem::operator().\n"}; std 65 old_code/2015linfty/iodof.h throw std::runtime_error {full_msg}; std 68 old_code/2015linfty/iodof.h const std::string curr_file {filename + std::to_string(f_no)}; std 72 old_code/2015linfty/iodof.h std::ofstream ofs {curr_file}; std 75 old_code/2015linfty/iodof.h ofs << *it << std::endl; std 78 old_code/2015linfty/iodof.h std::ifstream ifs {curr_file}; std 73 old_code/tumor_growth/heat.cc std::cerr << "Error: " << e << std::endl; std 36 old_code/tumor_growth/tumor_bdf.h EvoMapType(const std::string& input_filename); std 47 old_code/tumor_growth/tumor_bdf.h const std::array<double,3>& operator[] (const std::array<double,3>& vec) const; std 49 old_code/tumor_growth/tumor_bdf.h std::unordered_map<std::string, std::array<double, 3> > evoMap; std 51 old_code/tumor_growth/tumor_bdf.h std::vector<std::string> original_vertices; std 52 old_code/tumor_growth/tumor_bdf.h std::string filename; std 66 old_code/tumor_growth/tumor_bdf.h std::vector<std::string> getVertexList(const std::string& filename) std 71 old_code/tumor_growth/tumor_bdf.h auto error = [](const std::string& msg){ std 73 old_code/tumor_growth/tumor_bdf.h std::runtime_error("ERROR: BDF::getVertexList()\n" std 76 old_code/tumor_growth/tumor_bdf.h auto clean_whitespace = [](const std::string& line) std 79 old_code/tumor_growth/tumor_bdf.h std::string clean_line; std 80 old_code/tumor_growth/tumor_bdf.h std::istringstream iss {line}; std 81 old_code/tumor_growth/tumor_bdf.h for(std::string s; iss >> s; ) std 88 old_code/tumor_growth/tumor_bdf.h std::ifstream ifs {filename}; std 91 old_code/tumor_growth/tumor_bdf.h std::vector<std::string> lines; std 92 old_code/tumor_growth/tumor_bdf.h for(std::string line; getline(ifs, line); ) std 97 old_code/tumor_growth/tumor_bdf.h std::vector<std::string>::iterator first = std 98 old_code/tumor_growth/tumor_bdf.h std::find(lines.begin(), lines.end(), "VERTEX"); std 102 old_code/tumor_growth/tumor_bdf.h std::vector<std::string>::iterator last = std 103 old_code/tumor_growth/tumor_bdf.h std::find(first, lines.end(), "#"); std 107 old_code/tumor_growth/tumor_bdf.h return std::vector<std::string> {first, last}; std 111 old_code/tumor_growth/tumor_bdf.h size_t get_relevant_digits(const std::string& number){ std 120 old_code/tumor_growth/tumor_bdf.h size_t get_precision(const std::vector<std::string>& vertex_list){ std 123 old_code/tumor_growth/tumor_bdf.h std::vector<size_t> precision_vec; std 124 old_code/tumor_growth/tumor_bdf.h for(const std::string& line : vertex_list){ std 125 old_code/tumor_growth/tumor_bdf.h std::istringstream interpreter {line}; std 126 old_code/tumor_growth/tumor_bdf.h for(std::string s; interpreter >> s; ){ std 137 old_code/tumor_growth/tumor_bdf.h std::vector<size_t>::iterator p_precision = std 138 old_code/tumor_growth/tumor_bdf.h std::max_element(precision_vec.begin(), precision_vec.end()); std 143 old_code/tumor_growth/tumor_bdf.h std::string array_to_string(const std::array<T, N> ad, std 145 old_code/tumor_growth/tumor_bdf.h std::ostringstream oss; std 146 old_code/tumor_growth/tumor_bdf.h oss << std::scientific; std 149 old_code/tumor_growth/tumor_bdf.h for(std::size_t s = 0; s < ad.size() -1 ; ++s) std 150 old_code/tumor_growth/tumor_bdf.h oss << std::setprecision(precision) << ad[s] << ' '; std 156 old_code/tumor_growth/tumor_bdf.h std::array<T, N> string_to_array(const std::string s){ std 157 old_code/tumor_growth/tumor_bdf.h auto error = [] (const std::string& msg) { std 158 old_code/tumor_growth/tumor_bdf.h throw std::runtime_error("ERROR: BDF::string_to_array()\n" + msg); std 160 old_code/tumor_growth/tumor_bdf.h std::istringstream interpreter {s}; std 161 old_code/tumor_growth/tumor_bdf.h std::array<T, N> return_array; std 164 old_code/tumor_growth/tumor_bdf.h if(!(interpreter >> std::ws).eof()) // stuff left in stream? std 171 old_code/tumor_growth/tumor_bdf.h BDF::EvoMapType::EvoMapType( const std::string& input_filename ) : std 174 old_code/tumor_growth/tumor_bdf.h std::vector<std::string> lines = getVertexList(filename); std 182 old_code/tumor_growth/tumor_bdf.h const std::array<double, 3>& std 183 old_code/tumor_growth/tumor_bdf.h BDF::EvoMapType::operator[] (const std::array<double, 3>& node) const { std 184 old_code/tumor_growth/tumor_bdf.h const std::array<double, 3>* cpad3 = nullptr; std 189 old_code/tumor_growth/tumor_bdf.h catch(std::exception& e){ std 192 old_code/tumor_growth/tumor_bdf.h std::cout << "Your node\n" std 194 old_code/tumor_growth/tumor_bdf.h + "\nWas not good. Full node list." << std::endl; std 196 old_code/tumor_growth/tumor_bdf.h std::cout << it.first << std::endl; std 210 old_code/tumor_growth/tumor_bdf.h std::array<double, 3> it_new; // translate it back std 211 old_code/tumor_growth/tumor_bdf.h std::copy(x_0.data(), x_0.data() + x_0.size(), it_new.begin()); std 243 old_code/tumor_growth/tumor_bdf.h if(no_vertices * dim != dof_times_dim) throw std::runtime_error std 250 old_code/tumor_growth/tumor_bdf.h std::array<double, dim> array_new_vertex; std 40 old_code/tumor_growth/tumor_code.h DataOutputParameters(const std::string name, const int step) std 46 old_code/tumor_growth/tumor_code.h std::string prefix() const { std 47 old_code/tumor_growth/tumor_code.h std::stringstream s; std 53 old_code/tumor_growth/tumor_code.h std::string name_; std 70 old_code/tumor_growth/tumor_code.h const std::vector<double> bdf_alpha_coeff { NUMERIK::bdf_alpha_coeff(bdf_no) }; std 71 old_code/tumor_growth/tumor_code.h const std::vector<double> bdf_gamma_coeff { NUMERIK::bdf_gamma_coeff(bdf_no) }; std 83 old_code/tumor_growth/tumor_code.h const std::string gridkey = std 85 old_code/tumor_growth/tumor_code.h const std::string gridfile = std 86 old_code/tumor_growth/tumor_code.h Dune::Fem::Parameter::getValue< std::string >( gridkey ); std 88 old_code/tumor_growth/tumor_code.h std::cout << "Loading macro grid: " << gridfile << std::endl; std 110 old_code/tumor_growth/tumor_code.h std::vector<DiscreteFunctionType> prev_steps_U_nmk; // All previous U_{n-k} std 113 old_code/tumor_growth/tumor_code.h {"U_nm" + std::to_string(bdf_no - (i+1)), dfSpace} ); std 114 old_code/tumor_growth/tumor_code.h std::vector<DiscreteFunctionType> prev_steps_M_U_nmk; // All previous M_U_{n-k} std 117 old_code/tumor_growth/tumor_code.h {"M_U_nm" + std::to_string(bdf_no - (i+1)), dfSpace}); std 121 old_code/tumor_growth/tumor_code.h throw std::runtime_error("ERROR in bdf_cycle()."); std 126 old_code/tumor_growth/tumor_code.h std::ofstream l2h1error_ofs std 127 old_code/tumor_growth/tumor_code.h {Dune::Fem::Parameter::getValue<std::string>("fem.io.errorFile", std 129 old_code/tumor_growth/tumor_code.h std::ios_base::app}; std 135 old_code/tumor_growth/tumor_code.h DataOutputParameters(Dune::Fem::Parameter::getValue<std::string> std 153 old_code/tumor_growth/tumor_code.h auto write_error = [&](std::ostream& os) std 157 old_code/tumor_growth/tumor_code.h os << std::defaultfloat << timeProvider.deltaT() << ' ' std 158 old_code/tumor_growth/tumor_code.h << std::scientific std 161 old_code/tumor_growth/tumor_code.h << std::endl; std 186 old_code/tumor_growth/tumor_code.h time_step_no < std::min(prev_steps_U_nmk.size(), size_t(time_step_no_max)); std 390 old_code/tumor_growth/tumor_code.h const std::string gridkey = std 392 old_code/tumor_growth/tumor_code.h const std::string gridfile = std 393 old_code/tumor_growth/tumor_code.h Dune::Fem::Parameter::getValue< std::string >( gridkey ); std 395 old_code/tumor_growth/tumor_code.h std::cout << "Loading macro grid: " << gridfile << std::endl; std 494 old_code/tumor_growth/tumor_code.h DataOutputParameters(Dune::Fem::Parameter::getValue<std::string> std 500 old_code/tumor_growth/tumor_code.h DataOutputParameters(Dune::Fem::Parameter::getValue<std::string> std 44 old_code/tumor_growth/tumor_deformation.h std::array<double, 3> domain_node { {x[0], x[1], x[2]} }; std 46 old_code/tumor_growth/tumor_deformation.h const std::array<double, 3> range_node = (*evoMapPt)[domain_node]; std 47 old_code/tumor_growth/tumor_deformation.h for(std::size_t s = 0; s < range_node.size(); ++s) std 83 old_code/tumor_growth/tumor_growth.h : random_fun (std::bind(Random_dist {mean_value, variance}, std 94 old_code/tumor_growth/tumor_growth.h using Random_dist = std::normal_distribution<>; std 95 old_code/tumor_growth/tumor_growth.h using Random_engine = std::default_random_engine; std 96 old_code/tumor_growth/tumor_growth.h std::function<double()> random_fun; std 440 old_code/tumor_growth/tumor_growth.h return std::sqrt(vv[0] * vv[0] + vv[1] * vv[1] + vv[2] * vv[2]); std 475 old_code/tumor_growth/tumor_growth.h throw std::runtime_error {"Error in constructor of Brusselator_model.\n" std 72 src/brusselator_algo.cpp const std::string& parameter_fname) std 148 src/brusselator_algo.cpp << std::endl; std 170 src/brusselator_algo.cpp std::unique_ptr<SecOrd_op::vIdata> ex_ptr std 186 src/brusselator_algo.cpp << std::endl; std 193 src/brusselator_algo.cpp using std::unique_ptr; std 215 src/brusselator_algo.cpp << std::endl; std 222 src/brusselator_algo.cpp std::unique_ptr<vRhs> vRhs_ptr {vRhs::new_sd(fix_grid)}; std 223 src/brusselator_algo.cpp std::unique_ptr<vIdata> ex_ptr {vIdata::new_sd(fix_grid)}; std 237 src/brusselator_algo.cpp << std::endl; std 244 src/brusselator_algo.cpp using std::unique_ptr; std 375 src/brusselator_algo.cpp << norm.h1_err(fef.u.exact, fef.u.fun) << std::endl; std 378 src/brusselator_algo.cpp << norm.h1_err(fef.w.exact, fef.w.fun) << std::endl; std 384 src/brusselator_algo.cpp << std::endl; std 390 src/brusselator_algo.cpp << std::endl; std 404 src/brusselator_algo.cpp io.para << data << std::endl; std 416 src/brusselator_algo.cpp std::cout << "t0: " << data.start_time() << '\n' std 419 src/brusselator_algo.cpp << "last dT: " << data.last_step() << std::endl; std 420 src/brusselator_algo.cpp std::cout << "time: " << fix_grid.time_provider().time() << std::endl; std 423 src/brusselator_algo.cpp std::cout << "it < it_end-1" << std::endl; std 427 src/brusselator_algo.cpp std::cout << "it >= it_end-1" << std::endl; std 430 src/brusselator_algo.cpp std::cout << "time: " << fix_grid.time_provider().time() << std::endl; std 150 src/brusselator_algo.h const std::string& parameter_fname); std 355 src/brusselator_algo.h std::unique_ptr<SecOrd_op::vIdata> X_ptr; std 374 src/brusselator_algo.h const std::string tmpFile_path {FEF_PATH}; std 426 src/brusselator_algo.h << norm.h1_err(fem.fun, fem.exact) << std::endl; std 157 src/brusselator_algo_impl.cpp std::throw_with_nested(Bruss_error {"RhsAndSolve_helper()"}); std 232 src/brusselator_algo_impl.cpp << "H1err" << std::endl; std 233 src/brusselator_algo_impl.cpp file << std::scientific; std 8 src/brusselator_main.cpp void print_errMsg(const std::exception&); std 15 src/brusselator_main.cpp std::cerr << "Dune error: " << e << std::endl; std 18 src/brusselator_main.cpp catch(const std::exception& e){ std 26 src/brusselator_main.cpp void print_errMsg(const std::exception& e){ std 27 src/brusselator_main.cpp std::cerr << e.what() << '\n'; std 29 src/brusselator_main.cpp std::rethrow_if_nested(e); std 31 src/brusselator_main.cpp catch(const std::exception& nested){ std 453 src/config/config.h #define SHARED_PTR_NAMESPACE std std 453 src/config/config_osx.h #define SHARED_PTR_NAMESPACE std std 454 src/config/config_ubuntu.h #define SHARED_PTR_NAMESPACE std std 43 src/grid_FEfunSet.h r3fef(const std::string& name, const Grid_and_time& gt) std 52 src/grid_FEfunSet.h FEfun_set(const std::string& name, const Grid::Grid_and_time&); std 63 src/grid_FEfunSet.h const std::string& dir = "/tmp/"); std 66 src/grid_FEfunSet.h const std::string& dir = "/tmp/"); std 98 src/grid_FEfunSet.h inline std::string compose_dgfName(const std::string& fun_name, std 99 src/grid_FEfunSet.h const std::string& dir = "./"); std 115 src/grid_FEfunSet.h FEfun_set(const std::string& name, std 130 src/grid_FEfunSet.h const std::string& dir){ std 138 src/grid_FEfunSet.h read(const Esfem::Io::Dgf::Handler& h, const std::string& dir){ std 155 src/grid_FEfunSet.h inline std::string compose_dgfName(const std::string& fun_name, std 156 src/grid_FEfunSet.h const std::string& dir){ std 26 src/grid_GridAndTime.cpp using namespace std; std 53 src/grid_GridAndTime.cpp Data(const Io::Parameter&, const std::string& dgf_file, std 68 src/grid_GridAndTime.cpp Data(const Io::Parameter& p, const std::string& dgf_file, std 81 src/grid_GridAndTime.cpp try :d_ptr {std::make_unique<Data>(p)} std 86 src/grid_GridAndTime.cpp << "Dune error: " << e << std::endl; std 94 src/grid_GridAndTime.cpp Grid_and_time(const Io::Parameter& p, const std::string& dgf_file, const double t0) std 95 src/grid_GridAndTime.cpp try : d_ptr {std::make_unique<Data>(p, dgf_file, t0)} std 97 src/grid_GridAndTime.cpp catch(const std::exception&){ std 103 src/grid_GridAndTime.cpp << "Dune error: " << e << std::endl; std 20 src/grid_GridAndTime_impl.cpp using namespace std; std 31 src/grid_GridAndTime_impl.cpp Esfem::Impl::Evolving_grid::Evolving_grid(const std::string& filename) std 36 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 37 src/grid_GridAndTime_impl.cpp std::throw_with_nested std 53 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 54 src/grid_GridAndTime_impl.cpp std::throw_with_nested std 63 src/grid_GridAndTime_impl.cpp using std::cbegin; std 64 src/grid_GridAndTime_impl.cpp using std::cend; std 66 src/grid_GridAndTime_impl.cpp constexpr std::size_t dim = 3; std 74 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 75 src/grid_GridAndTime_impl.cpp std::throw_with_nested std 86 src/grid_GridAndTime_impl.cpp get_vertexList(const std::string& filename) try{ std 87 src/grid_GridAndTime_impl.cpp using std::cbegin; std 88 src/grid_GridAndTime_impl.cpp using std::cend; std 95 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 96 src/grid_GridAndTime_impl.cpp std::throw_with_nested std 108 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 109 src/grid_GridAndTime_impl.cpp std::throw_with_nested std 118 src/grid_GridAndTime_impl.cpp identity_map(const Nodes_key& nodes_lines, const std::size_t precision){ std 126 src/grid_GridAndTime_impl.cpp std::size_t Esfem::Impl::get_relevant_precision(const std::string& number) try{ std 133 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 134 src/grid_GridAndTime_impl.cpp std::throw_with_nested(logic_error std 155 src/grid_GridAndTime_impl.cpp std::size_t Esfem::Impl::max(const std::vector<std::size_t>& v){ std 161 src/grid_GridAndTime_impl.cpp std::string Esfem::Impl::node_to_string(const Node& node, const int precision) try{ std 163 src/grid_GridAndTime_impl.cpp std::ostringstream oss; std 164 src/grid_GridAndTime_impl.cpp oss << std::scientific; std 165 src/grid_GridAndTime_impl.cpp for(std::size_t i = 0; i < Node::dimension - 1; ++i) std 166 src/grid_GridAndTime_impl.cpp if( !(oss << std::setprecision(precision) << node[i] << ' ') ) std 168 src/grid_GridAndTime_impl.cpp oss << std::setprecision(precision) << node[Node::dimension - 1]; std 177 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 178 src/grid_GridAndTime_impl.cpp std::throw_with_nested(runtime_error {"Error in node_to_string()."}); std 184 src/grid_GridAndTime_impl.cpp Node Esfem::Impl::string_to_node(const std::string& str_node) try{ std 185 src/grid_GridAndTime_impl.cpp std::istringstream interpreter {str_node}; std 188 src/grid_GridAndTime_impl.cpp for(std::size_t it = 0; it < Node::dimension; ++it){ std 192 src/grid_GridAndTime_impl.cpp if( !(interpreter >> std::ws).eof()) // stuff left in stream std 196 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 197 src/grid_GridAndTime_impl.cpp std::throw_with_nested std 204 src/grid_GridAndTime_impl.cpp std::string Esfem::Impl:: std 205 src/grid_GridAndTime_impl.cpp normalize_node(const std::string& str_node, const int precision){ std 209 src/grid_GridAndTime_impl.cpp std::string Esfem::Impl::clean_whitespace(const std::string& line){ std 210 src/grid_GridAndTime_impl.cpp std::string clean_line {}; std 211 src/grid_GridAndTime_impl.cpp std::istringstream iss {line}; std 212 src/grid_GridAndTime_impl.cpp for(std::string s; iss >> s; ) std 218 src/grid_GridAndTime_impl.cpp std::vector<std::string> Esfem::Impl:: std 219 src/grid_GridAndTime_impl.cpp dgfFile_to_vec(const std::string& filename) try{ std 220 src/grid_GridAndTime_impl.cpp std::ifstream dgf_file {filename}; std 223 src/grid_GridAndTime_impl.cpp std::vector<std::string> lines; std 224 src/grid_GridAndTime_impl.cpp for(std::string line; std::getline(dgf_file, line); ) std 229 src/grid_GridAndTime_impl.cpp catch(const std::exception&){ std 230 src/grid_GridAndTime_impl.cpp std::throw_with_nested std 240 src/grid_GridAndTime_impl.cpp grid::grid(const std::string& fname){ std 305 src/grid_GridAndTime_impl.cpp std::string Esfem::Impl::hash::to_string(const range& k){ std 31 src/grid_GridAndTime_impl.h template<> struct std::hash<Esfem::Grid::Deformation::Domain>{ std 34 src/grid_GridAndTime_impl.h std::size_t operator()(const Esfem::Grid::Deformation::Domain& n) const{ std 35 src/grid_GridAndTime_impl.h using std::hash; std 49 src/grid_GridAndTime_impl.h template<> struct std::hash<array<int, Esfem::Grid::world_dim()> >{ std 54 src/grid_GridAndTime_impl.h std::size_t operator()(const Base& ai) const{ std 68 src/grid_GridAndTime_impl.h template<> struct std::equal_to<Esfem::Grid::Deformation::Domain>{ std 80 src/grid_GridAndTime_impl.h template<> struct std::equal_to<array<int, Esfem::Grid::world_dim()> >{ std 105 src/grid_GridAndTime_impl.h explicit Evolving_grid(const std::string& filename); std 113 src/grid_GridAndTime_impl.h using Nodes_key = std::vector<std::string>; std 115 src/grid_GridAndTime_impl.h using Map = std::unordered_map<std::string, Node>; std 120 src/grid_GridAndTime_impl.h std::size_t digit_precision {8}; std 135 src/grid_GridAndTime_impl.h using key = std::array<int, range::dimension>; std 143 src/grid_GridAndTime_impl.h using map = std::unordered_map<key, range>; std 145 src/grid_GridAndTime_impl.h using seq = std::vector<key>; std 152 src/grid_GridAndTime_impl.h struct bad : std::runtime_error{ std 154 src/grid_GridAndTime_impl.h explicit bad(const std::string& msg) :std::runtime_error {msg} {} std 159 src/grid_GridAndTime_impl.h explicit grid(const std::string& fname); std 170 src/grid_GridAndTime_impl.h std::string to_string(const range&); std 178 src/grid_GridAndTime_impl.h Evolving_grid::Nodes_key get_vertexList(const std::string& filename); std 181 src/grid_GridAndTime_impl.h std::size_t get_relevant_precision(const Evolving_grid::Nodes_key& vertex_list); std 184 src/grid_GridAndTime_impl.h identity_map(const Evolving_grid::Nodes_key&, const std::size_t precision); std 186 src/grid_GridAndTime_impl.h std::size_t get_relevant_precision(const std::string& number); std 189 src/grid_GridAndTime_impl.h std::vector<std::size_t> std 190 src/grid_GridAndTime_impl.h create_precision_vec(const std::vector<std::string>& vertex_list); std 194 src/grid_GridAndTime_impl.h std::size_t max(const std::vector<std::size_t>&); std 196 src/grid_GridAndTime_impl.h std::string node_to_string(const Evolving_grid::Node&, const int precision); std 197 src/grid_GridAndTime_impl.h Evolving_grid::Node string_to_node(const std::string&); std 198 src/grid_GridAndTime_impl.h std::string normalize_node(const std::string&, const int precision); std 200 src/grid_GridAndTime_impl.h std::string clean_whitespace(const std::string&); std 202 src/grid_GridAndTime_impl.h std::vector<std::string> dgfFile_to_vec(const std::string& filename); std 205 src/grid_GridAndTime_impl.h It dgf_find_vertex(It first, It last, const std::string& filename); std 207 src/grid_GridAndTime_impl.h It dgf_close_list(It first, It last, const std::string& filename); std 225 src/grid_GridAndTime_impl.h throw std::logic_error {err_msg}; std 232 src/grid_GridAndTime_impl.h It dgf_find_vertex(It first, It last, const std::string& filename){ std 233 src/grid_GridAndTime_impl.h std::string starting_point {"VERTEX"}; std 234 src/grid_GridAndTime_impl.h It rv = std::find(first, last, "VERTEX"); std 236 src/grid_GridAndTime_impl.h throw std::runtime_error std 243 src/grid_GridAndTime_impl.h It dgf_close_list(It first, It last, const std::string& filename){ std 244 src/grid_GridAndTime_impl.h std::string ending_point {"#"}; std 245 src/grid_GridAndTime_impl.h It rv = std::find(first, last, ending_point); std 247 src/grid_GridAndTime_impl.h throw std::runtime_error std 26 src/grid_deformation.cpp using namespace std; std 71 src/grid_deformation.cpp = std::inner_product(&x[0], &x[0] + Domain::dimension, &x[0], 0.); std 90 src/grid_deformation.cpp Data(const std::string& fname) :hg {fname} {} std 93 src/grid_deformation.cpp Esfem::Grid::Deformation::Deformation() :d_ptr {std::make_unique<Data>()} {} std 94 src/grid_deformation.cpp Esfem::Grid::Deformation::Deformation(const std::string& fname) std 95 src/grid_deformation.cpp :d_ptr {std::make_unique<Data>(fname)} {} std 26 src/grid_fef.cpp Scal_FEfun::Scal_FEfun(const std::string& fun_name, const Grid_and_time& gt) std 52 src/grid_fef.cpp Vec_FEfun::Vec_FEfun(const std::string& fun_name, const Grid_and_time& gt) std 23 src/io_dgf.cpp using namespace std; std 77 src/io_dgf.cpp Handler::Handler(const std::string& dgf_filename) std 90 src/io_dgf.cpp void Handler::write(const std::string& out_filename, const Scal_FEfun& fef) const try{ std 97 src/io_dgf.cpp dgf_file << "DGF" << std::endl; std 99 src/io_dgf.cpp dgf_file << "VERTEX" << std::endl; std 102 src/io_dgf.cpp << "SIMPLEX" << std::endl; std 104 src/io_dgf.cpp dgf_file << "#" << std::endl; std 107 src/io_dgf.cpp << "\n#\n#" << std::endl; std 113 src/io_dgf.cpp void Handler::write(const std::string& out_filename, const Vec_FEfun& vfef) const try{ std 120 src/io_dgf.cpp dgf_file << "DGF" << std::endl; std 122 src/io_dgf.cpp dgf_file << "VERTEX" << std::endl; std 125 src/io_dgf.cpp << "SIMPLEX" << std::endl; std 127 src/io_dgf.cpp dgf_file << "#" << std::endl; std 130 src/io_dgf.cpp << "\n#\n#" << std::endl; std 187 src/io_dgf.cpp Simplices Esfem::Io::Dgf::create_simplices(const std::string& filename) try{ std 30 src/io_errorStream.cpp void cannot_open_file(const std::string& fname); std 32 src/io_errorStream.cpp Esfem::Io::Error_stream::Error_stream(const std::string& fname) std 33 src/io_errorStream.cpp :ofs {fname, std::ios_base::app} std 36 src/io_errorStream.cpp std::clog << "Opening " << fname << " for Error_stream" << std::endl; std 43 src/io_errorStream.cpp Esfem::Io::Error_stream::Error_stream(const std::string& suffix, const Parameter& d) std 48 src/io_errorStream.cpp std::clog << "Using operator<<(StdManip)." << std::endl; std 61 src/io_errorStream.cpp void cannot_open_file(const std::string& fname){ std 62 src/io_errorStream.cpp std::ostringstream err_msg; std 65 src/io_errorStream.cpp throw std::runtime_error {err_msg.str()}; std 25 src/io_l2h1Calculator.cpp using namespace std; std 26 src/io_parameter.cpp using namespace std; std 37 src/io_parameter.cpp const std::string grid_dgf; std 38 src/io_parameter.cpp const std::string error_log; std 39 src/io_parameter.cpp const std::string paraview; std 41 src/io_parameter.cpp const std::vector<double> bdf_alphas; std 42 src/io_parameter.cpp const std::vector<double> bdf_gammas; std 55 src/io_parameter.cpp const std::string u_init_dof; std 56 src/io_parameter.cpp const std::string w_init_dof; std 66 src/io_parameter.cpp error_log {Dune::Fem::Parameter::getValue<std::string> std 68 src/io_parameter.cpp paraview {Dune::Fem::Parameter::getValue<std::string> std 92 src/io_parameter.cpp u_init_dof {Dune::Fem::Parameter::getValue<std::string> std 94 src/io_parameter.cpp w_init_dof {Dune::Fem::Parameter::getValue<std::string> std 110 src/io_parameter.cpp const std::string& parameter_file_name){ std 160 src/io_parameter.cpp const std::string& Esfem::Io::Parameter::grid() const noexcept{ std 163 src/io_parameter.cpp const std::string& Esfem::Io::Parameter::error_log() const noexcept{ std 166 src/io_parameter.cpp const std::string& Esfem::Io::Parameter::paraview() const noexcept{ std 197 src/io_parameter.cpp const std::vector<double>& Esfem::Io::Parameter::bdf_alphas() const noexcept{ std 200 src/io_parameter.cpp const std::vector<double>& Esfem::Io::Parameter::bdf_gammas() const noexcept{ std 236 src/io_parameter.cpp const std::string& Esfem::Io::Parameter::u_init_dof() const noexcept{ std 239 src/io_parameter.cpp const std::string& Esfem::Io::Parameter::w_init_dof() const noexcept{ std 243 src/io_parameter.cpp std::ostream& Esfem::Io::operator<<(std::ostream& os, const Parameter& d){ std 30 src/io_parameter_impl.cpp using namespace std; std 32 src/io_parameter_impl.cpp const std::string& Esfem::Impl::project_dir(){ std 34 src/io_parameter_impl.cpp static const std::string project_dir {"/home/power/cpp/DISS_surfaces/"}; std 38 src/io_parameter_impl.cpp void Esfem::Impl::dune_fem_parameter_append(int argc, char** argv, const std::string& file){ std 44 src/io_parameter_impl.cpp std::clog << "Using parameter file: " << file << std::endl; std 48 src/io_parameter_impl.cpp std::string Esfem::Impl::get_gridKey(){ std 52 src/io_parameter_impl.cpp std::clog << "grid_key: " << grid_key << std::endl; std 56 src/io_parameter_impl.cpp std::string Esfem::Impl::get_macroGrid(){ std 58 src/io_parameter_impl.cpp const auto macro_grid = Dune::Fem::Parameter::getValue<std::string>(grid_key); std 60 src/io_parameter_impl.cpp std::clog << "Loading macro grid: " << macro_grid << std::endl; std 63 src/io_parameter_impl.cpp std::string Esfem::Impl::doubleVector_to_string(const std::vector<double>& vd){ std 69 src/io_parameter_impl.cpp std::ostringstream oss; std 71 src/io_parameter_impl.cpp for(std::size_t it = 0; it < vd.size() - 1; ++it) std 76 src/io_parameter_impl.cpp void Esfem::Impl::file_check(const std::vector<std::string>& file_list){ std 25 src/io_parameter_impl.h const std::string& project_dir(); std 26 src/io_parameter_impl.h void dune_fem_parameter_append(int argc, char** argv, const std::string& file); std 27 src/io_parameter_impl.h std::string get_gridKey(); std 28 src/io_parameter_impl.h std::string get_macroGrid(); std 29 src/io_parameter_impl.h std::string doubleVector_to_string(const std::vector<double>&); std 30 src/io_parameter_impl.h void file_check(const std::vector<std::string>& file_list); std 23 src/io_paraview.cpp using namespace std; std 36 src/io_paraview.cpp std::string prefix() const; std 38 src/io_paraview.cpp std::string filename; std 64 src/io_paraview.cpp catch(const std::exception&){ std 65 src/io_paraview.cpp std::throw_with_nested(std::logic_error std 69 src/io_paraview.cpp throw std::logic_error{"Unkown error in constructor of Paraview."}; std 91 src/io_paraview.cpp inline std::string DataOutputParameters::prefix() const { std 92 src/io_paraview.cpp std::ostringstream s; std 42 src/linHeat_algo.h std::clog << data << std::endl; std 56 src/linHeat_algo.h err_log << std::scientific; std 64 src/linHeat_algo.h << err_cal.l2_err() << ' ' << err_cal.h1_err() << std::endl; std 77 src/linHeat_algo.h << err_cal.l2_err() << ' ' << err_cal.h1_err() << std::endl; std 8 src/linHeat_main.cpp void print_errMsg(const std::exception&); std 15 src/linHeat_main.cpp std::cerr << "Dune error: " << e << std::endl; std 18 src/linHeat_main.cpp catch(const std::exception& e){ std 26 src/linHeat_main.cpp void print_errMsg(const std::exception& e){ std 27 src/linHeat_main.cpp std::cerr << e.what() << '\n'; std 29 src/linHeat_main.cpp std::rethrow_if_nested(e); std 31 src/linHeat_main.cpp catch(const std::exception& nested){ std 23 src/maxH_main.cpp using namespace std; std 27 src/secOrd_op_brusselator.cpp using namespace std; std 60 src/secOrd_op_brusselator.cpp catch(const std::exception&){ std 61 src/secOrd_op_brusselator.cpp std::throw_with_nested(std::logic_error std 65 src/secOrd_op_brusselator.cpp throw std::logic_error {"Unkown error in constructor of Brusselator."}; std 74 src/secOrd_op_brusselator.cpp throw_with_nested(std::logic_error {"Error in constructor of Brusselator."}); std 109 src/secOrd_op_brusselator.cpp throw std::logic_error std 120 src/secOrd_op_brusselator.cpp throw std::logic_error std 26 src/secOrd_op_brusselator_impl.cpp using namespace std; std 48 src/secOrd_op_brusselator_impl.cpp inline Range massMatrix_loc(const std::size_t pt, const Quadrature& q, std 54 src/secOrd_op_brusselator_impl.cpp inline Jacobian_range stiffnessMatrix_loc(const std::size_t pt, const Quadrature& q, std 60 src/secOrd_op_brusselator_impl.cpp inline Range quad_massMatrix_loc(const std::size_t pt, const Quadrature& q, std 73 src/secOrd_op_brusselator_impl.cpp inline Range jacobian_quadMass_loc(const std::size_t pt, const Quadrature& q, std 85 src/secOrd_op_brusselator_impl.cpp static std::size_t calculate_matrix_row_size(const FE_function&); std 93 src/secOrd_op_brusselator_impl.cpp std::vector<Range> vec_range {}; std 94 src/secOrd_op_brusselator_impl.cpp std::vector<Jacobian_range> vec_jacRange {}; std 128 src/secOrd_op_brusselator_impl.cpp throw std::logic_error {"Error in constructor of Data. " std 142 src/secOrd_op_brusselator_impl.cpp catch(const std::exception&){ std 143 src/secOrd_op_brusselator_impl.cpp std::throw_with_nested(std::logic_error {"Error in constructor of " std 147 src/secOrd_op_brusselator_impl.cpp throw std::logic_error {"Unkown error in constructor of Brusselator_op."}; std 169 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 186 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 218 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 241 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 263 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 273 src/secOrd_op_brusselator_impl.cpp for(std::size_t localCol = 0; localCol < base_set.size(); ++localCol){ std 297 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 306 src/secOrd_op_brusselator_impl.cpp for(std::size_t localCol = 0; localCol < base_set.size(); ++localCol){ std 318 src/secOrd_op_brusselator_impl.cpp std::size_t calculate_matrix_row_size(const FE_function& fef){ std 56 src/secOrd_op_brusselator_impl.h std::unique_ptr<Data> d_ptr; std 30 src/secOrd_op_identity.cpp using namespace std; std 23 src/secOrd_op_initData.cpp using namespace std; std 46 src/secOrd_op_initData.cpp using std::begin; std 47 src/secOrd_op_initData.cpp using std::end; std 105 src/secOrd_op_initData.cpp :d_ptr {std::make_unique<Data>(gt)} std 47 src/secOrd_op_initData_impl.cpp using namespace std; std 95 src/secOrd_op_initData_impl.cpp r = x * y * std::exp(-6. * t); std 103 src/secOrd_op_initData_impl.cpp r = y * z * std::exp(-6. * t); std 120 src/secOrd_op_initData_impl.cpp std::cout << print_configuration(p, type) << std::endl; std 125 src/secOrd_op_initData_impl.cpp :random_fun {std::bind(Random_dist {hom_value, hom_value + pertubation}, std 286 src/secOrd_op_initData_impl.cpp std::string Esfem::Impl::print_configuration(const Esfem::Io::Parameter& p, std 288 src/secOrd_op_initData_impl.cpp std::ostringstream oss; std 296 src/secOrd_op_initData_impl.cpp std::string Esfem::Impl::dof_filename(const Io::Parameter& p, const Growth type){ std 297 src/secOrd_op_initData_impl.cpp std::string rv {}; std 317 src/secOrd_op_initData_impl.cpp :eid_ptr {std::make_unique<Explicit_initial_data>(gt,type)} std 322 src/secOrd_op_initData_impl.cpp rid_ptr {std::make_unique<Random_initial_data>(p, type)} std 69 src/secOrd_op_initData_impl.h std::function<void(const Domain&,Range&)> fun_impl; std 96 src/secOrd_op_initData_impl.h using Random_dist = std::uniform_real_distribution<>; std 98 src/secOrd_op_initData_impl.h using Random_engine = std::default_random_engine; std 101 src/secOrd_op_initData_impl.h std::function<double()> random_fun; std 362 src/secOrd_op_initData_impl.h std::string print_configuration(const Esfem::Io::Parameter&, const Esfem::Growth); std 366 src/secOrd_op_initData_impl.h std::string dof_filename(const Io::Parameter&, const Growth); std 373 src/secOrd_op_initData_impl.h const std::string dof_io_filename {}; std 375 src/secOrd_op_initData_impl.h std::unique_ptr<Impl::Explicit_initial_data> eid_ptr; std 377 src/secOrd_op_initData_impl.h std::unique_ptr<Impl::Random_initial_data> rid_ptr; std 59 src/secOrd_op_initData_u.cpp catch(const std::exception&){ std 60 src/secOrd_op_initData_u.cpp std::throw_with_nested(std::logic_error{"Error in constructor of Init_data_u."}); std 63 src/secOrd_op_initData_u.cpp throw std::logic_error{"Unknown error in constructor of Init_data_u."}; std 70 src/secOrd_op_initData_u.cpp std::cerr << "~Init_data_u(): delete d_ptr.\n"; std 90 src/secOrd_op_initData_u.cpp r = std::exp(-6.*t)*x*y; std 59 src/secOrd_op_initData_w.cpp catch(const std::exception&){ std 60 src/secOrd_op_initData_w.cpp std::throw_with_nested(std::logic_error{"Error in constructor of Init_data_w."}); std 63 src/secOrd_op_initData_w.cpp throw std::logic_error{"Unknown error in constructor of Init_data_w."}; std 70 src/secOrd_op_initData_w.cpp std::cerr << "~Init_data_w(): delete d_ptr.\n"; std 90 src/secOrd_op_initData_w.cpp r = std::exp(-6.*t)*y*z; std 26 src/secOrd_op_linearHeat.cpp using namespace std; std 60 src/secOrd_op_linearHeat.cpp Range mass_matrix(const std::size_t pt, const Quadrature& q, std 67 src/secOrd_op_linearHeat.cpp Jacobian_range stiffness_matrix(const std::size_t pt, const Quadrature& q, std 88 src/secOrd_op_linearHeat.cpp catch(const std::exception&){ std 89 src/secOrd_op_linearHeat.cpp std::throw_with_nested(std::logic_error{"Error in constructor of Linear_heat."}); std 92 src/secOrd_op_linearHeat.cpp throw std::logic_error{"Unkown error in constructor of Linear_heat."}; std 150 src/secOrd_op_linearHeat.cpp throw std::runtime_error{"Error in Linear_heat::mass_matrix()"}; std 172 src/secOrd_op_linearHeat.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 192 src/secOrd_op_linearHeat.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 38 src/secOrd_op_rhs.cpp :d_ptr {std::make_unique<Data>(gt, type)} std 51 src/secOrd_op_rhs.cpp :d_ptr {std::make_unique<Data>(gt)} std 47 src/secOrd_op_rhs_impl.cpp using namespace std; std 95 src/secOrd_op_rhs_impl.cpp void Rhs_fun::dassert(const bool assertion, const std::string& msg){ std 67 src/secOrd_op_rhs_impl.h std::function<void(const Domain&,Range&)> fun_impl; std 69 src/secOrd_op_rhs_impl.h void dassert(const bool assertion, const std::string& msg); std 325 src/secOrd_op_rhs_impl.h for(std::size_t pt = 0; pt < quad.nop(); ++pt){ std 78 src/secOrd_op_rhs_u.cpp catch(std::exception&){ std 79 src/secOrd_op_rhs_u.cpp std::throw_with_nested(std::logic_error std 83 src/secOrd_op_rhs_u.cpp throw std::logic_error {"Unknown error in constructor of Grid_and_time."}; std 90 src/secOrd_op_rhs_u.cpp std::cerr << "~Rhs_u(): delete d_ptr.\n"; std 119 src/secOrd_op_rhs_u.cpp const double u = std::exp(-6. * t) * x * y; std 152 src/secOrd_op_rhs_u.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 80 src/secOrd_op_rhs_w.cpp catch(std::exception&){ std 81 src/secOrd_op_rhs_w.cpp std::throw_with_nested(std::logic_error std 85 src/secOrd_op_rhs_w.cpp throw std::logic_error {"Unknown error in constructor of Grid_and_time."}; std 92 src/secOrd_op_rhs_w.cpp std::cerr << "~Rhs_w(): delete d_ptr.\n"; std 145 src/secOrd_op_rhs_w.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ std 70 src/secOrd_op_solutionDriven.cpp :d_ptr {std::make_unique<Data>(p, g, u)} std 27 src/secOrd_op_solutionDriven_impl.cpp using std::size_t; std 173 src/secOrd_op_solutionDriven_impl.cpp std::vector<MCF_op::Domain<Vector_fef> > std 187 src/secOrd_op_solutionDriven_impl.cpp (const std::vector<MCF_op::Domain<Vector_fef> >& basis){ std 161 src/secOrd_op_solutionDriven_impl.h std::vector<MCF_op::Domain<MCF_op::Vector_fef> > std 170 src/secOrd_op_solutionDriven_impl.h (const std::vector<MCF_op::Domain<MCF_op::Vector_fef> >& basis); std 180 src/secOrd_op_solutionDriven_impl.h return std::sqrt(v[0] * v[0] + v[1] * v[1] + v[2] * v[2]); std 200 src/secOrd_op_solutionDriven_impl.h evaluate(const std::size_t pt, std 209 src/secOrd_op_solutionDriven_impl.h evaluate(const std::size_t pt, std 218 src/secOrd_op_solutionDriven_impl.h jacobian(const std::size_t pt, std 20 test/io_dgf.cpp void print_errMsg(const std::exception&); std 22 test/io_dgf.cpp std::string name(const FEFun& fef){ std 27 test/io_dgf.cpp using namespace std; std 92 test/io_dgf.cpp void print_errMsg(const std::exception& e){ std 93 test/io_dgf.cpp std::cerr << e.what() << '\n'; std 95 test/io_dgf.cpp std::rethrow_if_nested(e); std 97 test/io_dgf.cpp catch(const std::exception& nested){