nonlininteg_vecconvection_pa_diag.cpp

Go to the documentation of this file.

// Copyright (c) 2010-2026, Lawrence Livermore National Security, LLC. Produced
// at the Lawrence Livermore National Laboratory. All Rights reserved. See files
// LICENSE and NOTICE for details. LLNL-CODE-806117.
//
// This file is part of the MFEM library. For more information and source code
// availability visit https://mfem.org.
//
// MFEM is free software; you can redistribute it and/or modify it under the
// terms of the BSD-3 license. We welcome feedback and contributions, see file
// CONTRIBUTING.md for details.
 
#include "../ceed/interface/util.hpp"
#include "./nonlininteg_vecconvection_pa_diag.hpp" // IWYU pragma: keep
 
namespace mfem
{
 
void VectorConvectionNLFIntegrator::AssembleGradDiagonalPA(Vector &de) const
{
   MFEM_VERIFY(!DeviceCanUseCeed(),
               "VectorConvectionNLFIntegrator PA gradients are not supported "
               "with the libCEED backend");
 
   if (dim == 2)
   {
      GradDiagPA2D::Run(d1d, q1d, ne,
                        maps->B.Read(),
                        maps->G.Read(),
                        pa_adj.Read(),
                        pa_u.Read(),
                        de.ReadWrite(),
                        d1d, q1d);
   }
   else if (dim == 3)
   {
      GradDiagPA3D::Run(d1d, q1d, ne,
                        maps->B.Read(),
                        maps->G.Read(),
                        pa_adj.Read(),
                        pa_u.Read(),
                        de.ReadWrite(),
                        d1d, q1d);
   }
   else
   {
      MFEM_ABORT("Unsupported dimension");
   }
}
 
} // namespace mfem