Implementation of changing lambdas and mus over subdomain interface

[larsson4]

void LinElastSolver::BuildDomainOperators()
{
   // SanityCheckOnCoeffs();
   as.SetSize(numSub);
 
   for (int m = 0; m < numSub; m++)
   {
      as[m] = new BilinearForm(fes[m]);
      as[m]->AddDomainIntegrator(new ElasticityIntegrator(*(lambda_c[m]), *(mu_c[m])));
 
      if (full_dg)
      {
         as[m]->AddInteriorFaceIntegrator(
             new DGElasticityIntegrator(*(lambda_c[m]), *(mu_c[m]), alpha, kappa));
      }
   }
 
   a_itf = new InterfaceForm(meshes, fes, topol_handler); 
   a_itf->AddIntefaceIntegrator(new InterfaceDGElasticityIntegrator(lambda_c[0], mu_c[0], alpha, kappa));
}

How should the last line be handled?