This is a static copy of a profile reportHome
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 Number | Code | Calls | Total Time | % Time | Time Plot |
33 | [Ttemp(i,j) Utemp(i,j)]=DBIEke... | 744 | 0.124 s | 44.4% |  |
6 | [gpt gwt]=lgwt(ngp,-1,1); | 31 | 0.041 s | 14.8% |  |
27 | [jacob_xi, normals, r, dr, drd... | 186 | 0.021 s | 7.4% |  |
14 | gamBar=nthroot( srcXi*xiStar +... | 31 | 0.021 s | 7.4% |  |
35 | end | 372 | 0.010 s | 3.7% |  |
All other lines | | | 0.062 s | 22.2% |  |
Totals | | | 0.280 s | 100% | |
Children (called functions)
M-Lint results
No M-Lint messages.Coverage results
[ Show coverage for parent directory ]
Total lines in function | 40 |
Non-code lines (comments, blank lines) | 11 |
Code lines (lines that can run) | 29 |
Code lines that did run | 29 |
Code lines that did not run | 0 |
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.