Diff of /trunk/finley/test/python/time_chunks.py

revision 1565 by gross, Thu May 22 10:19:47 2008 UTC revision 1566 by gross, Thu May 22 12:07:35 2008 UTC
# Line 108  if error>REL_TOL*Lsup(u_ex): raise Runti Line 108  if error>REL_TOL*Lsup(u_ex): raise Runti
108
109  for n in [1000, 10000, 50000, 100000]:  for n in [1000, 10000, 50000, 100000]:
110  # for n in [1000, 10000]:  # for n in [1000, 10000]:
111   for prop in [ (1,2), (2,2), (1,3), (2,3) ]:   #for prop in [ (1,2), (2,2), (1,3), (2,3) ]:
112     for prop in [ (1,2), (1,3) ]:
113     for tp in [ "s", "v" ]:     for tp in [ "s", "v" ]:
114        # create code:        # create code:
116        dim=prop[1]        dim=prop[1]
117        if isinstance(prop[0], int):        if isinstance(prop[0], int):
118            o=prop[0]            o=prop[0]
119            NE=int(float(n)**(1./dim)/o)+1            if tp=="s":
120            q=1
121          else:
122            q=dim
123              NE=int(float(n/dim)**(1./dim)/o)+1
124            prog+="NE=%d\n"%NE            prog+="NE=%d\n"%NE
125            if dim==2:            if dim==2:
126                if o==1:                if o==1:
# Line 138  for n in [1000, 10000, 50000, 100000]: Line 143  for n in [1000, 10000, 50000, 100000]:
143             prog+=TEST_3_s             prog+=TEST_3_s
144          else:          else:
145             prog+=TEST_3_v             prog+=TEST_3_v
146        print "l= %d, dim= %d, type=%s, order=%s"%((o*NE+1)**dim,dim,tp,o)        print "l= %d, dim= %d, type=%s, order=%s"%(q*(o*NE+1)**dim,dim,tp,o)
147
148        prog+=SOLVE_AND_TEST        prog+=SOLVE_AND_TEST
149        # run code:        # run code:

Legend:
 Removed from v.1565 changed lines Added in v.1566