x=schisq(vis,modvis,err) Make (data-model)/error for lsqnonlin
0001 function x=schisq(vis,modvis,err) 0002 % x=schisq(vis,modvis,err) 0003 % 0004 % Make (data-model)/error for lsqnonlin 0005 0006 % vectorize 0007 vis=cvec(vis); 0008 modvis=cvec(modvis); 0009 err=cvec(err); 0010 0011 % remove flagged data 0012 i=~isnan(vis); 0013 vis=vis(i); 0014 modvis=modvis(i); 0015 err=err(i); 0016 0017 D=vis-modvis; 0018 0019 if(~isreal(D)) 0020 D=[real(D);imag(D)]; 0021 0022 if(~isreal(err)) 0023 err=[real(err);imag(err)]; 0024 else 0025 err=[err;err]; 0026 end 0027 end 0028 0029 x=D./err; 0030 0031 return