beta.hpp

Go to the documentation of this file.

/*
 * This file is part of PLaSK (https://plask.app) by Photonics Group at TUL
 * Copyright (c) 2022 Lodz University of Technology
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, version 3.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU General Public License for more details.
 */
#ifndef PLASK__MODULE_ELECTRICAL_BETA2D_H
#define PLASK__MODULE_ELECTRICAL_BETA2D_H
 
#include "electr2d.hpp"
#include "electr3d.hpp"
 
namespace plask { namespace electrical { namespace shockley {
 
template <typename GeometryT>
struct PLASK_SOLVER_API BetaSolver : public std::conditional<std::is_same<GeometryT, Geometry3D>::value,
                                                             ElectricalFem3DSolver,
                                                             ElectricalFem2DSolver<GeometryT>>::type {
    typedef typename std::conditional<std::is_same<GeometryT, Geometry3D>::value,
                                      ElectricalFem3DSolver,
                                      ElectricalFem2DSolver<GeometryT>>::type BaseClass;
 
  protected:
    std::vector<double> js;    
    std::vector<double> beta;  
 
    Tensor2<double> activeCond(size_t n, double PLASK_UNUSED(U), double jy, double PLASK_UNUSED(T)) override {
        jy = abs(jy);
        return Tensor2<double>(0., 10. * jy * this->active[n].height * getBeta(n) / log(1e7 * jy / getJs(n) + 1.));
    }
 
  public:
    double getBeta(size_t n) const {
        if (beta.size() <= n) throw Exception("{0}: no beta given for junction {1}", this->getId(), n);
        return beta[n];
    }
    void setBeta(size_t n, double beta) {
        if (this->beta.size() <= n) {
            this->beta.reserve(n + 1);
            for (size_t s = this->beta.size(); s <= n; ++s) this->beta.push_back(NAN);
        }
        this->beta[n] = beta;
        this->invalidate();
    }
 
    double getJs(size_t n) const {
        if (js.size() <= n) throw Exception("{0}: no js given for junction {1}", this->getId(), n);
        return js[n];
    }
    void setJs(size_t n, double js) {
        if (this->js.size() <= n) {
            this->js.reserve(n + 1);
            for (size_t s = this->js.size(); s <= n; ++s) this->js.push_back(1.);
        }
        this->js[n] = js;
        this->invalidate();
    }
 
    void loadConfiguration(XMLReader& source, Manager& manager) override;
 
    BetaSolver(const std::string& name = "");
 
    std::string getClassName() const override;
 
    ~BetaSolver();
};
 
}}}  // namespace plask::electrical::shockley
 
#endif