Starting diary file diary-02-Apr-2018.txt

what

MATLAB Code files in the current folder /Users/gilbert/Documents/CS 219 2018 Spring/Matlab

clustereg  gploth     jacobi     misluby    misseq     startup    

MAT-files in the current folder /Users/gilbert/Documents/CS 219 2018 Spring/Matlab

clusters     heatproblem  meshes       


load heatproblem
whos
  Name           Size                Bytes  Class     Attributes

  A          10000x10000            873608  double    sparse    
  ans            1x50                  100  char                
  b          10000x1                 80000  double              
  d              1x25                   50  char                
  dfile          1x21                   42  char                

nnz(A)
ans =
       49600
A
A =
         (1,1)              4
         (2,1)             -1
       (101,1)             -1
      ...
      (9999,10000)         -1
     (10000,10000)          4

spy(A)
plot(b)

x = A \ b;

x
x =
    0.0588
    0.1181
    0.1784
    ...
    0.0938
    0.0751
    0.0564

clear

what

MATLAB Code files in the current folder /Users/gilbert/Documents/CS 219 2018 Spring/Matlab

clustereg  gploth     jacobi     misluby    misseq     startup    

MAT-files in the current folder /Users/gilbert/Documents/CS 219 2018 Spring/Matlab

clusters     heatproblem  meshes       


load clusters
whos
  Name          Size            Bytes  Class     Attributes

  A            50x50             9400  double    sparse    
  S             4x4                72  double    sparse    
  ans         100x1               800  double              
  d             1x25               50  char                
  dfile         1x21               42  char                
  ignore       50x1               400  double              
  lambda        2x1                16  double              
  p             1x50              400  double              
  pp           50x1               400  double              
  uv           50x2               800  double              
  xy           50x2               800  double              

A
A =
   (1,1)        9
   (4,1)       -1
   ...
  (49,50)      -1
  (50,50)      13

spy(A)
nnz(A - A')
ans =
     0

xy
xy =
         0    1.0000
    0.1253    0.9921
    0.2487    0.9686
    0.3681    0.9298
    ...
   -0.5878    0.8090
   -0.4818    0.8763
   -0.3681    0.9298
   -0.2487    0.9686
   -0.1253    0.9921

n = length(A)
n =
    50

gploth(A,xy,1:n)

help eigs
  ...

help fiedler
  ...

[uv, lambda] = fiedler(A,2);
size(uv)
ans =
    50     2
lambda
lambda =
    0.8595
    0.9956
uv
uv =
   -0.1861    0.0589
    0.0537   -0.1840
    0.1924    0.2106
    ...
   -0.1515    0.0652
   -0.1781    0.0492

u = uv(:,1);
v = uv(:,2);

u' * v
ans =
  -4.9960e-16

u' * u
ans =
    1.0000

gploth(A,uv,1:n)

spy(A)
shg

plot(u)
[u, perm]= sort(u)

u =
   -0.1894
   -0.1861
   -0.1858
   ...
    0.1871
    0.1924

perm =
     4
     1
    38
    22
    40
    ...
    37
    16
     8
     3

size(perm)
ans =
    50     1

plot(u)
Ap = A(perm,perm);
spy(Ap)

clear
load meshes

whos
  Name              Size               Bytes  Class     Attributes

  Airfoil        4253x4253            495328  double    sparse    
  Amap           4253x1                34024  double              
  Axy            4253x2                68048  double              
  Eppstein        547x547              63248  double    sparse    
  Exy             547x2                 8752  double              
  Rmat           1024x1024            349096  double    sparse    
  Rmat64           64x64                9096  double    sparse    
  Rxy            1024x2                16384  double              
  Rxy64            64x2                 1024  double              
  Smallmesh       136x136              14600  double    sparse    
  Sxy             136x2                 2176  double              
  Tapir          1024x1024            115656  double    sparse    
  Txy            1024x2                16384  double              
  d                 1x25                  50  char                
  dfile             1x21                  42  char                

A = Airfoil;
gplotg(A,Axy)
spy(A)

R = chol(A);

figure,spy(R)

n = length(A)
n =
        4253
perm = amd(A);

figure, spy(A(perm,perm))
figure, spy(chol(A(perm,perm)))