Search results for query: *

Hi, maki1337. I can see you sorted out whatever problem you had. But I'm just curious about why you should specify the roots of your quadratic as "integer" (instead of "real"). I don't think you can always know before hand if you'll get integer roots or otherwise. Unless you're only interested...

Hi mikrom, Thanks for the modifications and the extra index-finding functions you introduced. I did the adaptations to my problem and everything works perfectly well! Thank you so much again.

Hi mikrom, I have tried the idea you suggested above. But your post has highlighted the mistake I made of not using the appropriate indexing of the input data and the point of interest. I'll follow your suggestions and try to make my procedure work.

Hi mikrom, I am not able to attach the files themselves. So I'll just paste the data below. The file.dat.400 on is the higher resolution which is the main input file (I call the y-column there Y_hi). The file.dat.200 is the lower resolution file from which I only need the X's. The aim is to...

Hi all. I have a Lagrange program that uses a 4-point data set to interpolate some non-tabulated value in between any of those data points. I want to extend this script to act on a large input file, call it file1.dat, where the data is arranged in two columns Xin Yin. Then the script should...

...I am getting impression that you are confusing variable values with variable coefficients." Oh, I agree salgerman. My language use was very sloppy (or even inaccurate) and hence the confusion between matrix coefficients and variable. Also, the coefficient matrix I printed out in my original...

Yes, looking at the coefficient matrix, there equations are intriguing for sure. Just some brief background of the problem. This is part of an evolution problem which uses implicit time-stepping (integrator). So this involves a Newton-Raphson step at some point. It is in this step that we then...

mikrom,I'm sorry the example program you gave me (http://www.tek-tips.com/viewthread.cfm?qid=1678628 ) couldn't work because I had failed to properly install LAPACK on my system. These were the error messages from the two options you gave me: (1) /tmp/cczJi5YX.o: In function `MAIN__'...

Yes gummibaer, I have tried a countless number of times to get LAPACK onto my system but failed. My PC is running on Linux Ubuntu 11.10 & uses the Gnome desktop (&terminal). My compiler is "gfortran". Even though my professors don not fully approve of using "black-box" solvers, I always wish to...

Those numbers are practically zeros, even when printed to 16 decimal places. They are just floating point artifacts.

I'm trying to solve the following linear (6x6) system that contains a row of zeros in it using LU decomposition and LU back substitution. coefficient matrix = 0.07192 -0.05664 -0.22443 -159.36051 0.00000 -0.01644 0.00000 -0.00000 -0.00000 -0.00000 0.00000 -0.00000 -0.14818...

Has anyone ever used LAPACK in their fortran(90) codes? I believe this can help me compare the results from my hand-coded linear-algebra codes - or even replace these in case LAPACK gives far better performance. Small pieces of example code might be helpful. I've also tried installing this...

Here are the type declarations in "sf_1d_mod" below. The numerical values of a, b, x, y, p and q are computed in some other subroutine before they are "pointed" to in this "backEuler_time_step" one. Even though I haven't done a specific check, I believe they do find their ways into the arrays...

When I compile and link the following piece of code in the relevant program (main), all works perfectly and the executable is formed successfully. But when running the executable I find that the matrix "A" in "Ax=b" is singular (only zeros, in fact). I suspect this is due to the way I have used...

Thank you for the assurance and extra tip, gummibaer. When I said "retrieve", I actually meant to say "access". I need to access the 2D arrays (6x6 matrices) for each i, and use it as the matrix "A" in Ax=b and solve a system of equations. In the point you made concerning this, would I be...

Hi. I'm working with 1D and 2D f90 arrays. Actually, for a fixed value of a looping/iteration index, I have to form a 2D (dimension(6,6) ) array that holds some Jacobian matrix. Calling this iteration counter i, where i ranges from 0 to n (say n=5000), this means I need to compute the 6x6...

Oops, the central difference estimation is amazing! Thanks indeed, mikrom!!!

Hi. I'm implementing a forward finite difference approximation of the Jacobian matrix (2x2). I have adapted a subroutine (fdjac) from somewhere.The issue is that this subroutine requires explicit interfaces because of the deferred shapes of the involved arrays. 1. May I please have tips on...

That's the point. But I had one. Eventually, I have somehow managed to reproduce the error message. I noticed that if at some point you mess up slightly in the linking sequence e.g. using the " -o " option without supplying the name of the output executable, the main program object file gets...

This was the error message: /usr/lib/gcc/i686-linux-gnu/4.6.1/../../../i386-linux-gnu/crt1.o: In function `_start': (.text+0x18): undefined reference to `main' collect2: ld returned 1 exit status I'm sorry for not including it in my initial post. I only wrote it down somewhere. I believe it...