xxt_test2.m File Reference

Go to the source code of this file.

Functions
*Als	eye (ns)]
	norm (M-em *[inv(All) 0 *Als;0 *Als'Qs *inv(Qs'*dS *Qs)*Qs']*em') Rf0
	X (:, i)
pause end	norm (M-em *[inv(All) 0 *Als;0 *Als'Qs *X *X'*Qs']*em') norm(M-em *[inv(All) 0 *Als

Variables
	p
	M
end	Al0 =[]
	Ac0 =zeros(2)([],:)
	As0 =[1 -.5
	Al1
	Ac1 =[-.5 0
	As1 =[1 -.5
	A0 =[Al0 Ac0
	A1 =[Al1 Ac1
	Il =eye(4)
	Is =eye(2)
	gI =eye(6)
	Rl0 =Il([],:)
	Rl1 =Il([1 2 3 4],:)
	Rs0 =Is([5 6]-4,:)
	Rs1 =Is([5 6]-4,:)
	R0
	R1
	A =R0'*A0*R0+R1'*A1*R1
	All
Rs1	ns = size(Qs)(1)
	Q = bdiag(Il,Qs,[])
Qs *Als Qs *Ass *	Qs
	S0 = As0
	S1 = As1-Ac1'*(Al1\Ac1)
	dS = bdiag(S0,S1,[])
	em =[Il -All\Als
	Rf1 =Is([5 6]-4,:)
	Qf = [Rf0
	X0 = zeros(size(Rs0)(1),size(Rf0)(1))
	X1 = zeros(size(Rs1)(1),size(Rf1)(1))
	dX = bdiag(X0,X1,[])
	X = Qs'*dX*Qf
	n
	Ri = eye(n)([1:i-1],:)
	se

Function Documentation

* Als eye

(

ns

)

norm

(

M-em **em'

[inv(All) 0 *Als;0 *Als'Qs *inv(Qs'*dS *Qs)*Qs']

)

pause end norm

(

M-em **em'

[inv(All) 0 *Als;0 *Als'Qs *X *X'*Qs']

)

X	(	:	,
		i
	)

Variable Documentation

A =R0'*A0*R0+R1'*A1*R1

Definition at line 39 of file xxt_test2.m.

A0 =[Al0 Ac0

Definition at line 25 of file xxt_test2.m.

A1 =[Al1 Ac1

Definition at line 26 of file xxt_test2.m.

Ac0 =zeros(2)([],:)

Definition at line 15 of file xxt_test2.m.

Ac1 =[-.5 0

Definition at line 22 of file xxt_test2.m.

end Al0 =[]

Definition at line 14 of file xxt_test2.m.

Al1

Initial value:

=[  2 -.5  -1   0
     -.5   1   0 -.5
      -1   0   2 -.5
       0 -.5 -.5   1]

Definition at line 18 of file xxt_test2.m.

All

Initial value:

= bdiag(Al0,Al1,[])
Als = bdiag(Ac0,Ac1,[])
Ass = bdiag(As0,As1,[])
Qs = [Rs0

Definition at line 41 of file xxt_test2.m.

Ac0 As0 =[1 -.5

Definition at line 16 of file xxt_test2.m.

Ac1 As1 =[1 -.5

Definition at line 23 of file xxt_test2.m.

dS = bdiag(S0,S1,[])

Definition at line 56 of file xxt_test2.m.

dX = bdiag(X0,X1,[])

Definition at line 67 of file xxt_test2.m.

em =[Il -All\Als

Definition at line 58 of file xxt_test2.m.

gI =eye(6)

Definition at line 29 of file xxt_test2.m.

Il =eye(4)

Definition at line 28 of file xxt_test2.m.

Is =eye(2)

Definition at line 28 of file xxt_test2.m.

M

Initial value:

= [ A            zeros(ra,cb) zeros(ra,cc)
        zeros(rb,ca) B            zeros(rb,cc)
        zeros(rc,ca) zeros(rc,cb) C ]

Definition at line 9 of file xxt_test2.m.

n

Initial value:

= size(Qs)(2)
for i = [1:n]
  ei = eye(n)(:,i)

Definition at line 72 of file xxt_test2.m.

Rs1 ns = size(Qs)(1)

Definition at line 45 of file xxt_test2.m.

p

Initial value:

= [4 3 2 1 3 6 1 5 6  5 ]

%inv(A)(p,p)

function M=bdiag(A,B,C)
  [ra ca]=size(A)

Definition at line 1 of file xxt_test2.m.

Referenced by allreduce_map_setup(), apply_m_Asl(), apply_p_Als(), apply_QQt(), apply_S(), apply_S_col(), condense_matrix(), crs_solve(), discover_dofs(), discover_sep_ids(), factor_numeric(), factor_symbolic(), fcrystal_setup(), fgs_fields(), findpts(), findpts_el_2(), findpts_el_3(), findpts_local(), findpts_local_eval(), findpts_local_setup(), flagged_primaries_map(), gll_lag_setup(), gs_gather_many(), gs_gather_vec_to_many(), gs_init_many(), gs_scatter_many(), gs_scatter_many_to_vec(), gs_scatter_vec_to_many(), heap_sortv(), lagrange_setup(), legendre(), legendre_d1(), legendre_d2(), lob_bnd_setup(), local_map(), main(), make_topology_unique(), merge_sep_ids(), pw_comm_setup(), pw_exec_recvs(), pw_exec_sends(), pw_map_setup(), radix_passp_be(), seed(), shared_ids_aux(), sparse_cholesky_solve(), tensor_3t(), unique_ids(), unique_nonzero(), and unpack_int().

Q = bdiag(Il,Qs,[])

Definition at line 47 of file xxt_test2.m.

Qf = [Rf0

Definition at line 63 of file xxt_test2.m.

Qs* Als Qs* Ass* Qs

Definition at line 49 of file xxt_test2.m.

R0

Initial value:

=[Rl0 zeros(size(Rl0)(1),size(Rs0)(2))
    zeros(size(Rs0)(1),size(Rl0)(2)) Rs0]

Definition at line 34 of file xxt_test2.m.

R1

Initial value:

=[Rl1 zeros(size(Rl1)(1),size(Rs1)(2))
    zeros(size(Rs1)(1),size(Rl1)(2)) Rs1]

Definition at line 36 of file xxt_test2.m.

Rf1 =Is([5 6]-4,:)

Definition at line 62 of file xxt_test2.m.

Ri = eye(n)([1:i-1],:)

Definition at line 75 of file xxt_test2.m.

Rl0 =Il([],:)

Definition at line 30 of file xxt_test2.m.

Rl1 =Il([1 2 3 4],:)

Definition at line 31 of file xxt_test2.m.

Rs0 =Is([5 6]-4,:)

Definition at line 32 of file xxt_test2.m.

Rs1 =Is([5 6]-4,:)

Definition at line 33 of file xxt_test2.m.

S0 = As0

Definition at line 53 of file xxt_test2.m.

S1 = As1-Ac1'*(Al1\Ac1)

Definition at line 54 of file xxt_test2.m.

se

Initial value:

= dS*Qs*ei
  Xtse = dX'*se
  QQtXtse = Qf*Ri'*Ri*Qf'*Xtse
  Qy = Qs*ei - dX*QQtXtse
  ytsy = Qy'*dS*Qy
  Qx = Qy/sqrt(ytsy)
  xv = inv(Qs'*Qs)*Qs'*Qx

Definition at line 76 of file xxt_test2.m.

X = Qs'*dX*Qf

Definition at line 68 of file xxt_test2.m.

X0 = zeros(size(Rs0)(1),size(Rf0)(1))

Definition at line 65 of file xxt_test2.m.

X1 = zeros(size(Rs1)(1),size(Rf1)(1))

Definition at line 66 of file xxt_test2.m.