1#ifndef PSCF_FIELD_COMPARISON_TPP
2#define PSCF_FIELD_COMPARISON_TPP
11#include "FieldComparison.h"
36 for (
int i = begin_; i < n; ++i) {
37 diff = std::abs(a[i] - b[i]);
38 if (std::isnan(diff)) {
44 }
else if (diff > maxDiff_) {
47 rmsDiff_ += diff*diff;
49 rmsDiff_ = rmsDiff_/double(n);
50 rmsDiff_ = sqrt(rmsDiff_);
67 for (i = 0; i < m; ++i) {
71 for (j = begin_; j < n; ++j) {
72 diff = std::abs(a[i][j] - b[i][j]);
73 if (std::isnan(diff)) {
79 }
else if (diff > maxDiff_) {
82 rmsDiff_ += diff*diff;
85 rmsDiff_ = rmsDiff_/double(m*n);
86 rmsDiff_ = sqrt(rmsDiff_);
FieldComparison(int begin=0)
Default constructor.
double compare(FT const &a, FT const &b)
Compare individual fields.
int capacity() const
Return allocated size.
Dynamically allocatable contiguous array template.
#define UTIL_CHECK(condition)
Assertion macro suitable for serial or parallel production code.
C++ namespace for polymer self-consistent field theory (PSCF).
Utility classes for scientific computation.