Search:

Type: Posts; User: Richard Hanson

Page 1 of 6 1 2 3 4

Search: Search took 0.02 seconds.

  1. This is a software limitation for SURFND. For n...

    This is a software limitation for SURFND. For n > 7 your question remains open. You can look into an alternate code found in Chap. 12.2 of http://mathalacarte.com/c/math77_head.html. There is no...
  2. For n

    For n <= 7 see the routine in FNL 7, SURFND: Performs multidimensional interpolation and differentiation for up to 7 dimensions. This routine also can be used for numerical derivatives. If...
  3. Replies
    2
    Views
    2,266

    The Exponent of a matrix is not something that is...

    The Exponent of a matrix is not something that is recognizable to me. If you mean the Exponential of a matrix, then there is a nice paper by Cleve Moler and Charles van Loan, Nineteen Dubious Ways...
  4. Try this formulation so you can use a real...

    Try this formulation so you can use a real equation integrator:

    Define y(t)= u(t) + i v(t), where y is an n-vector, and so are u and v. But u and v are both real. Your problem is y'=f(y,t), say....
  5. Replies
    8
    Views
    7,859

    I am violating my admonition to the poster: This...

    I am violating my admonition to the poster: This posted code is an illustration only and shows how to call the C routine from Fortran. No support is intended. If you feel the C routine is...
  6. Replies
    8
    Views
    7,859

    The three files attached illustrate the...

    The three files attached illustrate the inter-language call of Fortran calling the IMSL C stat code for the multi-variate normal CDF. Understand that this is not supported by IMSL. Please do not...
  7. Replies
    7
    Views
    3,109

    As already mentioned, absorb the constraints in...

    As already mentioned, absorb the constraints in your mathematical model first. Then in the function evaluation routine, F, put the complex variables together using the symmetry conditions and the...
  8. Replies
    8
    Views
    7,859

    I will work up an example and post it later. ...

    I will work up an example and post it later. Other readers may benefit, so it is worth while to go to some work. What you will see has these elements:

    - You define a Fortran derived type that...
  9. Replies
    13
    Views
    5,438

    This is apparently a request for a vector...

    This is apparently a request for a vector integration. Each component function has the same domain. The short answer is that there is no direct support for this problem. You will have to call the...
  10. Replies
    8
    Views
    7,859

    Fortran does not have a multivariate normal cdf...

    Fortran does not have a multivariate normal cdf routine. But C does, as you note. If you have access to the C library, the Fortran 2003 standard supports an intrinsic module for inter-operating. ...
  11. Replies
    13
    Views
    5,438

    It is vague what you want to do here. Is the...

    It is vague what you want to do here. Is the routine called output2 to be a variable procedure itself? That sort of thing can be done using procedure pointers. If it is a fixed routine, you can...
  12. Replies
    7
    Views
    3,109

    OK - First thing is to mathematically simplify...

    OK - First thing is to mathematically simplify the problem. From what you say, this amounts to x(1)=conjg(x(5)), x(2)=conjg(x(4)), x(3) = real. Therefore you generally will have an odd number, n, ...
  13. Replies
    13
    Views
    5,438

    The value "a" is visible to your INTEGRATE...

    The value "a" is visible to your INTEGRATE routine, too. Try it. As long as you are happy with QDAGI, use it. Otherwise you can transform the integration domain.
  14. Replies
    7
    Views
    3,109

    More detail, please: 1. f(x) is a real...

    More detail, please:

    1. f(x) is a real function of the complex variables, x?
    2. A is a matrix of n complex columns vectors?
    3. What does Ax < b mean for the complex quantities? Real part of...
  15. Replies
    13
    Views
    5,438

    Put the type declaration for "a" right after the...

    Put the type declaration for "a" right after the IMPLICIT NONE in modulo1. Remove the typing for "a" from your main program and the function F. Assign the value as you have done. The key idea is...
  16. Replies
    13
    Views
    5,438

    With the organization you have, the "a" can be...

    With the organization you have, the "a" can be defined and given a type within the module (modulo1) and is thus visible to the function, F. You can chain a use-association of yet another module...
  17. Replies
    7
    Views
    3,109

    Best if you give an example of this using 5...

    Best if you give an example of this using 5 variables, say, and an objective function. Lots of possibilities here but it is good to be clear first. Put in everything you wish - bounds,...
  18. Replies
    13
    Views
    5,438

    Use-associate the module that describes the...

    Use-associate the module that describes the interface to QDAGI: use qdagi_int. Also remove the type declaration for your function, F. Its type is known by its placement as a double precision module...
  19. Replies
    1
    Views
    2,047

    There should be no problem for the larger size. ...

    There should be no problem for the larger size. Suggest you contact post-sales support. One candidate for the problem is that the environment is not set correctly. If you encountered a...
  20. Thread: Ivpag

    by Richard Hanson
    Replies
    2
    Views
    4,266

    Without seeing the behavior with numerical...

    Without seeing the behavior with numerical results it is difficult to suggest what might solve your problem. But for your sake please do not use IVPAG any longer. The code DAESL will integrate...
  21. Replies
    2
    Views
    5,275

    A statement of how you want to apply parallelism...

    A statement of how you want to apply parallelism is essential if a useful suggestion is to be given. This routine uses reverse communication for the matrix-vector products of the method. Presumably...
  22. Replies
    1
    Views
    2,016

    The names are generic. To illustrate, run the...

    The names are generic. To illustrate, run the document example but use double precision data types.

    (There may be a slight typo in the documentation.)

    The specific name is D_BSJ0() for...
  23. You are missing an important point. There is no...

    You are missing an important point. There is no need to do the algebraic preliminary such as you have written. Don't go there.

    Consider the example problem of computing the inverse sine of w=.49...
  24. Replies
    9
    Views
    4,454

    As suggested, this formulation will allow...

    As suggested, this formulation will allow combining linear and non-linear constraints. Call the combined constraints non-linear but take advantage of the linear parts during gradient evaluation. ...
  25. Replies
    1
    Views
    2,053

    I suggest you contact post-sales support. Or you...

    I suggest you contact post-sales support. Or you can post a snip of your application that shows the problem.
    This will require a running program with all parts included.
Results 1 to 25 of 141
Page 1 of 6 1 2 3 4