X-Recipient: archive-cygwin AT delorie DOT com
DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id
	:list-unsubscribe:list-subscribe:list-archive:list-post
	:list-help:sender:message-id:from:to:references:in-reply-to
	:subject:date:mime-version:content-type; q=dns; s=default; b=Hv+
	6/7xWAG+R6qexKCezo1I3xsi7iGt3ZE0GUzm224o+0PTL+db5LkUzyzIr5qXuTf+
	yppNVebyBx6JaX4K3b5WXEWgG4INj35sCxjEJkYhdn0yBbOPK455TlkOZU6a1iqZ
	v67FphHcnja92qd45f52pguSF0sCVkxKqBLYolZs=
DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id
	:list-unsubscribe:list-subscribe:list-archive:list-post
	:list-help:sender:message-id:from:to:references:in-reply-to
	:subject:date:mime-version:content-type; s=default; bh=MIYut3hSQ
	42p57If+PJS4ffZlVM=; b=R+VmplqbM6+1asDhikLFmtXWIIfnsrZ/ow2E8L6ej
	3zLYv+eq6cMz2UIMxHUIjGd6rqyaoxL3A2UUuLa+xZntoxXFQ6pWkIKtb3T+cqKk
	OkDFmB9nzE0zAuyOvcYPl6/yQ3d+/4xOSQSPFFClXOiFoTMG20qH9JIwgsBZPQK0
	wQ=
Mailing-List: contact cygwin-help AT cygwin DOT com; run by ezmlm
List-Id: <cygwin.cygwin.com>
List-Subscribe: <mailto:cygwin-subscribe AT cygwin DOT com>
List-Archive: <http://sourceware.org/ml/cygwin/>
List-Post: <mailto:cygwin AT cygwin DOT com>
List-Help: <mailto:cygwin-help AT cygwin DOT com>, <http://sourceware.org/ml/#faqs>
Sender: cygwin-owner AT cygwin DOT com
Mail-Followup-To: cygwin AT cygwin DOT com
Delivered-To: mailing list cygwin AT cygwin DOT com
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
X-TMN: [fdFPDIb3Z9mEbzVVGNklRnVvBbfy7axI]
Message-ID: <BAY169-DS11E052CA89ABFD18A7398A7070@phx.gbl>
From: "Tony Kelman" <tony AT kelman DOT net>
To: <cygwin AT cygwin DOT com>
References: <BAY169-DS260B02A7A24CFA816E17BA7070 AT phx DOT gbl> <5652B582 DOT 20405 AT gmail DOT 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 00:59:54 -0800
MIME-Version: 1.0
Content-Type: multipart/mixed;	boundary="----=_NextPart_000_0145_01D1258A.43CA6FF0"
X-IsSubscribed: yes

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

> 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"

#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

--
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--