cleanup code

Author: sbryngelson

NS_2d_cylinder_PHS.m

@@ -17,7 +17,7 @@
 % U0    =  W0(1:L_W2/2);
 % V0    =  W0(L_W2/2+1:end);
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
 
 %% advection: adams-bashforth method
 
@@ -46,8 +46,6 @@
  V1(end-L_B+1:end-0) = V(end-L_B+1:end-0);
 %  
 %  
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %% viscousity: Crank-Nicolson
 
 % 
@@ -72,8 +70,6 @@
 
 U2   =  L_u_inv(U2);
 V2   =  L_v_inv(V2);
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %% pressure correction 
 
    F0    =  (D0_12_x*(U2-0)+D0_12_y*(V2-0));
@@ -85,7 +81,7 @@
   p = p(1:(length(F0)+L_B_S));  % regularization
   F = F(1:(length(F0)+L_B_S));  % regularization
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
 
     U3   =   (U2-0) -  [(D0_21_x*p ); zeros(L_B,1)];
 
@@ -101,8 +97,6 @@
 
 % V3(L_W+1:end) = V(L_W+1:end);
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
 %% velocities at the following time-step
 
 W3   =   [U3;V3];

RBF_PHS_FD_all.m

@@ -1,16 +1,12 @@
 function [D_all] = RBF_PHS_FD_all(xy1,xy_s,Nearest_Idx,k,m,d)
 
 % RBF_PHS_FD_ALL  RBF-FD-based (PHS+poly) global differentiation matrices
-
-
 %   [D_ALL] = RBF_PHS_FD_ALL(XY1,XY_S,NEAREST_IDX,K,M,D,VARARGIN) returns
 %   the global differentiation matrices from the node set {xy_s} to node set
 %   {xy1}, including Dx, Dy, DL, Dxx, Dyy, Dxy.
 %   NEAREST_IDX is the index matrix return by NEAREST_INTERP.
 %   k is local RBF stencil size, m is the order of PHS-RBF, and d is the
 %   degree of the polynomial Augmentation.
-
-
 %   Reference:
 %     [1] T. Chu, O. T. Schmidt, RBF-FD discretization of the Navier-Stokes 
 %     equations on scattered but staggered nodes,
@@ -19,10 +15,6 @@
 %     Submitted to Journal of Computational Physics
 %
 %   T. Chu (tic173@ucsd.edu), O. T. Schmidt (oschmidt@ucsd.edu)
-
-
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%
 X1            =           xy1(:,1);
 NxNy1         =           length(X1);
@@ -35,8 +27,6 @@
 weight_yy     =           zeros(NxNy1,k);
 weight_xy     =           zeros(NxNy1,k);
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
  %%  Local RBF-FDs for each grid   
 
 
@@ -139,8 +129,6 @@
 % A_inv =    V*S*(U');
 % 
 % W = A_inv*L;
-
-
     w = W(1:k,:);      % Extragct the RBF-FD weights
 
     weight_x(m1,1:k)   =  w(:,1)/scale_x;
@@ -158,24 +146,16 @@
 
 
 end
-
-
 %%  Assigning RBF-FD weights to global matrices
-
-
 Idx_x                                 =       reshape((1:NxNy1).*ones(k,1),1,[]);
 
 Idx_y                                 =       reshape(Nearest_Idx(:,1:k)',1,[]);
-
-
 Wx                                    =       reshape(weight_x',1,[]);
 Wy                                    =       reshape(weight_y',1,[]);
 WL                                    =       reshape(weight_L',1,[]);
 Wxx                                   =       reshape(weight_xx',1,[]);
 Wyy                                   =       reshape(weight_yy',1,[]);
 Wxy                                   =       reshape(weight_xy',1,[]);
-
-
 Dx                                    =       sparse(Idx_x,Idx_y,Wx,NxNy1,NxNy2);
 Dy                                    =       sparse(Idx_x,Idx_y,Wy,NxNy1,NxNy2);
 L                                     =       sparse(Idx_x,Idx_y,WL,NxNy1,NxNy2);
@@ -184,6 +164,4 @@
 Dxy                                   =       sparse(Idx_x,Idx_y,Wxy,NxNy1,NxNy2);
 
 D_all                                 =       {Dx,Dy,L,Dxx,Dyy,Dxy};
-
-
 end