forums.silverfrost.com Forum Index forums.silverfrost.com
Welcome to the Silverfrost forums
 
 FAQFAQ   SearchSearch   MemberlistMemberlist   UsergroupsUsergroups   RegisterRegister 
 ProfileProfile   Log in to check your private messagesLog in to check your private messages   Log inLog in 

Search found 913 matches
forums.silverfrost.com Forum Index
Author Message
  Topic: Calculating the outer vector product
mecej4

Replies: 21
Views: 299

PostForum: Support   Posted: Sun Jun 24, 2018 12:16 pm   Subject: Re:
I understand that SPREAD is not elemental and that FTN95 is not doing what the programmer intended. At the moment I am not clear whether FTN95 can do anything to correct this. It may not even be able ...
  Topic: Calculating the outer vector product
mecej4

Replies: 21
Views: 299

PostForum: Support   Posted: Sat Jun 23, 2018 1:28 am   Subject:
The Fortran 95 standard documents the properties of the arguments and result value of intrinsic procedures in section 13.14. For each such procedure, the second attribute listed is "Class". ...
  Topic: Calculating the outer vector product
mecej4

Replies: 21
Views: 299

PostForum: Support   Posted: Fri Jun 22, 2018 12:23 pm   Subject: Poor performance of SPREAD function
Here is a simplified version of John Campbell's test program, with only the SPREAD based version of the outer product computation being included. I hope that this version succeeds in convincing reader ...
  Topic: Calculating the outer vector product
mecej4

Replies: 21
Views: 299

PostForum: Support   Posted: Thu Jun 21, 2018 11:48 am   Subject:
The computation of the outer product of two m X n matrices involves m.n multiplications plus 3 m.n load/store operations.

There is no need to do matrix multiplication in this. Firstly, the two inp ...
  Topic: Pure devilry catch in debugger
mecej4

Replies: 6
Views: 465

PostForum: 64-bit   Posted: Mon May 14, 2018 6:52 am   Subject:
Checking input data should include checking that data is input AND that the input values are correct. If a program uses undefined variables, those are variables that should have been given values earl ...
  Topic: Error Argument of ichar at (1) must be of length one
mecej4

Replies: 8
Views: 460

PostForum: Support   Posted: Sun May 13, 2018 11:35 am   Subject:
There are two places where ICHAR is used with an argument of length other than 1: lines 166 and 243. The program may have been developed and used with compilers that allowed such sloppy usage. Either ...
  Topic: Error Argument of ichar at (1) must be of length one
mecej4

Replies: 8
Views: 460

PostForum: Support   Posted: Sun May 13, 2018 10:04 am   Subject: Re:
Geomarine, your code extract still leaves out the declarations of the variables that occur in the line flagged by the compiler as erroneous.

However, if Paul's guess about E_det_strg is correct, t ...
  Topic: Error Argument of ichar at (1) must be of length one
mecej4

Replies: 8
Views: 460

PostForum: Support   Posted: Sun May 13, 2018 9:54 am   Subject:
Geomarine, your code extract still leaves out the declarations of the variables that occur in the line flagged by the compiler as erroneous.

However, if Paul's guess about E_det_strg is correct, t ...
  Topic: Error Argument of ichar at (1) must be of length one
mecej4

Replies: 8
Views: 460

PostForum: Support   Posted: Sat May 12, 2018 7:58 am   Subject:
The intrinsic function ICHAR takes as argument a character of length 1. Your code seems to be passing instead a REAL or double precision argument.

How to fix? That depends on what is intended to b ...
  Topic: Why we have to declare random@ ?
mecej4

Replies: 15
Views: 823

PostForum: Support   Posted: Tue May 08, 2018 5:00 pm   Subject: Re:
It is interesting that RANDOM@ works as both real*4 and real*8 for 32-bit ( if I tested it correctly)
That is the way that the 80X87 coprocessor behaves. You can often declare in the caller a wrong ...
  Topic: Why we have to declare random@ ?
mecej4

Replies: 15
Views: 823

PostForum: Support   Posted: Mon May 07, 2018 4:00 pm   Subject: Re:
It took all of my self-control not to start a thread entitled 'When is random not random', subtitled 'Some randoms are more random than other randoms'.
Eddie
Sorry, George Marsaglia beat you to it ...
  Topic: Why we have to declare random@ ?
mecej4

Replies: 15
Views: 823

PostForum: Support   Posted: Mon May 07, 2018 3:41 am   Subject:
Eddie, your points are well taken. Nonstandard functions, if provided by a vendor, need not be intrinsic. If they are not, and they are not also generic, the user must provide a type declaration if im ...
  Topic: Why we have to declare random@ ?
mecej4

Replies: 15
Views: 823

PostForum: Support   Posted: Sun May 06, 2018 2:15 pm   Subject:
program tbess
implicit none
! real, intrinsic :: bessel_j0 ! seems to have no effect
real :: x
! x = bessel_j0(1.0) ! works fine
x = bessel_j0() ...
  Topic: Why we have to declare random@ ?
mecej4

Replies: 15
Views: 823

PostForum: Support   Posted: Sun May 06, 2018 2:25 am   Subject:
You leave us guessing as to how you reached these conclusions. Example code, please. Do you need to use RANDOM@ at all, instead of calling the standard intrinsic subroutine RANDOM_NUMBER?

I do find ...
  Topic: FTN95 beta testing
mecej4

Replies: 18
Views: 853

PostForum: General   Posted: Sun May 06, 2018 1:48 am   Subject: Re:
Dan ... I DID use the code above !!!!
John, I think that there is something mixed up regarding which file you used to produce the plot and which compiler version you used. Add to that the chances for ...
 
Page 1 of 61 Goto page 1, 2, 3 ... 59, 60, 61  Next
All times are GMT + 1 Hour
Jump to:  


Powered by phpBB © 2001, 2005 phpBB Group