PLaSK library

Loading...
Searching...
No Matches
brent.hpp
Go to the documentation of this file.
1/*
2 * This file is part of PLaSK (https://plask.app) by Photonics Group at TUL
3 * Copyright (c) 2022 Lodz University of Technology
4 *
5 * This program is free software: you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation, version 3.
8 *
9 * This program is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU General Public License for more details.
13 */
14#ifndef PLASK__OPTICAL_EFFECTIVE_BRENT_H
15#define PLASK__OPTICAL_EFFECTIVE_BRENT_H
16
17#include "rootdigger.hpp"
18
19 namespace plask { namespace optical { namespace effective {
20
21 class RootBrent: public RootDigger {
22
23 // Write log message
24 template <typename... Args>
25 void writelog(LogLevel level, const std::string& msg, Args&&... args) const {
26 std::string prefix = solver.getId(); prefix += ": "; prefix += log_value.chartName(); prefix += ": ";
27 plask::writelog(level, prefix + msg, std::forward<Args>(args)...);
28 }
29
30 protected:
31
32 double axisBrent(dcomplex start, double& fx, bool real) const;
33
34 public:
35
36 // Constructor
37 RootBrent(Solver& solver, const function_type& val_fun, DataLog<dcomplex,dcomplex>& log_value,
38 const Params& pars): RootDigger(solver, val_fun, log_value, pars) {}
39
40
41 dcomplex find(dcomplex start) const override;
42};
43
44}}} // namespace plask::optical::effective
45#endif // PLASK__OPTICAL_EFFECTIVE_BRENT_H