def METHOD_NAME(self):
    self._update()
    has_z = self._ndim == 3
    n = self._ndim
    m = self.__len__()
    array_type = c_double * (m * n)
    data = array_type()
    temp = c_double()
    for i in range(m):
        lgeos.GEOSCoordSeq_getX(self._cseq, i, byref(temp))
        data[n*i] = temp.value
        lgeos.GEOSCoordSeq_getY(self._cseq, i, byref(temp))
        data[n*i+1] = temp.value
        if has_z:
            lgeos.GEOSCoordSeq_getZ(self._cseq, i, byref(temp))
            data[n*i+2] = temp.value
    return data