138 |
sp.initialize(f=F,fixed_u_mask=mask,eta=ETA) |
sp.initialize(f=F,fixed_u_mask=mask,eta=ETA) |
139 |
u0=(1-x[0])*x[0]*[0.,1.] |
u0=(1-x[0])*x[0]*[0.,1.] |
140 |
p0=Scalar(P1,ReducedSolution(self.domain)) |
p0=Scalar(P1,ReducedSolution(self.domain)) |
141 |
u,p=sp.solve(u0,p0,show_details=VERBOSE, verbose=VERBOSE,max_iter=100,solver="GMRES",iter_restart=12) |
u,p=sp.solve(u0,p0,show_details=VERBOSE, verbose=VERBOSE,max_iter=100,solver="GMRES",iter_restart=20) |
142 |
|
|
143 |
error_v0=Lsup(u[0]-u0[0]) |
error_v0=Lsup(u[0]-u0[0]) |
144 |
error_v1=Lsup(u[1]-u0[1])/0.25 |
error_v1=Lsup(u[1]-u0[1])/0.25 |
324 |
sp.initialize(f=F,fixed_u_mask=mask,eta=ETA) |
sp.initialize(f=F,fixed_u_mask=mask,eta=ETA) |
325 |
u0=(1-x[0])*x[0]*(1-x[1])*x[1]*[0.,0.,1.] |
u0=(1-x[0])*x[0]*(1-x[1])*x[1]*[0.,0.,1.] |
326 |
p0=Scalar(P1,ReducedSolution(self.domain)) |
p0=Scalar(P1,ReducedSolution(self.domain)) |
327 |
u,p=sp.solve(u0,p0,show_details=VERBOSE, verbose=VERBOSE,max_iter=100,solver="GMRES",iter_restart=14) |
u,p=sp.solve(u0,p0,show_details=VERBOSE, verbose=VERBOSE,max_iter=100,solver="GMRES",iter_restart=20) |
328 |
|
|
329 |
error_v0=Lsup(u[0]-u0[0]) |
error_v0=Lsup(u[0]-u0[0]) |
330 |
error_v1=Lsup(u[1]-u0[1]) |
error_v1=Lsup(u[1]-u0[1]) |