From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <cygwin-return-200204-listarch-cygwin=sourceware.org@cygwin.com>
Received: (qmail 15941 invoked by alias); 23 Nov 2015 09:00:46 -0000
Mailing-List: contact cygwin-help@cygwin.com; run by ezmlm
Precedence: bulk
List-Id: <cygwin.cygwin.com>
List-Subscribe: <mailto:cygwin-subscribe@cygwin.com>
List-Archive: <http://sourceware.org/ml/cygwin/>
List-Post: <mailto:cygwin@cygwin.com>
List-Help: <mailto:cygwin-help@cygwin.com>, <http://sourceware.org/ml/#faqs>
Sender: cygwin-owner@cygwin.com
Mail-Followup-To: cygwin@cygwin.com
Received: (qmail 15904 invoked by uid 89); 23 Nov 2015 09:00:45 -0000
Authentication-Results: sourceware.org; auth=none
X-Virus-Found: No
X-Spam-SWARE-Status: No, score=1.2 required=5.0 tests=AWL,BAYES_50,KAM_LAZY_DOMAIN_SECURITY,RCVD_IN_DNSWL_LOW autolearn=no version=3.3.2
X-HELO: BAY004-OMC1S14.hotmail.com
Received: from bay004-omc1s14.hotmail.com (HELO BAY004-OMC1S14.hotmail.com) (65.54.190.25) by sourceware.org (qpsmtpd/0.93/v0.84-503-g423c35a) with (AES256-SHA256 encrypted) ESMTPS; Mon, 23 Nov 2015 09:00:41 +0000
Received: from BAY169-DS1 ([65.54.190.60]) by BAY004-OMC1S14.hotmail.com over TLS secured channel with Microsoft SMTPSVC(7.5.7601.23008);	 Mon, 23 Nov 2015 01:00:39 -0800
X-TMN: [fdFPDIb3Z9mEbzVVGNklRnVvBbfy7axI]
Message-ID: <BAY169-DS11E052CA89ABFD18A7398A7070@phx.gbl>
From: "Tony Kelman" <tony@kelman.net>
To: <cygwin@cygwin.com>
References: <BAY169-DS260B02A7A24CFA816E17BA7070@phx.gbl> <5652B582.20405@gmail.com>
In-Reply-To: <5652B582.20405@gmail.com>
Subject: Re: Octave 4.0.0-3 x86_64 segfault (suspect openblas)
Date: Mon, 23 Nov 2015 09:00:00 -0000
MIME-Version: 1.0
Content-Type: multipart/mixed;	boundary="----=_NextPart_000_0145_01D1258A.43CA6FF0"
X-IsSubscribed: yes
X-SW-Source: 2015-11/txt/msg00327.txt.bz2

------=_NextPart_000_0145_01D1258A.43CA6FF0
Content-Type: text/plain;
	format=flowed;
	charset="Windows-1252";
	reply-type=response
Content-Transfer-Encoding: 7bit
Content-length: 278

> To test without openblas , you can renable netlib blas
> $ mv /usr/bin/cygblas-0.dll  /usr/bin/cygblas-0.dll_bk

Yep, works fine with netlib blas. Pure C repro case is attached,
gcc -g 5728.c -o 5728 -llapack && gdb ./5728

I can submit upstream, unless you'd like to.

-Tony

------=_NextPart_000_0145_01D1258A.43CA6FF0
Content-Type: text/plain;
	format=flowed;
	name="5728.c";
	reply-type=response
Content-Transfer-Encoding: quoted-printable
Content-Disposition: attachment;
	filename="5728.c"
Content-length: 1193

#include <stdint.h>
#include <stdlib.h>
#include <stdio.h>
#define BLASINT int32_t
void dgeev_(char*, char*, BLASINT*, double*, BLASINT*, double*, double*,
    double*, BLASINT*, double*, BLASINT*, double*, BLASINT*, BLASINT*);
int main() {
char jobvl =3D 'N';
char jobvr =3D 'N';
BLASINT n =3D 1000;
double *A =3D malloc(n*n*sizeof(double));
BLASINT lda =3D n;
double *WR =3D malloc(n*sizeof(double));
double *WI =3D malloc(n*sizeof(double));
double *VL;
BLASINT ldvl =3D n;
double *VR;
BLASINT ldvr =3D n;
double *work =3D malloc(1*sizeof(double));
BLASINT lwork =3D -1;
BLASINT info[1];
BLASINT i, j;
for (i=3D0; i<n; i++) {
    for (j=3D0; j<n; j++) {
        A[n*i+j] =3D n*(((double)rand()/(double)RAND_MAX) - 0.5);
    }
}
dgeev_(&jobvl, &jobvr, &n, A, &lda, WR, WI,
    VL, &ldvl, VR, &ldvr, work, &lwork, info);
printf("info[0] is %d\n", info[0]);
printf("lwork is %d\n", lwork);
printf("work[0] is %g\n", work[0]);
lwork =3D work[0];
free(work);
work =3D malloc(lwork*sizeof(double));
dgeev_(&jobvl, &jobvr, &n, A, &lda, WR, WI,
    VL, &ldvl, VR, &ldvr, work, &lwork, info);
printf("info[0] is %d\n", info[0]);
printf("lwork is %d\n", lwork);
printf("work[0] is %g\n", work[0]);
}


------=_NextPart_000_0145_01D1258A.43CA6FF0
Content-Type: text/plain; charset=us-ascii
Content-length: 218

--
Problem reports:       http://cygwin.com/problems.html
FAQ:                   http://cygwin.com/faq/
Documentation:         http://cygwin.com/docs.html
Unsubscribe info:      http://cygwin.com/ml/#unsubscribe-simple
------=_NextPart_000_0145_01D1258A.43CA6FF0--