Make Error Call to Registry

[TLWise]

Check these boxes if you have followed the posting rules.

• Q&A General is the most appropriate section for my question
• I have consulted the posting Guidelines on the Discussions front page
• I have searched the Discussions forum and my question has not been asked before
• I have searched the MOOSE website and the documentation does not answer my question
• I have formatted my post following the posting guidelines (screenshots as a last resort, triple back quotes around pasted text)

Question

How might I resolve this particular error?:

See complete error detailed below in comments.

[GiudGiud]

Hello
We ll need more context.
what were you trying to do? Which line is causing this?

[TLWise]

This is the rough code I started with prior to making modifications.:
Code.h

#pragma once

#include "SideIntegralPostprocessor.h"

class ADDragForcePostprocessor : public SideIntegralPostprocessor
{
public:
  static InputParameters validParams();
  ADDragForcePostprocessor(const InputParameters & params);

protected:
  virtual ADReal adComputeQpIntegral() override;

  // Coupled fields
  const ADVariableValue & _p;
  const ADVariableGradient & _grad_u;
  const ADVariableGradient & _grad_v;
  const ADVariableGradient * _grad_w; // optional (nullptr in 2D)

  // Viscosity (material property or constant)
  const ADMaterialProperty<Real> * _mu_prop;
  const Real _mu_param;

  // Options
  const RealVectorValue _drag_dir_raw; // user input, will be normalized
  RealVectorValue _drag_dir_unit;
  const bool _include_pressure;
  const bool _force_on_solid;
};

Code.C
#include "ADDragForcePostprocessor.h"
#include "MooseMesh.h"
#include "RankTwoTensor.h"

registerMooseObject("MooseApp", ADDragForcePostprocessor);

InputParameters
ADDragForcePostprocessor::validParams()
{
  InputParameters params = SideIntegralPostprocessor::validParams();
  params.addClassDescription(
      "Integrates the drag force F_d = ∫_S (σ·n)·e_d dS with "
      "σ = -p I + μ(grad u + grad u^T).");

  params.addRequiredCoupledVar("u", "x-velocity component");
  params.addRequiredCoupledVar("v", "y-velocity component");
  params.addCoupledVar("w", "z-velocity component (optional; omit in 2D)");
  params.addRequiredCoupledVar("p", "pressure");

  params.addParam<MaterialPropertyName>(
      "mu_property", "",
      "AD material property name for dynamic viscosity μ [Pa·s]. "
      "If empty, the 'mu' parameter is used.");

  params.addParam<Real>(
      "mu", -1.0,
      "Constant dynamic viscosity μ [Pa·s]; ignored if 'mu_property' is provided.");

  params.addParam<RealVectorValue>(
      "drag_direction", RealVectorValue(1.0, 0.0, 0.0),
      "Direction for the drag component (need not be unit).");

  params.addParam<bool>("include_pressure", true,
                        "Include the pressure contribution (-p I) in σ.");
  params.addParam<bool>("force_on_solid", true,
                        "Return the force ON THE SOLID (i.e., minus the fluid traction).");

  return params;
}

ADDragForcePostprocessor::ADDragForcePostprocessor(const InputParameters & params)
  : SideIntegralPostprocessor(params),
    _p(adCoupledValue("p")),
    _grad_u(adCoupledGradient("u")),
    _grad_v(adCoupledGradient("v")),
    _grad_w(isCoupled("w") ? &adCoupledGradient("w") : nullptr),
    _mu_prop(nullptr),
    _mu_param(getParam<Real>("mu")),
    _drag_dir_raw(getParam<RealVectorValue>("drag_direction")),
    _include_pressure(getParam<bool>("include_pressure")),
    _force_on_solid(getParam<bool>("force_on_solid"))
{
  const auto mu_name = getParam<MaterialPropertyName>("mu_property");
  if (!mu_name.empty())
    _mu_prop = &getADMaterialProperty<Real>(mu_name);

  if (_drag_dir_raw.norm() == 0.0)
    mooseError("drag_direction must be non-zero");
  _drag_dir_unit = _drag_dir_raw;
  _drag_dir_unit /= _drag_dir_raw.norm();

  if (!_mu_prop && _mu_param < 0.0)
    mooseError("Provide either a nonnegative 'mu' or a valid 'mu_property'.");
}

ADReal
ADDragForcePostprocessor::adComputeQpIntegral()
{
  const unsigned int dim = _mesh.dimension();

  // Build grad(u) as a rank-2 tensor with rows = component gradients
  ADRankTwoTensor gradV;
  gradV.zero();

  // Row for u-component
  for (unsigned int j = 0; j < dim; ++j)
    gradV(0, j) = _grad_u[_qp](j);

  // Row for v-component
  for (unsigned int j = 0; j < dim; ++j)
    gradV(1, j) = _grad_v[_qp](j);

  // Row for w-component if present (dim==3)
  if (_grad_w && dim == 3)
    for (unsigned int j = 0; j < dim; ++j)
      gradV(2, j) = (*_grad_w)[_qp](j);

  // Symmetric gradient and Cauchy stress
  const ADReal mu = _mu_prop ? (**_mu_prop)[_qp] : _mu_param;
  ADRankTwoTensor sigma = mu * (gradV + gradV.transpose());

  if (_include_pressure)
  {
    ADRankTwoTensor I;
    I.zero();
    I(0,0) = 1.0;
    I(1,1) = 1.0;
    if (dim == 3) I(2,2) = 1.0;
    sigma -= _p[_qp] * I;
  }

  // Traction and streamwise component
  const ADRealVectorValue n(_normals[_qp]);              // outward normal of the FLUID domain
  const ADRealVectorValue t = sigma * n;                 // traction on the fluid
  ADReal integrand = t * ADRealVectorValue(_drag_dir_unit);

  // Convert to force on the SOLID if requested (Newton's third law)
  if (_force_on_solid)
    integrand = -integrand;

  return integrand; // SideIntegralPostprocessor handles JxW*coord transform
}

[loganharbour]

Can you provide the rest of the error?

[TLWise]

@loganharbour, Please find the complete error below, thanks for taking a look:

(moose) twise@twise-GS73VR-6RF:~/projects/crab$ make -j4
Compiling C++ (in opt mode) /home/twise/projects/crab/build/unity_src/postprocessors_Unity.C...
In file included from /home/twise/projects/crab/src/postprocessors/ADDragForcePostprocessor.C:1,
from /home/twise/projects/crab/build/unity_src/postprocessors_Unity.C:2:
/home/twise/projects/crab/build/header_symlinks/ADDragForcePostprocessor.h:12:18: error: 'virtual ADReal ADDragForcePostprocessor::adComputeQpIntegral()' marked 'override', but does not override
12 | virtual ADReal adComputeQpIntegral() override;
| ^~~~~~~~~~~~~~~~~~~
/home/twise/projects/crab/src/postprocessors/ADDragForcePostprocessor.C: In member function 'virtual ADReal ADDragForcePostprocessor::adComputeQpIntegral()':
/home/twise/projects/crab/src/postprocessors/ADDragForcePostprocessor.C:89:33: error: no match for 'operator*' (operand type is 'const ADMaterialProperty')
89 | const ADReal mu = _mu_prop ? (_mu_prop)[_qp] : _mu_param;
| ^~~~~~~~~~
In file included from /home/twise/projects/moose/framework/build/header_symlinks/DataIO.h:18,
from /home/twise/projects/moose/framework/build/header_symlinks/RestartableData.h:13,
from /home/twise/projects/moose/framework/build/header_symlinks/ReporterData.h:12,
from /home/twise/projects/moose/framework/build/header_symlinks/ReporterInterface.h:12,
from /home/twise/projects/moose/framework/build/header_symlinks/UserObject.h:18,
from /home/twise/projects/moose/framework/build/header_symlinks/SideUserObject.h:13,
from /home/twise/projects/moose/framework/build/header_symlinks/SidePostprocessor.h:12,
from /home/twise/projects/moose/framework/build/header_symlinks/SideIntegralPostprocessor.h:12,
from /home/twise/projects/crab/build/header_symlinks/ADDragForcePostprocessor.h:3:
/home/twise/projects/moose/framework/build/header_symlinks/RankThreeTensor.h:293:1: note: candidate: 'template RankThreeTensorTempl operator*(T, const RankThreeTensorTempl&)'
293 | operator*(T a, const RankThreeTensorTempl & b)
| ^~~~~~~~
/home/twise/projects/moose/framework/build/header_symlinks/RankThreeTensor.h:293:1: note: template argument deduction/substitution failed:
/home/twise/projects/crab/src/postprocessors/ADDragForcePostprocessor.C:89:35: note: candidate expects 2 arguments, 1 provided
89 | const ADReal mu = _mu_prop ? (_mu_prop)[_qp] : _mu_param;
| ^~~~~~~~
/home/twise/projects/moose/framework/build/header_symlinks/RankThreeTensor.h:301:1: note: candidate: 'template RankTwoTensorTempl operator*(const libMesh::VectorValue&, const RankThreeTensorTempl&)'
301 | operator*(const libMesh::VectorValue & p, const RankThreeTensorTempl & b)
| ^~~~~~~~
/home/twise/projects/moose/framework/build/header_symlinks/RankThreeTensor.h:301:1: note: template argument deduction/substitution failed:
/home/twise/projects/crab/src/postprocessors/ADDragForcePostprocessor.C:89:35: note: 'const ADMaterialProperty' is not derived from 'const libMesh::VectorValue'
89 | const ADReal mu = _mu_prop ? (_mu_prop)[_qp] : _mu_param;
| ^~~~~~~~
In file included from /home/twise/projects/moose/framework/build/header_symlinks/DataIO.h:19:
/home/twise/projects/moose/framework/build/header_symlinks/RankFourTensor.h:562:1: note: candidate: 'template<class T1, class T2> typename std::enable_if<libMesh::ScalarTraits::value, RankFourTensorTempl<decltype ((T1() * T2()))> >::type operator*(const T1&, const RankFourTensorTempl&)'
562 | operator*(const T1 & a, const RankFourTensorTempl & b) ->
| ^~~~~~~~
/home/twise/projects/moose/framework/build/header_symlinks/RankFourTensor.h:562:1: note: template argument deduction/substitution failed:
/home/twise/projects/crab/src/postprocessors/ADDragForcePostprocessor.C:89:35: note: candidate expects 2 arguments, 1 provided
89 | const ADReal mu = _mu_prop ? (_mu_prop)[_qp] : _mu_param;
| ^~~~~~~~
In file included from /home/twise/miniforge/envs/moose/petsc/include/petscsystypes.h:584,
from /home/twise/miniforge/envs/moose/petsc/include/petscsys.h:44,
from /home/twise/miniforge/envs/moose/petsc/include/petscbag.h:3,
from /home/twise/miniforge/envs/moose/petsc/include/petsc.h:6,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/petsc_macro.h:49,
from /home/twise/projects/moose/framework/build/header_symlinks/MooseTypes.h:30,
from /home/twise/projects/moose/framework/build/header_symlinks/MooseUtils.h:13,
from /home/twise/projects/moose/framework/build/header_symlinks/InputParameters.h:13,
from /home/twise/projects/moose/framework/build/header_symlinks/DistributionInterface.h:12,
from /home/twise/projects/moose/framework/build/header_symlinks/UserObject.h:13:
/home/twise/miniforge/envs/moose/petsc/include/petsccxxcomplexfix.h:100:1: note: candidate: 'PetscComplex operator*(const PetscComplex&, const PetscInt&)'
100 | PETSC_CXX_COMPLEX_FIX(PetscInt)
| ^~~~~~~~~~~~~~~~~~~~~
/home/twise/miniforge/envs/moose/petsc/include/petsccxxcomplexfix.h:100:1: note: candidate expects 2 arguments, 1 provided
/home/twise/miniforge/envs/moose/petsc/include/petsccxxcomplexfix.h:100:1: note: candidate: 'PetscComplex operator*(const PetscInt&, const PetscComplex&)'
100 | PETSC_CXX_COMPLEX_FIX(PetscInt)
| ^~~~~~~~~~~~~~~~~~~~~
/home/twise/miniforge/envs/moose/petsc/include/petsccxxcomplexfix.h:100:1: note: candidate expects 2 arguments, 1 provided
In file included from /home/twise/miniforge/envs/moose/lib/gcc/x86_64-conda-linux-gnu/12.3.0/include/c++/bits/std_thread.h:43,
from /home/twise/miniforge/envs/moose/lib/gcc/x86_64-conda-linux-gnu/12.3.0/include/c++/thread:43,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/threads_pthread.h:31,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/threads.h:126,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/reference_counter.h:25,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/reference_counted_object.h:24,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/dof_object.h:29,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/elem.h:26,
from /home/twise/miniforge/envs/moose/libmesh/include/libmesh/remote_elem.h:24,
from /home/twise/projects/moose/framework/build/header_symlinks/libMeshReducedNamespace.h:17,
from /home/twise/projects/moose/framework/build/header_symlinks/Moose.h:12,
from /home/twise/projects/moose/framework/build/header_symlinks/MooseTypes.h:12:
/home/twise/miniforge/envs/moose/lib/gcc/x86_64-conda-linux-gnu/12.3.0/include/c++/bits/unique_ptr.h: In instantiation of 'std::__detail::__unique_ptr_t<_Tp> std::make_unique(_Args&& ...) [with _Tp = ADDragForcePostprocessor; _Args = {const InputParameters&}; __detail::__unique_ptr_t<_Tp> = __detail::__unique_ptr_t]':
/home/twise/projects/moose/framework/build/header_symlinks/Registry.h:337:31: required from 'std::unique_ptr RegistryEntry::build(const InputParameters&) [with T = ADDragForcePostprocessor]'
/home/twise/projects/moose/framework/build/header_symlinks/Registry.h:334:1: required from here
/home/twise/miniforge/envs/moose/lib/gcc/x86_64-conda-linux-gnu/12.3.0/include/c++/bits/unique_ptr.h:1065:30: error: invalid new-expression of abstract class type 'ADDragForcePostprocessor'
1065 | { return unique_ptr<_Tp>(new _Tp(std::forward<_Args>(__args)...)); }
| ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/home/twise/projects/crab/build/header_symlinks/ADDragForcePostprocessor.h:5:7: note: because the following virtual functions are pure within 'ADDragForcePostprocessor':
5 | class ADDragForcePostprocessor : public SideIntegralPostprocessor
| ^~~~~~~~~~~~~~~~~~~~~~~~
/home/twise/projects/moose/framework/build/header_symlinks/SideIntegralPostprocessor.h:35:16: note: 'virtual libMesh::Real SideIntegralPostprocessor::computeQpIntegral()'
35 | virtual Real computeQpIntegral() = 0;
| ^~~~~~~~~~~~~~~~~
make: *** [/home/twise/projects/moose/framework/build.mk:151: /home/twise/projects/crab/build/unity_src/postprocessors_Unity.x86_64-conda-linux-gnu.opt.lo] Error 1

[loganharbour]

This is the significant error:

In file included from /home/twise/projects/crab/src/postprocessors/ADDragForcePostprocessor.C:1,
from /home/twise/projects/crab/build/unity_src/postprocessors_Unity.C:2:
/home/twise/projects/crab/build/header_symlinks/ADDragForcePostprocessor.h:12:18: error: 'virtual ADReal ADDragForcePostprocessor::adComputeQpIntegral()' marked 'override', but does not override
12 | virtual ADReal adComputeQpIntegral() override;
| ^~~~~~~~~~~~~~~~~~~

Postprocessors don't produce AD values. They can consume AD variables, but no derivatives are produced. You should be overriding the method computeQpIntegral instead, returning a Real value.

Secondly, this part:

const ADReal mu = _mu_prop ? (**_mu_prop)[_qp] : _mu_param;

If you want to get the value from _mu_prop (which is a pointer), you would want to access it with:

(*_mu_prop)[_qp]

[TLWise]

@loganharbour, This is my remaining error after incorporating the modifications noted.

make -j4
Checking if header needs updating: /home/twise/projects/moose/framework/include/base/MooseRevision.h...
Checking if header needs updating: /home/twise/projects/moose/modules/combined/include/base/CombinedRevision.h...
Creating Unity /home/twise/projects/moose/modules/combined/build/unity_src/postprocessors_Unity.C
Compiling C++ (in opt mode) /home/twise/projects/moose/modules/combined/build/unity_src/postprocessors_Unity.C...
In file included from /home/twise/projects/moose/modules/combined/build/unity_src/postprocessors_Unity.C:2:
/home/twise/projects/moose/modules/combined/src/postprocessors/DragForcePostprocessor.C: In member function 'virtual libMesh::Real DragForcePostprocessor::computeQpIntegral()':
/home/twise/projects/moose/modules/combined/src/postprocessors/DragForcePostprocessor.C:111:10: error: cannot convert 'ADReal' {aka 'MetaPhysicL::DualNumber<double, MetaPhysicL::SemiDynamicSparseNumberArray<double, long unsigned int, MetaPhysicL::NWrapper<64> >, true>'} to 'libMesh::Real' {aka 'double'} in return
  111 |   return integrand; // SideIntegralPostprocessor handles JxW*coord transform
      |          ^~~~~~~~~
      |          |
      |          ADReal {aka MetaPhysicL::DualNumber<double, MetaPhysicL::SemiDynamicSparseNumberArray<double, long unsigned int, MetaPhysicL::NWrapper<64> >, true>}
/home/twise/projects/moose/modules/combined/src/postprocessors/DragForcePostprocessor.C:112:2: error: expected '}' at end of input
  112 | }
      |  ^
/home/twise/projects/moose/modules/combined/src/postprocessors/DragForcePostprocessor.C:67:1: note: to match this '{'
   67 | {
      | ^
make: *** [/home/twise/projects/moose/framework/build.mk:151: /home/twise/projects/moose/modules/combined/build/unity_src/postprocessors_Unity.x86_64-conda-linux-gnu.opt.lo] Error 1

[GiudGiud]

Hello
.
You cannot return an ADReal for a Real return type.
You ll need to return return integrand.value(); or use metaphysl's raw_value routine

/home/twise/projects/moose/modules/combined/src/postprocessors/DragForcePostprocessor.C:112:2: error: expected '}' at end of input

These two last errors are unbalanced brackets. You can just count the brackets to figure that one out

[TLWise]

I figured the problem was the declaration as an ADReal, I just wanted to make sure my assumption was confirmed by the experts. Thank you.