ceed-operator.c (381e65939e85104561074440c4dd3dd99bd0efff) ceed-operator.c (b8bf0bca2754e3c4d98d2d0a42f776e4e8976ea4)
1// Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors.
2// All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
3//
4// SPDX-License-Identifier: BSD-2-Clause
5//
6// This file is part of CEED: http://github.com/ceed
7
8#include <ceed/ceed.h>

--- 129 unchanged lines hidden (view full) ---

138**/
139static int CeedOperatorFieldView(CeedOperatorField field,
140 CeedQFunctionField qf_field,
141 CeedInt field_number, bool sub, bool input,
142 FILE *stream) {
143 const char *pre = sub ? " " : "";
144 const char *in_out = input ? "Input" : "Output";
145
1// Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors.
2// All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
3//
4// SPDX-License-Identifier: BSD-2-Clause
5//
6// This file is part of CEED: http://github.com/ceed
7
8#include <ceed/ceed.h>

--- 129 unchanged lines hidden (view full) ---

138**/
139static int CeedOperatorFieldView(CeedOperatorField field,
140 CeedQFunctionField qf_field,
141 CeedInt field_number, bool sub, bool input,
142 FILE *stream) {
143 const char *pre = sub ? " " : "";
144 const char *in_out = input ? "Input" : "Output";
145
146 fprintf(stream, "%s %s Field [%d]:\n"
146 fprintf(stream, "%s %s field %d:\n"
147 "%s Name: \"%s\"\n",
148 pre, in_out, field_number, pre, qf_field->field_name);
149
150 if (field->basis == CEED_BASIS_COLLOCATED)
151 fprintf(stream, "%s Collocated basis\n", pre);
152
153 if (field->vec == CEED_VECTOR_ACTIVE)
154 fprintf(stream, "%s Active vector\n", pre);

--- 21 unchanged lines hidden (view full) ---

176 ierr = CeedOperatorGetNumElements(op, &num_elem); CeedChk(ierr);
177 ierr = CeedOperatorGetNumQuadraturePoints(op, &num_qpts); CeedChk(ierr);
178
179 CeedInt total_fields = 0;
180 ierr = CeedOperatorGetNumArgs(op, &total_fields); CeedChk(ierr);
181 fprintf(stream, "%s %d elements with %d quadrature points each\n",
182 pre, num_elem, num_qpts);
183
147 "%s Name: \"%s\"\n",
148 pre, in_out, field_number, pre, qf_field->field_name);
149
150 if (field->basis == CEED_BASIS_COLLOCATED)
151 fprintf(stream, "%s Collocated basis\n", pre);
152
153 if (field->vec == CEED_VECTOR_ACTIVE)
154 fprintf(stream, "%s Active vector\n", pre);

--- 21 unchanged lines hidden (view full) ---

176 ierr = CeedOperatorGetNumElements(op, &num_elem); CeedChk(ierr);
177 ierr = CeedOperatorGetNumQuadraturePoints(op, &num_qpts); CeedChk(ierr);
178
179 CeedInt total_fields = 0;
180 ierr = CeedOperatorGetNumArgs(op, &total_fields); CeedChk(ierr);
181 fprintf(stream, "%s %d elements with %d quadrature points each\n",
182 pre, num_elem, num_qpts);
183
184 fprintf(stream, "%s %d Field%s\n", pre, total_fields,
184 fprintf(stream, "%s %d field%s\n", pre, total_fields,
185 total_fields>1 ? "s" : "");
186
185 total_fields>1 ? "s" : "");
186
187 fprintf(stream, "%s %d Input Field%s:\n", pre, op->qf->num_input_fields,
187 fprintf(stream, "%s %d input field%s:\n", pre, op->qf->num_input_fields,
188 op->qf->num_input_fields>1 ? "s" : "");
189 for (CeedInt i=0; i<op->qf->num_input_fields; i++) {
190 ierr = CeedOperatorFieldView(op->input_fields[i], op->qf->input_fields[i],
191 i, sub, 1, stream); CeedChk(ierr);
192 }
193
188 op->qf->num_input_fields>1 ? "s" : "");
189 for (CeedInt i=0; i<op->qf->num_input_fields; i++) {
190 ierr = CeedOperatorFieldView(op->input_fields[i], op->qf->input_fields[i],
191 i, sub, 1, stream); CeedChk(ierr);
192 }
193
194 fprintf(stream, "%s %d Output Field%s:\n", pre, op->qf->num_output_fields,
194 fprintf(stream, "%s %d output field%s:\n", pre, op->qf->num_output_fields,
195 op->qf->num_output_fields>1 ? "s" : "");
196 for (CeedInt i=0; i<op->qf->num_output_fields; i++) {
197 ierr = CeedOperatorFieldView(op->output_fields[i], op->qf->output_fields[i],
198 i, sub, 0, stream); CeedChk(ierr);
199 }
200 return CEED_ERROR_SUCCESS;
201}
202

--- 979 unchanged lines hidden (view full) ---

1182**/
1183int CeedOperatorView(CeedOperator op, FILE *stream) {
1184 int ierr;
1185
1186 if (op->is_composite) {
1187 fprintf(stream, "Composite CeedOperator\n");
1188
1189 for (CeedInt i=0; i<op->num_suboperators; i++) {
195 op->qf->num_output_fields>1 ? "s" : "");
196 for (CeedInt i=0; i<op->qf->num_output_fields; i++) {
197 ierr = CeedOperatorFieldView(op->output_fields[i], op->qf->output_fields[i],
198 i, sub, 0, stream); CeedChk(ierr);
199 }
200 return CEED_ERROR_SUCCESS;
201}
202

--- 979 unchanged lines hidden (view full) ---

1182**/
1183int CeedOperatorView(CeedOperator op, FILE *stream) {
1184 int ierr;
1185
1186 if (op->is_composite) {
1187 fprintf(stream, "Composite CeedOperator\n");
1188
1189 for (CeedInt i=0; i<op->num_suboperators; i++) {
1190 fprintf(stream, " SubOperator [%d]:\n", i);
1190 fprintf(stream, " SubOperator %d:\n", i);
1191 ierr = CeedOperatorSingleView(op->sub_operators[i], 1, stream);
1192 CeedChk(ierr);
1193 }
1194 } else {
1195 fprintf(stream, "CeedOperator\n");
1196 ierr = CeedOperatorSingleView(op, 0, stream); CeedChk(ierr);
1197 }
1198 return CEED_ERROR_SUCCESS;

--- 484 unchanged lines hidden --- 		1191 ierr = CeedOperatorSingleView(op->sub_operators[i], 1, stream);
1192 CeedChk(ierr);
1193 }
1194 } else {
1195 fprintf(stream, "CeedOperator\n");
1196 ierr = CeedOperatorSingleView(op, 0, stream); CeedChk(ierr);
1197 }
1198 return CEED_ERROR_SUCCESS;

--- 484 unchanged lines hidden ---