X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=blobdiff_plain;f=octave_packages%2Fsecs2d-0.0.8%2FQDDGOX%2FQDDGOXddcurrent.m;fp=octave_packages%2Fsecs2d-0.0.8%2FQDDGOX%2FQDDGOXddcurrent.m;h=c9603c97f87f36b35742fb6630bc18ebcca0cd77;hp=0000000000000000000000000000000000000000;hb=f5f7a74bd8a4900f0b797da6783be80e11a68d86;hpb=1705066eceaaea976f010f669ce8e972f3734b05

diff --git a/octave_packages/secs2d-0.0.8/QDDGOX/QDDGOXddcurrent.m b/octave_packages/secs2d-0.0.8/QDDGOX/QDDGOXddcurrent.m
new file mode 100644
index 0000000..c9603c9
--- /dev/null
+++ b/octave_packages/secs2d-0.0.8/QDDGOX/QDDGOXddcurrent.m
@@ -0,0 +1,39 @@
+function [current,divrg]=QDDGOXddcurrent(mesh,Sinodes,data,contacts,options);
+
+% [current,divrg]=QDDGOXddcurrent(mesh,Sinodes,data,contacts);
+constants
+Nelements = size(mesh.t,2);
+
+
+mob = Ufielddepmob(mesh,data.un,data.Fn,data.vsatn,data.mubn);
+
+if (~isfield("options","FD")|(options.FD==0))
+data.FDn=0;
+data.FDp=0;
+end
+
+An  = Uscharfettergummel(mesh,data.V(Sinodes)+data.G+data.FDn,mob);
+
+if (options.holes==1)
+mob = Ufielddepmob(mesh,data.up,data.Fp,data.vsatp,data.mubp);
+Ap  = Uscharfettergummel(mesh,-data.V(Sinodes)-data.Gp-data.FDp,mob);
+divrg     = An * data.n + Ap * data.p;
+else
+divrg     = An * data.n;
+end
+
+for con = 1:length(contacts)
+
+    cedges = [];
+    cedges=[cedges,find(mesh.e(5,:)==contacts(con))];
+    cnodes = mesh.e(1:2,cedges);
+    cnodes = [cnodes(1,:) cnodes(2,:)];
+    cnodes = unique(cnodes);
+
+    current(con) = sum(divrg(cnodes));
+
+end
+
+Is = q*data.us*data.Vs*data.ns;
+current = current * Is;
+