Commit 694e2cbf authored by matecellab's avatar matecellab
Browse files

updated comments after consistency check by Willeke

parent afe3fc09
......@@ -106,7 +106,7 @@ for b = 1:numel(basil)
% OK, let's not optimise at all, and rely on the mmmoi for now
% com and eigenvectors of the mass moment of inertia tensor
[com,ev] = bas_getIndDir(X(Li),Y(Li),Z(Li)); %
[com, ev] = bas_getIndDir(X(Li),Y(Li),Z(Li)); %
% indenter is mostly in the xy-plane, I hope
idir = ev(:,ev(3,:) == min(ev(3,:)) );
% indenter always points in the positive y-direction, I hope
......@@ -124,7 +124,7 @@ for b = 1:numel(basil)
dxyz = [ix(dp == min(dp)), iy(dp == min(dp)), iz(dp == min(dp))];
top = com - idir'*sqrt(sum(dxyz.^2));
% output: [indenter centre of mass (com); cylinder axis vector]
% output: [indenter tip estimate; cylinder axis vector]
tipNR = [tip; R];
if param.figs % gimme a plot
......
......@@ -27,7 +27,7 @@ indSE = zeros(2,3);
for b = 1:numel(basil)
tipNR = load(fullfile(param.basilhome, num2str(basil(b)),strcat('Indenter_tipNR.txt')));
tip = tipNR(1,:); % indenter centre of mass
tip = tipNR(1,:); % indenter tip estimate
R = tipNR(2:4,:); % indenter rotation matrix
idir = -R(3,:)'; % indenter long axis orientation
......@@ -42,7 +42,7 @@ for b = 1:numel(basil)
%% indenter start position
% offset undeformed skin to indenter tip
% offset undeformed skin to indenter tip (estimate)
V = Vs - repmat(tip,size(Vs,1),1);
% find point 'most' in the opposite direction of indenter movement
Vl = sqrt(sum(V.^2,2));
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment