costFunction.m File Reference

Go to the source code of this file.

Functions
Brian	Staber (brian.staber @gmail.com) % function output
	writeXMLParameterList (xmlfilename, modelParameters.mu,... modelParameters.beta,... modelParameters.lc,... modelParameters.delta,... optimParameters.nmc,... optimParameters.station)
	if (optimParameters.station==44) executable
	if (s~=0) fprintf('Trilinos program failed.\n')
	Y (:, j+1)
	Yexp (:, k)
P	sum ()
output	m (i)
	if (isempty(pts)==0) [pdf
end end end	unix (strcat('rm/home/s/staber/Trilinos_results/nrl/random_generator_for_pca_likeliehood/station', num2str(optimParameters.station),'/*'))

Variables
end	xmlfilename = strcat('nrl.msme.station',num2str(optimParameters.station),'.xml')
else	executable
end	s = unix(executable)
	return
end	theta = ['15'
	theta_to_id = [5,6
output	eta = cell(4,1)
output	etaExp = cell(4,1)
output	m = zeros(4,1)
output	fval = 0
for	i
	Yexp = zeros(2355,2)
for	j
	filename = strcat(path, '/RandomVariableY_angle=',num2str(theta(i,:)),'_nmc=',num2str(j),'.mtx')
end	meanY = mean(Y,2)
	meanYexp = mean(Yexp,2)
output	Y {i} = Y
for	k
	exx = dlmread(filename)
	eyy = dlmread(filename)
	exy = dlmread(filename)
	CovarianceMatrix = cov(Y')
	L = L(:,idx)
	err = [1
for	l
	pts = find(output.etaExp{i}(l,:)>=min(supp) & output.etaExp{i}(l,:)<=max(supp))

Function Documentation

if() [1/3]

if

(

optimParameters.

station = =44

)

if() [2/3]

if

(

s~

= 0

)

if() [3/3]

if

(

isempty(pts)

= =0

)

m()

output m

(

i

)

Staber()

Brian Staber

(

brian.staber @gmail.

com

)

sum()

P sum ( )

virtual

unix()

end end end unix

(

strcat( 'rm/home/s/staber/Trilinos_results/nrl/random_generator_for_pca_likeliehood/station', num2str(optimParameters.station),'/*')

)

writeXMLParameterList()

writeXMLParameterList	(	xmlfilename	,
		modelParameters.	mu,
		... modelParameters.	beta,
		... modelParameters.	lc,
		... modelParameters.	delta,
		... optimParameters.	nmc,
		... optimParameters.	station
	)

Y()

Y	(	:	,
		j+	1
	)

Yexp()

Yexp	(	:	,
		k
	)

Variable Documentation

CovarianceMatrix

CovarianceMatrix = cov(Y')

Definition at line 63 of file costFunction.m.

err

err = [1

Definition at line 67 of file costFunction.m.

eta

output eta = cell(4,1)

Definition at line 33 of file costFunction.m.

etaExp

output etaExp = cell(4,1)

Definition at line 34 of file costFunction.m.

executable

else executable

Initial value:

= strcat('/usr/local/apps/petsc-3.7.6/gnu-mpich/bin/mpirun -np', 32, num2str(optimParameters.np), 32, './trilinos_mpi_', ...
                           num2str(optimParameters.station), 32, '--xml-in-file="',xmlfilename,'"')

Definition at line 21 of file costFunction.m.

exx

exx = dlmread(filename)

Definition at line 54 of file costFunction.m.

exy

exy = dlmread(filename)

Definition at line 58 of file costFunction.m.

eyy

eyy = dlmread(filename)

Definition at line 56 of file costFunction.m.

filename

filename = strcat(path, '/RandomVariableY_angle=',num2str(theta(i,:)),'_nmc=',num2str(j),'.mtx')

Definition at line 44 of file costFunction.m.

fval

output fval = 0

Definition at line 36 of file costFunction.m.

i

for i

Initial value:

= 1:4
        Y    = zeros(2355,optimParameters.nmc)

Definition at line 38 of file costFunction.m.

j

for j

Initial value:

= 0:optimParameters.nmc-1
            path     = strcat('/home/s/staber/Trilinos_results/nrl/random_generator_for_pca_likelihood/station',num2str(optimParameters.station))

Definition at line 42 of file costFunction.m.

k

for k

Initial value:

= 1:2
%             id = theta_to_id(i,k)

Definition at line 51 of file costFunction.m.

L

L = L(:,idx)

Definition at line 66 of file costFunction.m.

l

for l

Initial value:

= 1:output.m(i)
            output.eta{i}(l,:)    = ( (Y    - repmat(meanY,1,optimParameters.nmc))'*L(:,l) )/sqrt(P(l))

Definition at line 73 of file costFunction.m.

m

output m = zeros(4,1)

Definition at line 35 of file costFunction.m.

meanY

end meanY = mean(Y,2)

Definition at line 47 of file costFunction.m.

meanYexp

meanYexp = mean(Yexp,2)

Definition at line 48 of file costFunction.m.

pts

pts = find(output.etaExp{i}(l,:)>=min(supp) & output.etaExp{i}(l,:)<=max(supp))

Definition at line 77 of file costFunction.m.

return

return

Definition at line 27 of file costFunction.m.

s

end s = unix(executable)

Definition at line 24 of file costFunction.m.

theta

end theta = ['15'

Definition at line 30 of file costFunction.m.

theta_to_id

theta_to_id = [5,6

Definition at line 31 of file costFunction.m.

xmlfilename

end xmlfilename = strcat('nrl.msme.station',num2str(optimParameters.station),'.xml')

Definition at line 10 of file costFunction.m.

Y

output Y {i} = Y

Definition at line 49 of file costFunction.m.

Yexp

end output Yexp = zeros(2355,2)

Definition at line 40 of file costFunction.m.