fields *should* be ok with this check

CEED · Nov 8, 2023 · 9f69eb1 · 9f69eb1
1 parent 6429868
commit 9f69eb1
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/interface/ceed-operator.c b/interface/ceed-operator.c
@@ -636,9 +636,10 @@ int CeedOperatorCreateAtPoints(Ceed ceed, CeedQFunction qf, CeedQFunction dqf, C
 
   CeedCall(CeedCalloc(1, op));
   CeedCall(CeedReferenceCopy(ceed, &(*op)->ceed));
-  (*op)->ref_count   = 1;
-  (*op)->input_size  = -1;
-  (*op)->output_size = -1;
+  (*op)->ref_count    = 1;
+  (*op)->is_at_points = true;
+  (*op)->input_size   = -1;
+  (*op)->output_size  = -1;
   CeedCall(CeedQFunctionReferenceCopy(qf, &(*op)->qf));
   if (dqf && dqf != CEED_QFUNCTION_NONE) CeedCall(CeedQFunctionReferenceCopy(dqf, &(*op)->dqf));
   if (dqfT && dqfT != CEED_QFUNCTION_NONE) CeedCall(CeedQFunctionReferenceCopy(dqfT, &(*op)->dqfT));
@@ -801,7 +802,7 @@ int CeedOperatorSetField(CeedOperator op, const char *field_name, CeedElemRestri
     op->has_restriction = true;  // Restriction set, but num_elem may be 0
   }
   CeedCall(CeedBasisReferenceCopy(b, &(*op_field)->basis));
-  if (op->num_qpts == 0) op->num_qpts = num_qpts;
+  if (op->num_qpts == 0 && !op->is_at_points) op->num_qpts = num_qpts;  // no consistent number of qpts for OperatorAtPoints
   op->num_fields += 1;
   CeedCall(CeedStringAllocCopy(field_name, (char **)&(*op_field)->field_name));
   return CEED_ERROR_SUCCESS;