specfem::medium_physics::compute_cosserat_couple_stress

template<typename T, typename PointPropertiesType, typename PointStressType, typename PointAccelerationType>
void specfem::medium_physics::compute_cosserat_couple_stress(const PointPropertiesType &point_properties, const T factor, const PointStressType &point_stress, PointAccelerationType &acceleration)

Compute Cosserat couple stress contribution for micropolar elastic media.

Provides compile-time dispatch to medium-specific implementations. Adds the moment equilibrium contribution to acceleration based on stress tensor asymmetry. Only modifies the spin component for Cosserat media; for all other media this is a no-op.

Note

The caller is responsible for zero-initialising acceleration before calling this function. The contribution is accumulated with +=, consistent with the scatter pattern where the epilogue applies the global sign negation.

Template Parameters:

  • T – Scalar type for the weight product

  • PointPropertiesType – Point-wise material properties

  • PointStressType – Point-wise physical stress tensor

  • PointAccelerationType – Point-wise acceleration field

Parameters:

  • point_properties – Cosserat material properties

  • factor – Integration factor w(iz)*w(ix)*J for this GLL point

  • point_stress – Physical stress tensor at this GLL point

  • acceleration[in, out] – Acceleration delta (spin component modified)

Implementation Details