Mail Archives: cygwin/2015/11/23/04:01:03
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
|
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--
- Raw text -