annotate d2_4th.txt @ 129:1aaeb46ba5f4 cell_based_test

Improve efficiency of apply by the following: - Remove divisions in interior loop by storing and multiplying by the reciprocal of grid spacing instead. - Add @inline to apply(::Laplace - Remove initialization of w = 0 in apply(::Stencil) by manually unrolling first iteration of the loop.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 14 Feb 2019 16:25:22 +0100
parents 8a30ea842c8e
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
36
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
1 # D2 order 4
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
2
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
3 boundary_stencils
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
4 2 -5 4 -1 0 0
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
5 1 -2 1 0 0 0
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
6 -4/43 59/43 -110/43 59/43 -4/43 0
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
7 -1/49 0 59/49 -118/49 64/49 -4/49
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
8
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
9 inner_stencil
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
10 -1/12 4/3 -5/2 4/3 -1/12
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
11
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
12 e
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
13 1
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
14
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
15 d1
8a30ea842c8e add fourth order stencil
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
16 -11/6 3 -3/2 1/3