# Search:

Type: Posts; User: Richard Hanson

Page 1 of 6 1

1. ## Thread: Interpolation in more than 2 dimensions?

by Richard Hanson
Replies
3
Views
1,642

### 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. ## Thread: Interpolation in more than 2 dimensions?

by Richard Hanson
Replies
3
Views
1,642

### 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. ## Thread: exponents of a matrix

by Richard Hanson
Replies
2
Views
1,644

### 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. ## Thread: IMSL for solving ODE with complex coefficients

by Richard Hanson
Replies
1
Views
2,624

### 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. ## Thread: multivariate normal cdf?

by Richard Hanson
Replies
8
Views
4,969

### 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. ## Thread: multivariate normal cdf?

by Richard Hanson
Replies
8
Views
4,969

### 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. ## Thread: Constraints on the solution

by Richard Hanson
Replies
7
Views
1,978

### 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. ## Thread: multivariate normal cdf?

by Richard Hanson
Replies
8
Views
4,969

### 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. ## Thread: QDAGI using a module?

by Richard Hanson
Replies
13
Views
3,536

### 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. ## Thread: multivariate normal cdf?

by Richard Hanson
Replies
8
Views
4,969

### 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. ## Thread: QDAGI using a module?

by Richard Hanson
Replies
13
Views
3,536

### 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. ## Thread: Constraints on the solution

by Richard Hanson
Replies
7
Views
1,978

### 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. ## Thread: QDAGI using a module?

by Richard Hanson
Replies
13
Views
3,536

### 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. ## Thread: Constraints on the solution

by Richard Hanson
Replies
7
Views
1,978

### 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. ## Thread: QDAGI using a module?

by Richard Hanson
Replies
13
Views
3,536

### 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. ## Thread: QDAGI using a module?

by Richard Hanson
Replies
13
Views
3,536

### 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. ## Thread: Constraints on the solution

by Richard Hanson
Replies
7
Views
1,978

### 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. ## Thread: QDAGI using a module?

by Richard Hanson
Replies
13
Views
3,536

### 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. ## Thread: Error with the linear operator .i

by Richard Hanson
Replies
1
Views
1,550

### 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
3,363

### 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. ## Thread: Parallel Programming with IMSL?

by Richard Hanson
Replies
2
Views
4,758

### 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. ## Thread: IMSL 7 Double precision

by Richard Hanson
Replies
1
Views
1,555

### 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. ## Thread: Can IMSL Solve a system of n nonlinear complex equations

by Richard Hanson
Replies
5
Views
6,035

### 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. ## Thread: ¿fmincon in IMSL?

by Richard Hanson
Replies
9
Views
3,419

### 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. ## Thread: Problem with LCONF

by Richard Hanson
Replies
1
Views
1,612

### 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