Displaying 3 results from an estimated 3 matches for "eigensolut".
1999 Mar 31
0
ace () - again
...ading because I thought I might finish the current
project before it but if this is suspicious I go for the most recent
version now).
Any pointers or ideas are appreciated. Regards, Lorenz
g77 -O2 -fpic -c ace.f -o ace.o
ace.f: In subroutine `mace':
ace.f:322: warning:
590 format( 15h0eigensolution i2, 1h:)
^
Missing comma in FORMAT statement at (^)
ace.f:323: warning:
600 format( 15h eigensolution i2, 23h r**2 = 1 - e**2 =g12.4)
^
Missing comma in FORMAT statement at (^)
ace.f:323: warning:
600...
2010 Jul 10
1
Not nice behaviour of nlminb (windows 32 bit, version, 2.11.1)
...ut that does not mean that the authors ever thought anyone would use their large
package to solve a small 1D problem, or if they did, whether they seriously tested for the
awkward things that happen when one only has one parameter. I've seen similar things go
wrong with matrix packages e.g., eigensolutions of 1 by 1 matrices.
So this msg is to watch carefully when problems are n=1 and put in some checks that
answers make sense.
JN
2011 May 11
1
Problem with constrained optimization with maxBFGS
Dear all,
I need to maximize the v:
v= D' W D
D is a column vector ( n , 1)
W is a given matrix (n, n)
subject to:
sum D= 1
(BTW, n is less than 300)
I´ve tried to use maxBFGS, as follows:
#####################################
objectiveFunction<-function(x)
{
return(t(D)%*%W%*%D)
}
Amat<-diag(nrow(D))
Amat<-rbind((rep(-1, nrow(D))), Amat)
bvec<-matrix( c(0), nrow(D)+1,