--- trunk/doc/examples/heatedblock.py 2008/10/17 00:15:17 1891 +++ trunk/doc/examples/heatedblock.py 2008/10/17 00:55:59 1892 @@ -50,8 +50,9 @@ +whereZero(x[2])*[0.,0.,1.] sigma0=(lam+2./3.*mu)*alpha*(T-T_ref)*kronecker(mydomain) mypde.setValue(A=C,X=sigma0,q=msk) +mypde.setSolverMethod(mypde.PCG,mypde.AMG) #... solve pde ... -u=mypde.getSolution() +u=mypde.getSolution(verbose=True) #... calculate von-Misses g=grad(u) sigma=mu*(g+transpose(g))+lam*trace(g)*kronecker(mydomain)-sigma0