| MATLAB Function Reference | |
Syntax
Description
D, N0
R,V,Fpo,Ft,F,v,r
Rth
dvdf: coefficient of proportionality: incr. of v vs. F
dt: r=r+v*dt;
define two functions: dforce(RR)-takes array, gives force;
cappa(RR)-takes array, outputs friction
dfdv: cappa vs. v; %coefficients
rscal=2;
dt=0.1;
dvdf=0.1*dt;
dfdv=1;
ons=ones(N,1);
Nn=200*N;
[min(sort(ij)) max(sort(ij))]
return
Ilen=length(I)
eps=zeros(Lii,1);
eps1(eps1<1)=1;
V(:,i)=v(ii,i)-v(jj,i);
Vn=Vn+(V(:,i).*nn(:,i)); %normal velocity
capp=cappa(RR(I),dfdv);
cappPA=capp; % parallel capp;
find(iii==0)
Vt=V-Vn*ones(1,D).*nn; %transversal velocity
Ft=-V.*(capp*ones(1,D));
Ft=-Vt.*(capp*ones(1,D)); %if cappa is a function of RR
Ft=-Vt*capp; %else: if cappa is constant;
FF{i}=sparse(ii,jj,Fpo(:,i),N,N);
dd(1,i)=F(I(l));
ss=full(ss);
max(abs(v(:)))
v=reshape(v,[N,D])
dd=[FF{1}(ii(l),jj(l)),FF{2}(ii(l),jj(l)),FF{3}(ii(l),jj(l))];
sclr=dd*nn(l,:)'
norms=norm(dd)
dVmin=min(dv(:))
v=v+dv+ RandPart*randn(size(v));
Remarks
See
Also
| |
diffuse1 | difmat | |
| Main Page - Institute for Genomics and Bioinformatics - Computableplant Project - MATLAB |