fortran

问题 How would I insert a vector b into a matrix at column col? I cannot find any syntax for and insert or append function in Fortran. So far all I have done is reassigned the values in the column, but I only want to insert the vector. real :: M(n,n) integer :: n, col real :: b(n) M(n:col) = b(:) 回答1: If I understood of your problem, you want to: increase the number n of columns of the matrix m by 1; insert the content of a vector b in m as a new column, at index col ; right-shift the remaining

问题 How would I insert a vector b into a matrix at column col? I cannot find any syntax for and insert or append function in Fortran. So far all I have done is reassigned the values in the column, but I only want to insert the vector. real :: M(n,n) integer :: n, col real :: b(n) M(n:col) = b(:) 回答1: If I understood of your problem, you want to: increase the number n of columns of the matrix m by 1; insert the content of a vector b in m as a new column, at index col ; right-shift the remaining

问题 I am trying to read a data file written by a Fortran program, in which every once in a while there is a very small float like 0.3299880-104 . The error message is: >np.loadtxt(filename, usecols = (1,)) File "/home/anaconda2/lib/python2.7/site-packages/numpy/lib/npyio.py", line 928, in loadtxt items = [conv(val) for (conv, val) in zip(converters, vals)] File "/home/anaconda2/lib/python2.7/site-packages/numpy/lib/npyio.py", line 659, in floatconv return float(x) ValueError: invalid literal for

问题 I am trying to read a data file written by a Fortran program, in which every once in a while there is a very small float like 0.3299880-104 . The error message is: >np.loadtxt(filename, usecols = (1,)) File "/home/anaconda2/lib/python2.7/site-packages/numpy/lib/npyio.py", line 928, in loadtxt items = [conv(val) for (conv, val) in zip(converters, vals)] File "/home/anaconda2/lib/python2.7/site-packages/numpy/lib/npyio.py", line 659, in floatconv return float(x) ValueError: invalid literal for

问题 In C I can easily return a value from a function: int foo(int b) { if (b == 0) return 42; int a; // calculate a return a; } But in Fortran the RETURN statement serves error handling. I could do integer function foo(b) integer :: b, a if (b == 0) ! what should I enter here?? // calculate a foo = a end function How do I do this in modern Fortran? I know that in this case and if-then-else-endif would suffice. But there are cases when it wouldn't and I don't want to make an overly complex example

问题 This is my code: Program Output_Format Implicit none Integer::k Integer,parameter:: Br_nn_mre = 5 Character(56),parameter:: FMT_01 = '(1x,"NN_DM:",1x,*("NN_",i2.2,1x))' Open( 15 , File = 'Output.txt' , Status = 'Unknown' , Action = 'Write' ) Write( 15 , FMT_01 ) ( k , k = 1 , Br_nn_mre ) Close( 15 , Status = 'Keep' ) End Program Output_Format The content of Output.txt file is: NN_DM: NN_01 NN_02 NN_03 NN_04 NN_05 NN_ I want to get this content in Output.txt : NN_DM: NN_01 NN_02 NN_03 NN_04 NN

问题 I want to define a pointer to a subarray. For a simple range this is easily done by pointer => array(i:j) , but I can't figure out how to do this for a map like k=[k1,k2,k3] . If I would define another array I could use a loop like array2=[(array1(k(j)),j=1,size(k,1))] . But it isn't possible to assign a pointer in a similar way ( pointer => [(array1(k(j)),j=1,size(k,1))] ) since the r.h.s. of the expression seems to define another variabel which then not even has the target attribute. For

问题 I want to run a Fortran program within Octave. I would like to do this for automation purposes and use Octave for all the data processing. Is it possible to run a Fortran program from octave using cygwin, if so, could you provide me some pointers along that direction? Moreover, I have a gfortran compiler installed in my system, Is there a way I could make use of it to complete my task mentioned above? Furthermore, I tried to use mex to perform the same: mckoctfile --mex HelloWorld.f I got the

问题 I want to run a Fortran program within Octave. I would like to do this for automation purposes and use Octave for all the data processing. Is it possible to run a Fortran program from octave using cygwin, if so, could you provide me some pointers along that direction? Moreover, I have a gfortran compiler installed in my system, Is there a way I could make use of it to complete my task mentioned above? Furthermore, I tried to use mex to perform the same: mckoctfile --mex HelloWorld.f I got the

问题 I am using Fortran 90 with the gfortran compiler as part of cygwin. I want to write a function that will create a series of new folders into a directory that is also passed as a parameter along with a number that is the maximum number of new consecutively numbered folders. Since I have to declare the length of the characters (ie strings) but also want to universally be able to pass different paths, I tried to pass the trimmed strings to the function. program main implicit none character(len =