This is a static copy of a profile report

Home

integrateGsubmatrix_Telles (31 calls, 0.280 sec)
Generated 18-Mar-2011 23:31:44 using cpu time.
M-function in file /Users/Robert/Documents/Work/Cardiff/postdoc/myrepo/NURBS/isoBEM/integrateGsubmatrix_Telles.m
Copy to new window for comparing multiple runs

Parents (calling functions)
No parent
Lines where the most time was spent

Line NumberCodeCallsTotal Time% TimeTime Plot
33
[Ttemp(i,j) Utemp(i,j)]=DBIEke...
7440.124 s44.4%
6
[gpt gwt]=lgwt(ngp,-1,1);
310.041 s14.8%
27
[jacob_xi, normals, r, dr, drd...
1860.021 s7.4%
14
gamBar=nthroot( srcXi*xiStar +...
310.021 s7.4%
35
end
3720.010 s3.7%
All other lines  0.062 s22.2%
Totals  0.280 s100% 
Children (called functions)

Function NameFunction TypeCallsTotal Time% TimeTime Plot
DBIEkernelsM-function7440.104 s37.0%
lgwtM-function310.031 s11.1%
getKernelParametersM-function1860.021 s7.4%
nthrootM-function620.021 s7.4%
NURBSbasisMEX-function5580 s0%
convertToParamSpaceM-function1860 s0%
convertToParentCoordSpaceM-function310 s0%
Self time (built-ins, overhead, etc.)  0.104 s37.0%
Totals  0.280 s100% 
M-Lint results
No M-Lint messages.
Coverage results
[ Show coverage for parent directory ]
Total lines in function40
Non-code lines (comments, blank lines)11
Code lines (lines that can run)29
Code lines that did run29
Code lines that did not run0
Coverage (did run/can run)100.00 %
Function listing
   time   calls  line
1 function [Gsubmatrix] = integrateGsubmatrix_Telles(ngp, elcoords, bsFnConn, collocCoords, srcXi_param, range)
2 % integrate the G submatrix using Telles' transformation
3
0.01 31 4 global controlPts p knotVec
5
0.04 31 6 [gpt gwt]=lgwt(ngp,-1,1);
7
31 8 numBasisFns=length(bsFnConn);
0.01 31 9 N=zeros(1,numBasisFns); dN=zeros(1,numBasisFns);
31 10 Gsubmatrix=zeros(2,6);
11
31 12 srcXi=convertToParentCoordSpace(srcXi_param, range);
31 13 xiStar=srcXi^2 - 1;
0.02 31 14 gamBar=nthroot( srcXi*xiStar + abs(xiStar) ,3) + nthroot(srcXi*xiStar - abs(xiStar),3) + srcXi;
31 15 jacob_param=(range(2)-range(1)) / 2; % jacobian from parent to parameter space
16
31 17 for pt=1:ngp % integrate using Gaussian quadrature
186 18 gXi=gpt(pt);
186 19 xi=((gXi-gamBar).^3 + gamBar*(gamBar.^2+3)) / (1+3*gamBar.^2);
186 20 xi_param=convertToParamSpace( xi, range); % get the gauss point in parameter space
21
0.01 186 22 for lclBasis=1:numBasisFns
558 23 i=bsFnConn(lclBasis);
558 24 [N(lclBasis) dN(lclBasis)]=NURBSbasis(i, p, xi_param, knotVec, controlPts(:,3)' );
558 25 end
26
0.02 186 27 [jacob_xi, normals, r, dr, drdn] = getKernelParameters( elcoords, collocCoords, N, dN );
186 28 jacob=jacob_xi*jacob_param; % the final jacobian we use
29
0.01 186 30 Ttemp=zeros(2,2); Utemp=zeros(2,2);
186 31 for i=1:2
372 32 for j=1:2
0.12 744 33 [Ttemp(i,j) Utemp(i,j)]=DBIEkernels(i,j,r,dr,drdn,normals);
0.01 744 34 end
0.01 372 35 end
186 36 jacobTelles=3* ((gXi-gamBar).^2) / (1+3*gamBar.^2);
186 37 Gsubmatrix=Gsubmatrix + [N(1)*Utemp N(2)*Utemp N(3)*Utemp]*jacob*gwt(pt)*jacobTelles;
38
186 39 end
31 40 end

Other subfunctions in this file are not included in this listing.