Source_Generique_base.cpp

/****************************************************************************
* Copyright (c) 2026, CEA
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
* 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
* 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
* 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
* IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
* OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
*****************************************************************************/
 
 
#include <Source_Generique_base.h>
#include <Equation_base.h>
#include <Postraitement.h>
 
Implemente_base(Source_Generique_base,"Source_Generique_base",Source_base);
// XD source_generique source_base source_generique NO_BRACE to define a source term depending on some discrete fields
// XD_CONT of the problem and (or) analytic expression. It is expressed by the way of a generic field usually used for
// XD_CONT post-processing.
// XD attr champ champ_generique_base champ REQ the source field
 
/*! @brief DOES NOTHING - to override in derived classes.
 *
 *     Prints the source term to an output stream.
 *
 * @param (Sortie& os) the output stream
 * @return (Sortie&) the modified output stream
 */
Sortie& Source_Generique_base::printOn(Sortie& os) const
{
  return os;
}
 
/*! @brief DOES NOTHING - to override in derived classes.
 *
 *     Reads a source term from an input stream.
 *
 * @param (Entree& is) the input stream
 * @return (Entree&) the modified input stream
 */
Entree& Source_Generique_base::readOn(Entree& is)
{
  is >> ch_source_;
  return is;
}
 
void Source_Generique_base::associer_pb(const Probleme_base&)
{
 
}
 
DoubleTab& Source_Generique_base::calculer(DoubleTab& resu) const
{
  resu=0;
  return ajouter(resu);
}
 
void Source_Generique_base::completer()
{
  Source_base::completer();
 
  Postraitement post_bidon;
  post_bidon.associer_nom_et_pb_base("", equation().probleme());
  ch_source_->completer(post_bidon);
  Entity loc = ch_source_->get_localisation();
  if (((loc == Entity::ELEMENT) && (localisation_source() != "elem"))
      || ((loc == Entity::FACE) && (localisation_source() != "faces")))
    {
      Nom nom_loc;
      if (loc == Entity::ELEMENT)
        nom_loc = "elem";
      else if (loc == Entity::FACE)
        nom_loc = "faces";
      Cerr << "----------------------------------------------------------------" << finl;
      Cerr << "Error during TRUST calculation :" << finl;
      Cerr << "----------------------------------------------------------------" << finl;
      Cerr << "The generic source field associated to the " << que_suis_je() << " source term " << finl;
      Cerr << "of the equation " << equation().que_suis_je() << " will be evaluated at localisation " << nom_loc
           << finl;
      Cerr << "while the source term itself must be computed at location " << localisation_source() << finl;
      Cerr << "Please contact TRUST support." << finl;
      exit();
    }
}
 
void Source_Generique_base::mettre_a_jour(double temps)
{
 
}