Mercurial > repos > public > sbplib_julia

annotate src/SbpOperators/volumeops/inference_trouble.txt @ 959:e9752c1e92f8 feature/variable_derivatives

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix conflict from merge
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 14 Mar 2022 09:46:22 +0100
parents b41180efb6c2
children f1c2a4fa0ee1
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
919
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 Innan ändringarna på den här branchen:
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 begin
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 using Sbplib
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 using Sbplib.Grids
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 using Sbplib.SbpOperators
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 g = EquidistantGrid((10,10),(0.,0.), (1.,1.))
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 v = evalOn(g, (x,y)->x^2+y^2+1)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 H = inner_product(g, 1., [1/2])
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 end
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 # Not type stable
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 LazyTensors.apply(H.t1, H.t2*v, 1,2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15 @code_warntype LazyTensors.apply(H.t1, H.t2*v, 1,2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
16 @code_warntype LazyTensors.apply(H.t1.tm, view(H.t2*v,:,1), 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
17
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18 # Nedan är halvdåliga
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
19 @code_warntype LazyTensors.apply(H.t1.tm, view(H.t2*v,1,:), 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
20 @code_warntype LazyTensors.apply(H.t1.tm, view(v,1,:), 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
21 @code_warntype LazyTensors.apply(H.t1.tm, view(v,:,1), 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
22 @code_warntype LazyTensors.apply(H.t1.tm, v[:,1], 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
23
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
24
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
25
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
26
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
27
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
28
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
29
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32 begin
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
33 using Sbplib
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
34 using Sbplib.Grids
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
35 using Sbplib.SbpOperators
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
36 import Sbplib.SbpOperators: Stencil
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
37 using Sbplib.RegionIndices
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
38
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
39 g = EquidistantGrid(10,0., 1.)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
40 v = evalOn(g, (x)->x^2+1)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
41 H = inner_product(g, 1., [1/2])
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
42 V = SbpOperators.volume_operator(g, Stencil(1.,center=1), (Stencil(1/2,center=1),), SbpOperators.even,1)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
43 b = SbpOperators.boundary_operator(g, Stencil(1/2,center=1), CartesianBoundary{1,Lower}())
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
44 end
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
45
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
46 @code_warntype LazyTensors.apply(H, H*v, 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
47 @code_warntype LazyTensors.apply(V, V*v, 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
48 @code_warntype LazyTensors.apply(b, b*v, 2)
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
49
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
50
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
51
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
52 begin
b41180efb6c2 Start refactor to improve type stability for apply(op,::TensorApplication,...)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
53 end