62 |
# Set Boundary Mask / pdelib Template "q" Parameter Vector |
# Set Boundary Mask / pdelib Template "q" Parameter Vector |
63 |
bndryMask = Vector(value = 0, what = contfunc) |
bndryMask = Vector(value = 0, what = contfunc) |
64 |
for i in range(meshDim): |
for i in range(meshDim): |
65 |
bndryMask += (x[i].whereZero() + (x[i]-recDim[i]).whereZero()) \ |
bndryMask += (whereZero(x[i]) + whereZero(x[i]-recDim[i])) \ |
66 |
* ones((numDim,)) |
* ones((numDim,)) |
67 |
|
|
68 |
# Set True Solution / pdelib Template "r" Parameter Vector |
# Set True Solution / pdelib Template "r" Parameter Vector |
132 |
# Set Boundary Mask / pdelib Template "q" Parameter Vector |
# Set Boundary Mask / pdelib Template "q" Parameter Vector |
133 |
bndryMask = Scalar(value = 0, what = contfunc) |
bndryMask = Scalar(value = 0, what = contfunc) |
134 |
for i in range(meshDim): |
for i in range(meshDim): |
135 |
bndryMask += (x[i].whereZero() + (x[i]-recDim[i]).whereZero()) * 1.0 |
bndryMask += (whereZero(x[i]) + whereZero(x[i]-recDim[i])) * 1.0 |
136 |
|
|
137 |
# Set True Solution / pdelib Template "r" Parameter Vector |
# Set True Solution / pdelib Template "r" Parameter Vector |
138 |
u = Scalar(value = 0, what = contfunc) |
u = Scalar(value = 0, what = contfunc) |