libMesh/libmesh: coverage diff

Base bc1d1a Head #4503 40b389
Total Total +/- New
Rate 65.66% 65.68% +0.01% 100.00%
Hits 78815 78832 +17 10
Misses 41214 41198 -16 0
Filename Stmts Miss Cover
src/base/dof_map.C 0 -6 +0.47%
src/mesh/mesh_triangle_holes.C 0 +1 -0.27%
src/systems/equation_systems.C 0 -6 +0.74%
src/systems/system_projection.C +1 -5 +0.63%
TOTAL +1 -16 +0.01%
code
coverage unchanged
code
coverage increased
code
coverage decreased
+
line added or modified

src/base/dof_map.C

2621  
2622  
2623  
2624  
2625  
2626  
2627  

#ifdef LIBMESH_ENABLE_AMR
  // If we're asking for old dofs then we'd better have some
  if (old_dofs)
    libmesh_assert_greater_equal(n_old_dofs(), n_SCALAR_dofs());

  dof_id_type my_idx = old_dofs ?
    this->_first_old_scalar_df[vn] : this->_first_scalar_df[vn];

2750  
2751  
2752  
2753  
2754  
2755  
2756  
2757  
2758  
2759  
2760  
2761  
2762  

          if ((vn == v) || (vn == libMesh::invalid_uint))
            {
              if (var.type().family == SCALAR &&
                  (!elem ||
                   var.active_on_subdomain(elem->subdomain_id())))
                {
                  // We asked for this variable, so add it to the vector.
                  std::vector<dof_id_type> di_new;
                  this->SCALAR_dof_indices(di_new,v,true);
                  di.insert( di.end(), di_new.begin(), di_new.end());
                }
              else
                if (var.active_on_subdomain(elem->subdomain_id()))

src/mesh/mesh_triangle_holes.C

270  
271  
272  
273  
274  
275  
276  

    {
      ray_target = inside - Point(1);
      intersection_distances =
        this->find_ray_intersections(inside, ray_target);
    }

  // I'd make this an assert, but I'm not 100% confident we can't

src/systems/equation_systems.C

217  
218  
219  
220  
221  
222  
223  
224  
225  
226  
227  
228  
229  
230  
231  
232  
233  
234  

      // if necessary
      if (mesh_refine.coarsen_elements())
        {
          for (auto i : make_range(this->n_systems()))
            {
              System & sys = this->get_system(i);
              sys.get_dof_map().distribute_dofs(_mesh);
              sys.reinit_constraints();

              // Even if there weren't any constraint changes,
              // reinit_constraints() did prepare_send_list() for us.

              sys.restrict_vectors();
            }
          mesh_changed = true;
        }

      // Once vectors are all restricted, we can delete

src/systems/system_projection.C

385  
386  
387  
388 +
389  
390  
391 +
392 +
393 +
394  
395  
396  

          std::iota(vars.begin(), vars.end(), 0);
        }

      std::vector<unsigned int> regular_vars, vector_vars, scalar_vars;
      for (auto var : vars)
      {
        if (this->variable(var).type().family == SCALAR)
          scalar_vars.push_back(var);
        else if (FEInterface::field_type(this->variable_type(var)) == TYPE_SCALAR)
          regular_vars.push_back(var);
        else
          vector_vars.push_back(var);

435  
436  
437  
438 +
439  
440  
441 +
442 +
443 +
444 +
445 +
446  
447  
448  

      if (this->processor_id() == (this->n_processors()-1))
        {
          const DofMap & dof_map = this->get_dof_map();
          for (auto var : scalar_vars)
            {
              // We can just map SCALAR dofs directly across
              std::vector<dof_id_type> new_SCALAR_indices, old_SCALAR_indices;
              dof_map.SCALAR_dof_indices (new_SCALAR_indices, var, false);
              dof_map.SCALAR_dof_indices (old_SCALAR_indices, var, true);
              for (auto i : index_range(new_SCALAR_indices))
                new_vector.set(new_SCALAR_indices[i], old_vector(old_SCALAR_indices[i]));
            }
        }
    }