This is a static copy of a profile report

Home

integrateHGsubmatrices_GLQ (255 calls, 1.960 sec)
Generated 18-Mar-2011 23:31:44 using cpu time.
M-function in file /Users/Robert/Documents/Work/Cardiff/postdoc/myrepo/NURBS/isoBEM/integrateHGsubmatrices_GLQ.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
6
[gpt gwt]=lgwt(ngp,-1,1);
2550.705 s36.0%
31
[Ttemp(i,j) Utemp(i,j)]=DBIEke...
61200.653 s33.3%
22
[N(lclBasis) dN(lclBasis)]=NUR...
45900.238 s12.2%
25
[jacob_xi, normals, r, dr, drd...
15300.156 s7.9%
36
Gsubmatrix=Gsubmatrix + [N(1)*...
15300.031 s1.6%
All other lines  0.176 s9.0%
Totals  1.960 s100% 
Children (called functions)

Function NameFunction TypeCallsTotal Time% TimeTime Plot
lgwtM-function2550.695 s35.4%
DBIEkernelsM-function61200.529 s27.0%
getKernelParametersM-function15300.135 s6.9%
NURBSbasisMEX-function45900.083 s4.2%
convertToParamSpaceM-function15300 s0%
Self time (built-ins, overhead, etc.)  0.518 s26.5%
Totals  1.960 s100% 
M-Lint results
No M-Lint messages.
Coverage results
[ Show coverage for parent directory ]
Total lines in function39
Non-code lines (comments, blank lines)14
Code lines (lines that can run)25
Code lines that did run25
Code lines that did not run0
Coverage (did run/can run)100.00 %
Function listing
   time   calls  line
1 function [ Hsubmatrix, Gsubmatrix ] = integrateHGsubmatrices_GLQ( ngp, elcoords, bsFnConn, collocCoords, range)
2 % integrate the H and G submatrices using Gauss Legendre quadrature
3
0.02 255 4 global controlPts p knotVec
5
0.71 255 6 [gpt gwt]=lgwt(ngp,-1,1);
7
255 8 Hsubmatrix=zeros(2,6);
255 9 Gsubmatrix=zeros(2,6);
10
255 11 numBasisFns=length(bsFnConn);
255 12 N=zeros(1,numBasisFns); dN=zeros(1,numBasisFns);
13
255 14 jacob_param=(range(2)-range(1)) / 2; % jacobian from parent to parameter space
15
255 16 for pt=1:ngp % integrate using Gaussian quadrature
17
0.02 1530 18 xi_param=convertToParamSpace( gpt(pt), range); % get the gauss point in parameter space
19
1530 20 for lclBasis=1:numBasisFns
4590 21 i=bsFnConn(lclBasis);
0.24 4590 22 [N(lclBasis) dN(lclBasis)]=NURBSbasis(i, p, xi_param, knotVec, controlPts(:,3)' );
0.01 4590 23 end
24
0.16 1530 25 [jacob_xi, normals, r, dr, drdn] = getKernelParameters( elcoords, collocCoords, N, dN );
1530 26 jacob=jacob_xi*jacob_param; % the final jacobian we use
27
1530 28 Ttemp=zeros(2,2); Utemp=zeros(2,2);
0.01 1530 29 for i=1:2
0.02 3060 30 for j=1:2
0.65 6120 31 [Ttemp(i,j) Utemp(i,j)]=DBIEkernels(i,j,r,dr,drdn,normals);
0.03 6120 32 end
0.03 3060 33 end
34
1530 35 Hsubmatrix=Hsubmatrix + [N(1)*Ttemp N(2)*Ttemp N(3)*Ttemp]*jacob*gwt(pt);
0.03 1530 36 Gsubmatrix=Gsubmatrix + [N(1)*Utemp N(2)*Utemp N(3)*Utemp]*jacob*gwt(pt);
37
0.02 1530 38 end
255 39 end

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