Tek-Tips is the largest IT community on the Internet today!
Members share and learn making Tek-Tips Forums the best source of peer-reviewed technical information on the Internet!
-
Congratulations Chriss Miller on being selected by the Tek-Tips community for having the most helpful posts in the forums last week. Way to Go!
-
-
Vegas integration problem
Thanks! Actually I found a strange thing: If I compile with gfortran actually looks that it work only with iseed=0; but if before call vegas i put a call srand(iseed) and in the vegas program i take of the argument of Ran(iseed), so if I write Ran() looks like that I can change seed, that is... -
Vegas integration problem
Hi, now should be ok :) program test implicit real*8 (a-h,o-z) external fxn COMMON/BVEG1/NCALL,ITMX,NPRN,NDEV,XL(100),XU(100),ACC COMMON/BVEG2/IT,NDO,SI,SWGT,SCHI,XI(50,100) COMMON/BVEG3/ALPH,NDMX,MDS COMMON/RNDM/ISEED ndim=1 ncall=10 itmx=1 nprn=-1... -
Vegas integration problem
Hi, thanks! Do you mean formatted in this way? Ah, I use gfortran to compile the code Thanks ************** *** MAIN **** ************** program test implicit real*8 (a-h,o-z) external fxn COMMON/BVEG1/NCALL,ITMX,NPRN,NDEV,XL(100),XU(100),ACC COMMON/BVEG2/IT,NDO,SI,SWGT,SCHI,XI(50,100)... -
Vegas integration problem
Hi, I have just started to use Vegas for numerical integration and I've found a problem that I can't get rid of... I'm trying, kust for test, to integrate a naive function, say x**2. The main function of my program and the vegas function are here in the post. My problem is the following: with... -
Matrix Diagonalization with constraint
Thanks! I'll try as soon as I can. It looks quite general (maybe I'm wrong... ) so I think I can adapt this event to shift in non diagonal elements! Hope so! Thanks again, I'll try with a 2x2 and the try to adapt to my problem!! -
Matrix Diagonalization with constraint
@ gummibaer Yes, is possible. I've done this with Mathematica, but I need to write the same program in fortran in order to connetc my output with some other programs. @ melmacianalf Can I easily find these funcions "min" and "eval" on the internete? Sorry maybe for these trivial questions... -
Matrix Diagonalization with constraint
Hi, I'm new in this forum and I hope I can find some help :) I'm dealing with the following problem. I've got a matrix that I'm diagonalizing with the subroutine Jacobi and Eigsrt that I've found online, and everithing works very fine. But now I need to impose some constraint on the...
