Recent content by Angelos K

  1. Angelos K

    Is LAPACK slow on my system? And if so why?

    I apologize. I do not know. How can I check?
  2. Angelos K

    Is LAPACK slow on my system? And if so why?

    I am a programming newbie, I hope I don't give you too much trouble. I use LAPACK on my 2.7GHz CPU. My expectations on the running times are based on Figure 4.1 in this dedicated publication , as well as on Fanfan's answer to this Stackoverflow Post : Both computers mentioned (and especially...
  3. Angelos K

    Cannot use simplest ifstream with LAPACK

    All errors disappeared by using C linkage. Every occurrence of extern void function {...} was replaced by extern "C" { void function {...} } Specifically, I now use: /* DSYEV prototype */ extern "C"{ void dsyev( char* jobz, char* uplo, int* n, double* a, int* lda, double* w, double* work...
  4. Angelos K

    Cannot use simplest ifstream with LAPACK

    Thank you so much! What you say is so logical and resolves many errors. But there are still some errors left: ludi@ludi-M17xR4:~/Desktop/tests$ g++ -Ddsyev=dsyev_ -o combo.x combo.cc -L/usr/local/lib -llapack -lblas && ./combo.x combo.cc: In function ‘int main()’: combo.cc:107:74: warning...
  5. Angelos K

    Cannot use simplest ifstream with LAPACK

    I have an official LAPACK example from here: https://software.intel.com/sites/products/documentation/doclib/mkl_sa/11/mkl_lapack_examples/dsyev_ex.c.htm It compiles and runs without errors or warnings after using ludi@ludi-M17xR4:~/Desktop/tests$ gcc -Ddsyev=dsyev_ -o sylapack sylapack.c...
  6. Angelos K

    LAPACK dgeev: parameter had illegal value

    Thank you very much for the clarification about [in] parameters! As you resolved all problems: Should I mark this thread as solved, and if so, how?
  7. Angelos K

    LAPACK dgeev: parameter had illegal value

    Thank you. It appears to work, but I have only tested trivial cases. I saw your warning earlier, but understood it just now (blush). It seems I am not able to modify my post anymore. So I should replace [ i ] by [k] evereywhere and then make the corresponding logical adjustments, right? What...
  8. Angelos K

    LAPACK dgeev: parameter had illegal value

    Thanks. This code appears to work. #include<stdio.h> #include<math.h> #include <stdlib.h> //................... void dgeTranspose( double *Transposed, double *M ,int n) { int i,j; for(i=0;i<n;i++) for(j=0;j<n;j++) Transposed[i+n*j] = M[i*n+j]; } //................... //...
  9. Angelos K

    LAPACK dgeev: parameter had illegal value

    Mod note: I revised the code below slightly, changing the loop control variable i to either j or k. The reason for this is that the browser mistakes the letter i in brackets for the BBCode italics tag, which causes some array expressions to partially disappear. Hello, I am trying for the first...
  10. Angelos K

    MATLAB Why are there saw-teeth on my Matlab eps prints?

    Thank you, unstable. I am using Linux. Release 2010b. How can I send you the .fig file? Greetings, Angelos
  11. Angelos K

    MATLAB Why are there saw-teeth on my Matlab eps prints?

    Dear all, I generate some MATLAB figures that I want to print as eps. The figures look fine in the .fig format. When I open the eps, there is a problem. Wherever there is a vertical line (or part of function) the line shows tiny saw-teeth. I would like to attach I file but it seems one can't...
  12. Angelos K

    Is the Error Building Cuba-2.1 Causing Faulty Outputs?

    Hi, all! Unfortunately I'm not much of a computer guy, so I'm sorry if I leave out anything important. I'm attempting to build Cuba-2.1. The distribution I'm using is: DISTRIB_ID=Ubuntu DISTRIB_RELEASE=10.04 DISTRIB_CODENAME=lucid DISTRIB_DESCRIPTION="Ubuntu 10.04.2 LTS" I go...
  13. Angelos K

    Multiple images in gravitational lensing

    Dear all, so far I have only studied weak lensing. Therefore I understand little about multiple images.The light paths, by Fermat's principle, have to correspond to stationary points of the (proper) time delay. Minima, maxima saddle points. Basics are nicely reviewed here [p.14]...
  14. Angelos K

    What are entropy perturbations?

    Hi, I can only offer my own thoughts based on some thermodynasmics, but the following seems very reasonable to me: The class of pertubation you allready know is the adiabatic one. IF we assume reversibility, adiabatic processes are exactly isentropic (fixed entropy) ones. It makes...
  15. Angelos K

    Is time measured at receding objects dilated?

    Related paper An interesting paper on the subject is this one: Davis and Lineweaver 2004 in Publications of the Astronomical Society of Australia, Volume 21, Issue 1, pp. 97-109 It also deals with the question why light from such superluminarily receding objects could reach us - which I...
Back
Top