本文整理汇总了Python中problem.Problem.solve方法的典型用法代码示例。如果您正苦于以下问题:Python Problem.solve方法的具体用法?Python Problem.solve怎么用?Python Problem.solve使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类problem.Problem
的用法示例。
在下文中一共展示了Problem.solve方法的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: testKnownProblem
# 需要导入模块: from problem import Problem [as 别名]
# 或者: from problem.Problem import solve [as 别名]
def testKnownProblem(self):
def f(x,y):
return -2*x*y*(2*exp(2*x)*(x**2+y**2)+2*exp(2*x)*x)-(1+x**2*y)*(4*exp(2*x)*(x**2+y**2)+8*exp(2*x)*x+2*exp(2*x))-2*x**2*exp(2*x)*y-2*(1+x**2*y)*exp(2*x)
def k(x,y):
return 1+x**2*y
def u(x,y):
return exp(2*x)*(x**2+y**2)
def ux(x,y):
return 2*exp(2*x)*(x**2+y**2)+2*exp(2*x)*x
def uy(x,y):
return 2*exp(2*x)*y
ewidth = 30
eheight = 30
xwidth = 2.
ywidth = 2.
ll = (-1.,-1.)
[nodes, boundary_nodes, tris] = generateRectangularMesh((ewidth, eheight), ll, (xwidth/(ewidth-1),ywidth/(eheight-1)))
p = Problem(nodes, boundary_nodes, tris)
# U, g = p.solve(lambda x,y: cos(x*y*pi) - 0.5, lambda x,y: 1, lambda x,y: 0, ux)
Uest, g = p.solve(f, k, u, ux)
Utrue = [u(n.pos[0],n.pos[1]) for n in p.free_nodes]
# view_rectangular_grid_solution((ewidth, eheight), p, Uest, g)
# view_rectangular_grid_solution((ewidth, eheight), p, Utrue, g)
error = np.sqrt(np.sum((Utrue - Uest) ** 2))
print error
self.assertTrue(error < 0.5)
示例2: main
# 需要导入模块: from problem import Problem [as 别名]
# 或者: from problem.Problem import solve [as 别名]
def main():
args = argparser().parse_args()
p = Problem(args.input.read(), debug=args.debug)
p.solve()
args.output.write(str(p))