Hi Sam - interesting glitch. To make this easier to debug, can you make a smaller (less refined) version of your mesh, say with just a few thousand quads? Current version (367,418 quads) is a bit heavy for my debugger.
Ideally, if you could make a gmsh geometry input file, I could spin my own mesh, then trace what happens as mesh resolution increases. Nigel On Tuesday, December 19, 2017 at 12:39:42 PM UTC+11, SAM WANG wrote: > > Dear all, >> > Sorry I post the thread again. > >> >> I use Pointwise to do a 2D mesh, but it gives me this following error. >> Can anyone help me look at it? >> >> By the way, for compressible flow, how can I obtain temperature data from >> Pyfr? >> > > Can Pyfr do large eddy simulation (LES) or RANS (k-epsilon)? If it does, > how can I set those additional parameters? > >> >> Please help me out. >> >> The ini and mesh files are attached as follows, > > https://www.dropbox.com/sh/ky4bmkpoqz8nb5q/AAACNxSNKgkfOkyUNMU7MlXja?dl=0 > <https://www.google.com/url?q=https%3A%2F%2Fwww.dropbox.com%2Fsh%2Fky4bmkpoqz8nb5q%2FAAACNxSNKgkfOkyUNMU7MlXja%3Fdl%3D0&sa=D&sntz=1&usg=AFQjCNEUYIQlCWPlLz7TwRYrg1ALBelR5g> > > > > Really thank you all in advance. > > Kind regards, > Sam > > > Traceback (most recent call last): > File "/usr/local/bin/pyfr", line 11, in <module> > sys.exit(main()) > File "/usr/local/lib/python3.5/ > dist-packages/pyfr/__main__.py", line 110, in main > args.process(args) > File "/usr/local/lib/python3.5/dist-packages/pyfr/__main__.py", line > 235, in process_run > args, NativeReader(args.mesh), None, Inifile.load(args.cfg) > File "/usr/local/lib/python3.5/dist-packages/pyfr/__main__.py", line > 216, in _process_common > solver = get_solver(backend, rallocs, mesh, soln, cfg) > File "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/__init__.py", > line 16, in get_solver > return get_integrator(backend, systemcls, rallocs, mesh, initsoln, cfg) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/__init__.py", line > 46, in get_integrator > return integrator(backend, systemcls, rallocs, mesh, initsoln, cfg) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/std/controllers.py", > line 80, in __init__ > super().__init__(*args, **kwargs) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/std/controllers.py", > line 14, in __init__ > super().__init__(*args, **kwargs) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/std/steppers.py", > line 165, in __init__ > super().__init__(*args, **kwargs) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/std/steppers.py", > line 8, in __init__ > super().__init__(*args, **kwargs) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/std/base.py", line > 10, in __init__ > super().__init__(*args, **kwargs) > File "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/base.py", > line 58, in __init__ > self._init_system(systemcls, backend, rallocs, mesh, initsoln) > File "/usr/local/lib/python3.5/dist-packages/pyfr/integrators/base.py", > line 88, in _init_system > self.system = systemcls(*args, nreg=self.nreg, cfg=self.cfg) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/system.py", line > 58, in __init__ > int_inters = self._load_int_inters(rallocs, mesh, elemap) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/system.py", line > 116, in _load_int_inters > self.cfg) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/navstokes/inters.py", > line 13, in __init__ > super().__init__(be, lhs, rhs, elemap, cfg) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/baseadvecdiff/inters.py", > > line 12, in __init__ > super().__init__(be, lhs, rhs, elemap, cfg) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/baseadvec/inters.py", > line 16, in __init__ > self._gen_perm(lhs, rhs) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/baseadvecdiff/inters.py", > > line 36, in _gen_perm > self._elemap) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/inters.py", line > 15, in get_opt_view_perm > vm = _get_inter_objs(interside, mat, elemap) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/inters.py", line > 11, in _get_inter_objs > return [emap[type](eidx, fidx) for type, eidx, fidx, flags in > interside] > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/inters.py", line > 11, in <listcomp> > return [emap[type](eidx, fidx) for type, eidx, fidx, flags in > interside] > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/elements.py", > line 346, in get_scal_fpts_for_inter > rmap = self._srtd_face_fpts[fidx][eidx] > File "/usr/local/lib/python3.5/dist-packages/pyfr/util.py", line 138, in > __get__ > value = self.fn(instance) > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/elements.py", > line 109, in _srtd_face_fpts > for ffpts in self.basis.facefpts] > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/elements.py", > line 109, in <listcomp> > for ffpts in self.basis.facefpts] > File > "/usr/local/lib/python3.5/dist-packages/pyfr/solvers/base/elements.py", > line 108, in <listcomp> > return [[np.array(fuzzysort(pts.tolist(), ffpts)) for pts in plocfpts] > File "/usr/local/lib/python3.5/dist-packages/pyfr/nputil.py", line 96, > in fuzzysort > srtdidx[i:] = fuzzysort(arr, srtdidx[i:], dim + 1, tol) > File "/usr/local/lib/python3.5/dist-packages/pyfr/nputil.py", line 96, > in fuzzysort > srtdidx[i:] = fuzzysort(arr, srtdidx[i:], dim + 1, tol) > File "/usr/local/lib/python3.5/dist-packages/pyfr/nputil.py", line 85, > in fuzzysort > arrd = arr[dim] > IndexError: list index out of range > -- You received this message because you are subscribed to the Google Groups "PyFR Mailing List" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send an email to [email protected]. Visit this group at https://groups.google.com/group/pyfrmailinglist. For more options, visit https://groups.google.com/d/optout.
