40 |
import esys.finley |
import esys.finley |
41 |
from esys.escript.linearPDEs import LinearPDE |
from esys.escript.linearPDEs import LinearPDE |
42 |
import os |
import os |
43 |
|
import math |
44 |
|
|
45 |
class FinleyFilter(BenchmarkFilter): |
class FinleyFilter(BenchmarkFilter): |
46 |
""" |
""" |
280 |
u*=(x[i]-i) |
u*=(x[i]-i) |
281 |
|
|
282 |
gamma_rad=self.gamma/360.*8*math.atan(1.) |
gamma_rad=self.gamma/360.*8*math.atan(1.) |
283 |
cg=maths.cos(self.gamma_rad) |
cg=math.cos(self.gamma_rad) |
284 |
sg=maths.sin(self.gamma_rad) |
sg=math.sin(self.gamma_rad) |
285 |
C=kronecker(domain) |
C=kronecker(domain) |
286 |
C[0,0]=cg**2+self.c*sg**2 |
C[0,0]=cg**2+self.c*sg**2 |
287 |
C[1,0]=(self.c-1.)*cg*sg |
C[1,0]=(self.c-1.)*cg*sg |