https://gitweb.dealii.org/

Fix a parenthesis.

I see a bunch of warnings of the form

/Users/runner/work/dealii/dealii/include/deal.II/matrix_free/tools.h:497:35:
error: taking the absolute value of unsigned type 'bool' has no effect [-Werror,-Wabsolute-value]
854
std::abs(values_dofs[j][0] &&

on the CI tester.

commit | commitdiff | tree

David Wells [Thu, 29 Jul 2021 16:21:54 +0000 (12:21 -0400)]

Add a lot more Witherden-Vincent rules.

This includes most of the even-degree rules.

For posterity: here's the script I used to print quadpy rules rather than
copying numbers individually:

import quadpy as qp

def format_16(x):
    return '{:.16e}'.format(float(x))

def print_scheme(scheme):
    sym_data = scheme.symmetry_data
    print("// WV-{}, {}D".format(scheme.degree, scheme.dim))
    for key in sym_data.keys():
        if key == 'centroid' or key == 's4':
            print("b_point_permutations.push_back({centroid});")
            print("b_weights.push_back({});".format(format_16(sym_data[key][0][0])))
        elif key == 'd3_aa' or key == 's31':
            weights = [format_16(w) for w in sym_data[key][0]]
            points = [format_16(p) for p in sym_data[key][1]]
            assert len(weights) == len(points)
            for i in range(len(weights)):
                print("process_point_1({}, {});".format(points[i], weights[i]))
        elif key == 'd3_ab' or key == 's211':
            weights = [format_16(w) for w in sym_data[key][0]]
            points1 = [format_16(p) for p in sym_data[key][1]]
            points2 = [format_16(p) for p in sym_data[key][2]]
            assert len(weights) == len(points1)
            assert len(points2) == len(points1)
            for i in range(len(weights)):
                print(("process_point_3({},"
                       "\n                {},"
                       "\n                {});")
                      .format(points1[i], points2[i], weights[i]))
        elif key == 's22':
            assert scheme.dim == 3
            weights = [format_16(w) for w in sym_data[key][0]]
            points = [format_16(p) for p in sym_data[key][1]]
            assert len(weights) == len(points)
            for i in range(len(weights)):
                print(("process_point_2({},"
                       "\n                {});")
                      .format(points[i], weights[i]))
        else:
            assert False

commit | commitdiff | tree

Martin Kronbichler [Sun, 1 Aug 2021 11:40:40 +0000 (13:40 +0200)]

Merge pull request #12613 from drwells/fix-structdim-comments

Fix two comments in TriaAccessor.

commit | commitdiff | tree

Marc Fehling [Sat, 31 Jul 2021 21:06:00 +0000 (15:06 -0600)]

Merge pull request #12559 from tamiko/modernize_step_69

Step-69: Use SolutionTransfer instead of low-level boost archive

commit | commitdiff | tree

Matthias Maier [Sat, 31 Jul 2021 14:25:52 +0000 (09:25 -0500)]

Step-69: Apply review suggestion

commit | commitdiff | tree

Matthias Maier [Sat, 31 Jul 2021 14:23:37 +0000 (09:23 -0500)]

Step-69: Use SolutionTransfer instead of low-level boost archive