problem calling sparse MKL library in a user macro

[M. Kurdi]

Hello FEAP Users,
I would like to use the sparse routines from Intel MKL library. There seems to be memory overwriting error that occurs when the routine is called from FEAP. To help resolve this problem a small routine is attached to illustrate. The input data file calls the routine for matrix vector operation. The result retrieved is incorrect. Attached the same routine but in standalone form that performs same operation. Now the result is correct.
I am wondering if I am missing something in regards to the user macro or perhaps if anyone encountered similar problem. Please let me know.

[FEAP_Admin]

Your standalone program is not a standalone program.  How did you run this?

That said, the codes are essentially the same but in one you initialize Y and in the other you do not.
Is this the problem?

Otherwise I see no reason for the failure.

[M. Kurdi]

Both programs compiled using gfortran and linked to
-L /usr/global/intel/mkl/11.0.0.079/lib/intel64/ -lmkl_rt -lX11 -lm -g. In FEAP this is just inserted in the makefile.in

 

[M. Kurdi]

The problem is due to the -i8 flag for "ifort" or -fdefault-integer-8 flag in gfortran. If these are removed then the program can be executed correctly within gdb but crashes when run separately. This could be due to the fact that the mkl library is pre-built based on short integers. However I am puzzled why running without gdb crashes.

[FEAP_Admin]

Why are you using -i8?  This something you only need to do if you have extremely large problems.  Otherwise you should be using 32-bit integers (note this is independent of whether or not you have 64-bit pointers and are using include/integer8).

[M. Kurdi]

Correct, I'm not using this anymore. In addition I had erroneously set ipr=1 in feap84.f Now life is good

[Brian]

How to link umacr19 to feap by visual studio?

[FEAP_Admin]

add the files to your project and rebuild.

[Brian]

umacr19 has been added to project in visual studio and feap has been rebuiled but there is the following error when Ifile is run. How to manage the problem, please?

  Inconsistency occurred from PMACIO in data file Ifile                                                                                                                           
  at or near record number    27.  Input record is:

  mumv conv !END OPTIMIZATION LOOP                                             

  If this record is correct error may result from
  missing blank record before new command type.

[Prof. R.L. Taylor]

If your umacr19 is the  mumv command then that is what is causing the error --  feap does not see your command?  If you use
the following command while in an interactive mode you can see all the available commands that feap knows

> MANUal,,5
> HELP

The first sets the command to show all manual commands to a maxiumum,  otherwise you only see a subset of the available commands. 

[Brian]

Calling manual,,5 and help in an interactive mode does not show mumv command because umacr19 cannot be linked to feap. How can umacr19 be linked to feap in visual studio? I did the following steps but it does not work.

1. Add umacr19 to feap project
2. Rebuild feap

[pjarz]

Are you sure that you have removed (or replaced) the default file umacr19.f from "user" folder?

[Brian]

Yes, I am sure. I did all steps as mentioned in post http://feap.berkeley.edu/forum/index.php?board=7.0.