コード例 #1
0
ファイル: cart3d.py プロジェクト: hurlei/pyNastran
    def _read_header_binary(self):
        data = self.infile.read(4)
        size_little, = unpack(b'<i', data)
        size_big, = unpack(b'>i', data)
        if size_big in [12, 8]:
            self._endian = b'>'
            size = size_big
        elif size_little in [8, 12]:
            self._endian = b'<'
            size = size_little
        else:
            self._rewind()
            self.show(100)
            raise RuntimeError('unknown endian')

        self.n += 4
        data = self.infile.read(size)
        self.n += size

        so4 = size // 4  # size over 4
        if so4 == 3:
            (npoints, nelements, nresults) = unpack(self._endian + b('iii'), data)
            self.log.info("npoints=%s nelements=%s nresults=%s" % (npoints, nelements, nresults))
        elif so4 == 2:
            (npoints, nelements) = unpack(self._endian + b('ii'), data)
            nresults = 0
            self.log.info("npoints=%s nelements=%s" % (npoints, nelements))
        else:
            self._rewind()
            self.show(100)
            raise RuntimeError('in the wrong spot...endian...size/4=%s' % so4)
        self.infile.read(8)  # end of first block, start of second block
        return (npoints, nelements, nresults)
コード例 #2
0
    def _read_header_binary(self):
        data = self.infile.read(4)
        size_little, = unpack(b'<i', data)
        size_big, = unpack(b'>i', data)
        if size_big in [12, 8]:
            self._endian = b'>'
            size = size_big
        elif size_little in [8, 12]:
            self._endian = b'<'
            size = size_little
        else:
            self._rewind()
            self.show(100)
            raise RuntimeError('unknown endian')

        self.n += 4
        data = self.infile.read(size)
        self.n += size

        so4 = size // 4  # size over 4
        if so4 == 3:
            (npoints, nelements, nresults) = unpack(self._endian + b('iii'),
                                                    data)
            self.log.info("npoints=%s nelements=%s nresults=%s" %
                          (npoints, nelements, nresults))
        elif so4 == 2:
            (npoints, nelements) = unpack(self._endian + b('ii'), data)
            nresults = 0
            self.log.info("npoints=%s nelements=%s" % (npoints, nelements))
        else:
            self._rewind()
            self.show(100)
            raise RuntimeError('in the wrong spot...endian...size/4=%s' % so4)
        self.infile.read(8)  # end of first block, start of second block
        return (npoints, nelements, nresults)
コード例 #3
0
ファイル: cart3d.py プロジェクト: mnekkach/pyNastran
    def _write_regions(self, outfile, regions, is_binary):
        """writes the regions"""
        if is_binary:
            fmt = self._endian + b('i')
            four = pack(fmt, 4)
            outfile.write(four)

            nregions = len(regions)
            fmt = self._endian + b('%ii' % nregions)
            ints = pack(fmt, *regions)
            outfile.write(ints)

            outfile.write(four)
        else:
            fmt = b'%i'
            np.savetxt(outfile, regions, fmt)
コード例 #4
0
ファイル: cart3d.py プロジェクト: hurlei/pyNastran
    def _write_regions(self, outfile, regions, is_binary):
        """writes the regions"""
        if is_binary:
            fmt = self._endian + b('i')
            four = pack(fmt, 4)
            outfile.write(four)

            nregions = len(regions)
            fmt = self._endian + b('%ii' % nregions)
            ints = pack(fmt, *regions)
            outfile.write(ints)

            outfile.write(four)
        else:
            fmt = b'%i'
            np.savetxt(outfile, regions, fmt)
コード例 #5
0
ファイル: cart3d.py プロジェクト: mnekkach/pyNastran
    def _write_points(self, outfile, points, is_binary, float_fmt='%6.6f'):
        """writes the points"""
        if is_binary:
            four = pack(self._endian + b('i'), 4)
            outfile.write(four)

            npoints = points.shape[0]
            fmt = self._endian + b('%if' % (npoints * 3))
            floats = pack(fmt, *np.ravel(points))

            outfile.write(floats)
            outfile.write(four)
        else:
            if isinstance(float_fmt, bytes_type):
                fmt_ascii = float_fmt
            else:
                fmt_ascii = float_fmt.encode('latin1')
            np.savetxt(outfile, points, fmt_ascii)
コード例 #6
0
ファイル: cart3d.py プロジェクト: hurlei/pyNastran
    def _write_points(self, outfile, points, is_binary, float_fmt='%6.6f'):
        """writes the points"""
        if is_binary:
            four = pack(self._endian + b('i'), 4)
            outfile.write(four)

            npoints = points.shape[0]
            fmt = self._endian + b('%if' % (npoints * 3))
            floats = pack(fmt, *np.ravel(points))

            outfile.write(floats)
            outfile.write(four)
        else:
            if isinstance(float_fmt, bytes_type):
                fmt_ascii = float_fmt
            else:
                fmt_ascii = float_fmt.encode('latin1')
            np.savetxt(outfile, points, fmt_ascii)
コード例 #7
0
ファイル: cart3d.py プロジェクト: mnekkach/pyNastran
    def _read_points_binary(self, npoints):
        """reads the xyz points"""
        size = npoints * 12  # 12=3*4 all the points
        data = self.infile.read(size)

        dtype = np.dtype(self._endian + b('f4'))
        points = np.frombuffer(data, dtype=dtype).reshape((npoints, 3)).copy()

        self.infile.read(8)  # end of second block, start of third block
        return points
コード例 #8
0
ファイル: cart3d.py プロジェクト: mnekkach/pyNastran
    def _write_header(self,
                      outfile,
                      points,
                      elements,
                      is_loads,
                      is_binary=False):
        """
        writes the cart3d header

        Without results
        ---------------
        npoints nelements

        With results
        ------------
        npoints nelements nresults

        """
        npoints = points.shape[0]
        nelements = elements.shape[0]

        if is_binary:
            if is_loads:
                fmt = self._endian + b('iiiii')
                msg = pack(fmt, 3 * 4, npoints, nelements, 6, 4)
            else:
                fmt = self._endian + b('iiii')
                msg = pack(fmt, 2 * 4, npoints, nelements, 4)

            int_fmt = None
        else:
            # this is ASCII data
            if is_loads:
                msg = b("%i %i 6\n" % (npoints, nelements))
            else:
                msg = b("%i %i\n" % (npoints, nelements))

            # take the max value, string it, and length it
            # so 123,456 is length 6
            int_fmt = b('%%%si' % len(str(nelements)))
        outfile.write(msg)
        return int_fmt
コード例 #9
0
ファイル: cart3d.py プロジェクト: mnekkach/pyNastran
    def _write_elements(self, outfile, elements, is_binary, int_fmt='%6i'):
        """writes the triangles"""
        min_e = elements.min()
        assert min_e == 0, 'min(elements)=%s' % min_e
        if is_binary:
            fmt = self._endian + b('i')
            four = pack(fmt, 4)
            outfile.write(four)
            nelements = elements.shape[0]
            fmt = self._endian + b('%ii' % (nelements * 3))
            ints = pack(fmt, *np.ravel(elements + 1))

            outfile.write(ints)
            outfile.write(four)
        else:
            if isinstance(int_fmt, bytes_type):
                fmt_ascii = int_fmt
            else:
                fmt_ascii = int_fmt.encode('latin1')
            np.savetxt(outfile, elements + 1, fmt_ascii)
コード例 #10
0
ファイル: cart3d.py プロジェクト: hurlei/pyNastran
    def _write_elements(self, outfile, elements, is_binary, int_fmt='%6i'):
        """writes the triangles"""
        min_e = elements.min()
        assert min_e == 1, 'min(elements)=%s' % min_e
        if is_binary:
            fmt = self._endian + b('i')
            four = pack(fmt, 4)
            outfile.write(four)
            nelements = elements.shape[0]
            fmt = self._endian + b('%ii' % (nelements * 3))
            ints = pack(fmt, *np.ravel(elements))

            outfile.write(ints)
            outfile.write(four)
        else:
            if isinstance(int_fmt, bytes_type):
                fmt_ascii = int_fmt
            else:
                fmt_ascii = int_fmt.encode('latin1')
            np.savetxt(outfile, elements, fmt_ascii)
コード例 #11
0
ファイル: cart3d.py プロジェクト: upendra117/pyNastran
    def _read_elements_binary(self, nelements):
        """reads the triangles"""
        size = nelements * 12  # 12=3*4 all the elements
        data = self.infile.read(size)

        dtype = np.dtype(self._endian + b('i4'))
        elements = np.fromstring(data, dtype=dtype).reshape((nelements, 3))

        self.infile.read(
            8)  # end of third (element) block, start of regions (fourth) block
        assert elements.min() == 1, elements.min()
        return elements - 1
コード例 #12
0
ファイル: cart3d.py プロジェクト: hurlei/pyNastran
    def _write_header(self, outfile, points, elements, is_loads, is_binary=False):
        """
        writes the cart3d header

        Without results
        ---------------
        npoints nelements

        With results
        ------------
        npoints nelements nresults

        """
        npoints = points.shape[0]
        nelements = elements.shape[0]

        if is_binary:
            if is_loads:
                fmt = self._endian + b('iiiii')
                msg = pack(fmt, 3*4, npoints, nelements, 6, 4)
            else:
                fmt = self._endian + b('iiii')
                msg = pack(fmt, 2*4, npoints, nelements, 4)

            int_fmt = None
        else:
            # this is ASCII data
            if is_loads:
                msg = b("%i %i 6\n" % (npoints, nelements))
            else:
                msg = b("%i %i\n" % (npoints, nelements))

            # take the max value, string it, and length it
            # so 123,456 is length 6
            int_fmt = b('%%%si' % len(str(nelements)))
        outfile.write(msg)
        return int_fmt