<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Hi Rafael Guerra, et al.<br>
<br>
Thank you for clarifying this to me - that interp1 uses not-a-knot
as default. IMHO Scilab documentation should be clear about this.<font
face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif"><br>
<br>
>What do you mean by your result is “OK, but it's not
perfect, compared to a Fortran script</span></font><font
face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">”?</span></font><br>
<br>
I can try to make a snapshot ... what kind do you need? Like a
difference in the resulting vector, a piece of source code ... or
maybe a graph plot illustrating the difference?<br>
<br>
Since the default acc. to your script is not-a-knot ... returning
to my original question, is it true I can alter this to 'natural'
cubic spline by writing:<font face="Times New Roman" color="black"
size="3"><span style="font-size:12.0pt"><br>
<br>
a</span></font><font color="#4a55db"><span
style="color:#4A55DB">(</span></font><font color="#ffaa00"><span
style="color:#FFAA00">:</span></font>,<font color="rosybrown"><span
style="color:rosybrown">1</span></font><font color="#4a55db"><span
style="color:#4A55DB">)</span></font>
<font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font>
<u><font color="#ae5cb0"><span style="color:#AE5CB0">interp1</span></font></u><font
color="#4a55db"><span style="color:#4A55DB">(</span></font><font
color="#32b9b9"><span style="color:#32B9B9">log</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font>f3<font
color="#4a55db"><span style="color:#4A55DB">)</span></font>,a1,<font
color="#32b9b9"><span style="color:#32B9B9">log</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font>f<font
color="#4a55db"><span style="color:#4A55DB">)</span></font>,<font
color="rosybrown"><span style="color:rosybrown">'<span
class="SpellE">spline'<font color="black"><span
style="color:black">,</span></font>'n</span>atural'</span></font><font
color="#4a55db"><span style="color:#4A55DB">)</span></font>;<o:p></o:p><br>
<br>
<font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">?</span></font><br>
<br>
... The documentation for interp1 spline isn't clear whether
'natural' is an option for the extrapolation. The documentation in
Scilab ( <a class="moz-txt-link-freetext" href="https://help.scilab.org/docs/6.0.0/en_US/interp1.html">https://help.scilab.org/docs/6.0.0/en_US/interp1.html</a> )
only says:<br>
<br>
<p class="para" style="color: rgb(0, 0, 0); font-family:
sans-serif; font-size: 16px; font-style: normal;
font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: normal; letter-spacing: normal; orphans: 2;
text-align: start; text-indent: 0px; text-transform: none;
white-space: normal; widows: 2; word-spacing: 0px;
-webkit-text-stroke-width: 0px;">The<span
class="Apple-converted-space"> </span><code class="literal"
style="font-size: 13px;">extrapolation</code><span
class="Apple-converted-space"> </span>parameter sets the
evaluation rule for extrapolation, i.e for<span
class="Apple-converted-space"> </span><code class="literal"
style="font-size: 13px;">xp(i)</code>not in [x1,xn] interval</p>
<dl style="color: rgb(0, 0, 0); font-family: sans-serif;
font-size: 16px; font-style: normal; font-variant-ligatures:
normal; font-variant-caps: normal; font-weight: normal;
letter-spacing: normal; orphans: 2; text-align: start;
text-indent: 0px; text-transform: none; white-space: normal;
widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px;">
<dt><span class="term" style="color: rgb(128, 0, 0); font-size:
16px;">"extrap"</span></dt>
<dd>
<p class="para">the extrapolation is performed by the defined
method. yp=interp1(x,y,xp,method,"extrap")</p>
</dd>
</dl>
<br>
... And the documentation provides no examples of adding this
"extrap" parameter to the spline fitting (or any other fitting).<br>
<br>
<br>
???<br>
<br>
Anyway, the Fortran code uses 'natural' cubic spline and (as
mentioned) <font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt"><span class="SpellE">DGTTRF</span>
and <span class="SpellE">
DGTTRS </span></span></font>- are there any Scilab
equivalent functions? ... I might wish to try to code (exactly)
the same spline functionality, just to eliminate a potential error
here.<br>
<br>
DGTTRF computes an LU factorization of a real tridiagonal matrix A
<br>
* using elimination with partial pivoting and row interchanges.<br>
<br>
DGTTRS solves one of the systems of equations<br>
* A*X = B or A'*X = B,<br>
* with a tridiagonal matrix A using the LU factorization computed<br>
* by DGTTRF.<br>
<br>
( Source:
<a class="moz-txt-link-freetext" href="https://www.gfd-dennou.org/arch/ruby/products/ruby-lapack/doc/dgt.html">https://www.gfd-dennou.org/arch/ruby/products/ruby-lapack/doc/dgt.html</a>
)<br>
<br>
P.S. Sorry for the not-so-fast response time, 1) I'm thinking a
lot about it, 2) I shall try not to make too many false statements
... but it's a bit difficult because I don't understand why I get
the differences that I observe - and I feel a bit like I'm
searching for the reason (the solution) while being blind-folded.<br>
<br>
Best regards,<br>
Claus<br>
<br>
On 11-09-2016 18:48, Rafael Guerra wrote:<br>
</div>
<blockquote
cite="mid:DB5PR04MB1575A0BE875B5956644D6175CCFC0@DB5PR04MB1575.eurprd04.prod.outlook.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<meta name="ProgId" content="Word.Document">
<meta name="Generator" content="Microsoft Word 15">
<meta name="Originator" content="Microsoft Word 15">
<link rel="File-List" href="cid:filelist.xml@01D20C5D.1413F5F0">
<!--[if gte mso 9]><xml>
<o:OfficeDocumentSettings>
<o:AllowPNG/>
<o:DoNotRelyOnCSS/>
</o:OfficeDocumentSettings>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:WordDocument>
<w:Zoom>180</w:Zoom>
<w:SpellingState>Clean</w:SpellingState>
<w:TrackMoves/>
<w:TrackFormatting/>
<w:EnvelopeVis/>
<w:ValidateAgainstSchemas/>
<w:SaveIfXMLInvalid>false</w:SaveIfXMLInvalid>
<w:IgnoreMixedContent>false</w:IgnoreMixedContent>
<w:AlwaysShowPlaceholderText>false</w:AlwaysShowPlaceholderText>
<w:DoNotPromoteQF/>
<w:LidThemeOther>EN-US</w:LidThemeOther>
<w:LidThemeAsian>X-NONE</w:LidThemeAsian>
<w:LidThemeComplexScript>X-NONE</w:LidThemeComplexScript>
<w:Compatibility>
<w:DoNotExpandShiftReturn/>
<w:BreakWrappedTables/>
<w:SplitPgBreakAndParaMark/>
<w:EnableOpenTypeKerning/>
</w:Compatibility>
<m:mathPr>
<m:mathFont m:val="Cambria Math"/>
<m:brkBin m:val="before"/>
<m:brkBinSub m:val="--"/>
<m:smallFrac m:val="off"/>
<m:dispDef/>
<m:lMargin m:val="0"/>
<m:rMargin m:val="0"/>
<m:defJc m:val="centerGroup"/>
<m:wrapIndent m:val="1440"/>
<m:intLim m:val="subSup"/>
<m:naryLim m:val="undOvr"/>
</m:mathPr></w:WordDocument>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:LatentStyles DefLockedState="false" DefUnhideWhenUsed="false" DefSemiHidden="false" DefQFormat="false" DefPriority="99" LatentStyleCount="371">
<w:LsdException Locked="false" Priority="0" QFormat="true" Name="Normal"/>
<w:LsdException Locked="false" Priority="9" QFormat="true" Name="heading 1"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 2"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 3"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 4"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 5"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 6"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 7"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 8"/>
<w:LsdException Locked="false" Priority="9" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="heading 9"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 5"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 6"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 7"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 8"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index 9"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 1"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 2"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 3"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 4"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 5"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 6"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 7"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 8"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" Name="toc 9"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Normal Indent"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="footnote text"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="annotation text"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="header"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="footer"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="index heading"/>
<w:LsdException Locked="false" Priority="35" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="caption"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="table of figures"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="envelope address"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="envelope return"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="footnote reference"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="annotation reference"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="line number"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="page number"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="endnote reference"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="endnote text"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="table of authorities"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="macro"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="toa heading"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Bullet"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Number"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List 5"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Bullet 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Bullet 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Bullet 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Bullet 5"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Number 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Number 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Number 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Number 5"/>
<w:LsdException Locked="false" Priority="10" QFormat="true" Name="Title"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Closing"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Signature"/>
<w:LsdException Locked="false" Priority="1" SemiHidden="true" UnhideWhenUsed="true" Name="Default Paragraph Font"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text Indent"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Continue"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Continue 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Continue 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Continue 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="List Continue 5"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Message Header"/>
<w:LsdException Locked="false" Priority="11" QFormat="true" Name="Subtitle"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Salutation"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Date"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text First Indent"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text First Indent 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Note Heading"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text Indent 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Body Text Indent 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Block Text"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Hyperlink"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="FollowedHyperlink"/>
<w:LsdException Locked="false" Priority="22" QFormat="true" Name="Strong"/>
<w:LsdException Locked="false" Priority="20" QFormat="true" Name="Emphasis"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Document Map"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Plain Text"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="E-mail Signature"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Top of Form"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Bottom of Form"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Normal (Web)"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Acronym"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Address"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Cite"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Code"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Definition"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Keyboard"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Preformatted"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Sample"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Typewriter"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="HTML Variable"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Normal Table"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="annotation subject"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="No List"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Outline List 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Outline List 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Outline List 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Simple 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Simple 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Simple 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Classic 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Classic 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Classic 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Classic 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Colorful 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Colorful 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Colorful 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Columns 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Columns 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Columns 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Columns 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Columns 5"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 5"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 6"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 7"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Grid 8"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 4"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 5"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 6"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 7"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table List 8"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table 3D effects 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table 3D effects 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table 3D effects 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Contemporary"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Elegant"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Professional"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Subtle 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Subtle 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Web 1"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Web 2"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Web 3"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Balloon Text"/>
<w:LsdException Locked="false" Priority="39" Name="Table Grid"/>
<w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true" Name="Table Theme"/>
<w:LsdException Locked="false" SemiHidden="true" Name="Placeholder Text"/>
<w:LsdException Locked="false" Priority="1" QFormat="true" Name="No Spacing"/>
<w:LsdException Locked="false" Priority="60" Name="Light Shading"/>
<w:LsdException Locked="false" Priority="61" Name="Light List"/>
<w:LsdException Locked="false" Priority="62" Name="Light Grid"/>
<w:LsdException Locked="false" Priority="63" Name="Medium Shading 1"/>
<w:LsdException Locked="false" Priority="64" Name="Medium Shading 2"/>
<w:LsdException Locked="false" Priority="65" Name="Medium List 1"/>
<w:LsdException Locked="false" Priority="66" Name="Medium List 2"/>
<w:LsdException Locked="false" Priority="67" Name="Medium Grid 1"/>
<w:LsdException Locked="false" Priority="68" Name="Medium Grid 2"/>
<w:LsdException Locked="false" Priority="69" Name="Medium Grid 3"/>
<w:LsdException Locked="false" Priority="70" Name="Dark List"/>
<w:LsdException Locked="false" Priority="71" Name="Colorful Shading"/>
<w:LsdException Locked="false" Priority="72" Name="Colorful List"/>
<w:LsdException Locked="false" Priority="73" Name="Colorful Grid"/>
<w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 1"/>
<w:LsdException Locked="false" Priority="61" Name="Light List Accent 1"/>
<w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 1"/>
<w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 1"/>
<w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 1"/>
<w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 1"/>
<w:LsdException Locked="false" SemiHidden="true" Name="Revision"/>
<w:LsdException Locked="false" Priority="34" QFormat="true" Name="List Paragraph"/>
<w:LsdException Locked="false" Priority="29" QFormat="true" Name="Quote"/>
<w:LsdException Locked="false" Priority="30" QFormat="true" Name="Intense Quote"/>
<w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 1"/>
<w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 1"/>
<w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 1"/>
<w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 1"/>
<w:LsdException Locked="false" Priority="70" Name="Dark List Accent 1"/>
<w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 1"/>
<w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 1"/>
<w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 1"/>
<w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 2"/>
<w:LsdException Locked="false" Priority="61" Name="Light List Accent 2"/>
<w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 2"/>
<w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 2"/>
<w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 2"/>
<w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 2"/>
<w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 2"/>
<w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 2"/>
<w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 2"/>
<w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 2"/>
<w:LsdException Locked="false" Priority="70" Name="Dark List Accent 2"/>
<w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 2"/>
<w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 2"/>
<w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 2"/>
<w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 3"/>
<w:LsdException Locked="false" Priority="61" Name="Light List Accent 3"/>
<w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 3"/>
<w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 3"/>
<w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 3"/>
<w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 3"/>
<w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 3"/>
<w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 3"/>
<w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 3"/>
<w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 3"/>
<w:LsdException Locked="false" Priority="70" Name="Dark List Accent 3"/>
<w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 3"/>
<w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 3"/>
<w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 3"/>
<w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 4"/>
<w:LsdException Locked="false" Priority="61" Name="Light List Accent 4"/>
<w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 4"/>
<w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 4"/>
<w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 4"/>
<w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 4"/>
<w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 4"/>
<w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 4"/>
<w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 4"/>
<w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 4"/>
<w:LsdException Locked="false" Priority="70" Name="Dark List Accent 4"/>
<w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 4"/>
<w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 4"/>
<w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 4"/>
<w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 5"/>
<w:LsdException Locked="false" Priority="61" Name="Light List Accent 5"/>
<w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 5"/>
<w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 5"/>
<w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 5"/>
<w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 5"/>
<w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 5"/>
<w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 5"/>
<w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 5"/>
<w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 5"/>
<w:LsdException Locked="false" Priority="70" Name="Dark List Accent 5"/>
<w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 5"/>
<w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 5"/>
<w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 5"/>
<w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 6"/>
<w:LsdException Locked="false" Priority="61" Name="Light List Accent 6"/>
<w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 6"/>
<w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 6"/>
<w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 6"/>
<w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 6"/>
<w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 6"/>
<w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 6"/>
<w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 6"/>
<w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 6"/>
<w:LsdException Locked="false" Priority="70" Name="Dark List Accent 6"/>
<w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 6"/>
<w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 6"/>
<w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 6"/>
<w:LsdException Locked="false" Priority="19" QFormat="true" Name="Subtle Emphasis"/>
<w:LsdException Locked="false" Priority="21" QFormat="true" Name="Intense Emphasis"/>
<w:LsdException Locked="false" Priority="31" QFormat="true" Name="Subtle Reference"/>
<w:LsdException Locked="false" Priority="32" QFormat="true" Name="Intense Reference"/>
<w:LsdException Locked="false" Priority="33" QFormat="true" Name="Book Title"/>
<w:LsdException Locked="false" Priority="37" SemiHidden="true" UnhideWhenUsed="true" Name="Bibliography"/>
<w:LsdException Locked="false" Priority="39" SemiHidden="true" UnhideWhenUsed="true" QFormat="true" Name="TOC Heading"/>
<w:LsdException Locked="false" Priority="41" Name="Plain Table 1"/>
<w:LsdException Locked="false" Priority="42" Name="Plain Table 2"/>
<w:LsdException Locked="false" Priority="43" Name="Plain Table 3"/>
<w:LsdException Locked="false" Priority="44" Name="Plain Table 4"/>
<w:LsdException Locked="false" Priority="45" Name="Plain Table 5"/>
<w:LsdException Locked="false" Priority="40" Name="Grid Table Light"/>
<w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light"/>
<w:LsdException Locked="false" Priority="47" Name="Grid Table 2"/>
<w:LsdException Locked="false" Priority="48" Name="Grid Table 3"/>
<w:LsdException Locked="false" Priority="49" Name="Grid Table 4"/>
<w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark"/>
<w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful"/>
<w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful"/>
<w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light Accent 1"/>
<w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 1"/>
<w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 1"/>
<w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 1"/>
<w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 1"/>
<w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful Accent 1"/>
<w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful Accent 1"/>
<w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light Accent 2"/>
<w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 2"/>
<w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 2"/>
<w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 2"/>
<w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 2"/>
<w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful Accent 2"/>
<w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful Accent 2"/>
<w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light Accent 3"/>
<w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 3"/>
<w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 3"/>
<w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 3"/>
<w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 3"/>
<w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful Accent 3"/>
<w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful Accent 3"/>
<w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light Accent 4"/>
<w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 4"/>
<w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 4"/>
<w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 4"/>
<w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 4"/>
<w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful Accent 4"/>
<w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful Accent 4"/>
<w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light Accent 5"/>
<w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 5"/>
<w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 5"/>
<w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 5"/>
<w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 5"/>
<w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful Accent 5"/>
<w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful Accent 5"/>
<w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light Accent 6"/>
<w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 6"/>
<w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 6"/>
<w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 6"/>
<w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 6"/>
<w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful Accent 6"/>
<w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful Accent 6"/>
<w:LsdException Locked="false" Priority="46" Name="List Table 1 Light"/>
<w:LsdException Locked="false" Priority="47" Name="List Table 2"/>
<w:LsdException Locked="false" Priority="48" Name="List Table 3"/>
<w:LsdException Locked="false" Priority="49" Name="List Table 4"/>
<w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark"/>
<w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful"/>
<w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful"/>
<w:LsdException Locked="false" Priority="46" Name="List Table 1 Light Accent 1"/>
<w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 1"/>
<w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 1"/>
<w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 1"/>
<w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 1"/>
<w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful Accent 1"/>
<w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful Accent 1"/>
<w:LsdException Locked="false" Priority="46" Name="List Table 1 Light Accent 2"/>
<w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 2"/>
<w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 2"/>
<w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 2"/>
<w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 2"/>
<w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful Accent 2"/>
<w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful Accent 2"/>
<w:LsdException Locked="false" Priority="46" Name="List Table 1 Light Accent 3"/>
<w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 3"/>
<w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 3"/>
<w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 3"/>
<w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 3"/>
<w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful Accent 3"/>
<w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful Accent 3"/>
<w:LsdException Locked="false" Priority="46" Name="List Table 1 Light Accent 4"/>
<w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 4"/>
<w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 4"/>
<w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 4"/>
<w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 4"/>
<w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful Accent 4"/>
<w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful Accent 4"/>
<w:LsdException Locked="false" Priority="46" Name="List Table 1 Light Accent 5"/>
<w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 5"/>
<w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 5"/>
<w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 5"/>
<w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 5"/>
<w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful Accent 5"/>
<w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful Accent 5"/>
<w:LsdException Locked="false" Priority="46" Name="List Table 1 Light Accent 6"/>
<w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 6"/>
<w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 6"/>
<w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 6"/>
<w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 6"/>
<w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful Accent 6"/>
<w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful Accent 6"/>
</w:LatentStyles>
</xml><![endif]-->
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;
mso-font-charset:1;
mso-generic-font-family:roman;
mso-font-pitch:variable;
mso-font-signature:-536870145 1107305727 0 0 415 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;
mso-font-charset:0;
mso-generic-font-family:swiss;
mso-font-pitch:variable;
mso-font-signature:-536870145 1073786111 1 0 415 0;}
@font-face
{font-family:Monospaced;
panose-1:0 0 0 0 0 0 0 0 0 0;
mso-font-alt:"Times New Roman";
mso-font-charset:0;
mso-generic-font-family:roman;
mso-font-format:other;
mso-font-pitch:auto;
mso-font-signature:0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{mso-style-unhide:no;
mso-style-qformat:yes;
mso-style-parent:"";
margin:0in;
margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:12.0pt;
font-family:"Times New Roman",serif;
mso-fareast-font-family:Calibri;
color:black;}
a:link, span.MsoHyperlink
{mso-style-noshow:yes;
mso-style-priority:99;
color:#0563C1;
text-decoration:underline;
text-underline:single;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-noshow:yes;
mso-style-priority:99;
color:#954F72;
text-decoration:underline;
text-underline:single;}
p
{mso-style-noshow:yes;
mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
mso-pagination:widow-orphan;
font-size:12.0pt;
font-family:"Times New Roman",serif;
mso-fareast-font-family:Calibri;
color:black;}
pre
{mso-style-noshow:yes;
mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0in;
margin-bottom:.0001pt;
mso-pagination:widow-orphan;
tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt;
font-size:10.0pt;
font-family:"Courier New";
mso-fareast-font-family:"Times New Roman";
color:windowtext;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-style-unhide:no;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
mso-pagination:widow-orphan;
font-size:12.0pt;
font-family:"Times New Roman",serif;
mso-fareast-font-family:Calibri;
color:black;}
span.EmailStyle19
{mso-style-type:personal-reply;
mso-style-noshow:yes;
mso-style-unhide:no;
mso-ansi-font-size:11.0pt;
mso-bidi-font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-ascii-font-family:Calibri;
mso-hansi-font-family:Calibri;
color:black;
font-weight:normal;
font-style:normal;
text-decoration:none;
text-underline:none;
text-decoration:none;
text-line-through:none;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-noshow:yes;
mso-style-priority:99;
mso-style-unhide:no;
mso-style-locked:yes;
mso-style-link:"HTML Preformatted";
font-family:"Courier New";
mso-ascii-font-family:"Courier New";
mso-hansi-font-family:"Courier New";
mso-bidi-font-family:"Courier New";}
span.SpellE
{mso-style-name:"";
mso-spl-e:yes;}
.MsoChpDefault
{mso-style-type:export-only;
mso-default-props:yes;
font-size:10.0pt;
mso-ansi-font-size:10.0pt;
mso-bidi-font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;
mso-header-margin:.5in;
mso-footer-margin:.5in;
mso-paper-source:0;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 10]><style>/* Style Definitions */
table.MsoNormalTable
{mso-style-name:"Table Normal";
mso-tstyle-rowband-size:0;
mso-tstyle-colband-size:0;
mso-style-noshow:yes;
mso-style-priority:99;
mso-style-parent:"";
mso-padding-alt:0in 5.4pt 0in 5.4pt;
mso-para-margin:0in;
mso-para-margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:10.0pt;
font-family:"Times New Roman",serif;}
</style><![endif]--><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">Hi,<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif"><o:p> </o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">The
code line using
<b style="mso-bidi-font-weight:normal"><span
style="font-weight:bold;mso-bidi-font-weight:normal">interp1</span></b>
and the spline method seems to perform "<span
class="SpellE">not_a_knot</span>" cubic spline
interpolation, as demonstrated by the simple test here
below:<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></font></p>
<pre><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">x0</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <font color="#4a55db"><span style="color:#4A55DB">[</span></font><font color="rosybrown"><span style="color:rosybrown">0</span></font> <font color="rosybrown"><span style="color:rosybrown">1</span></font> <font color="rosybrown"><span style="color:rosybrown">2</span></font> <font color="rosybrown"><span style="color:rosybrown">3</span></font> <font color="rosybrown"><span style="color:rosybrown">4</span></font><font color="#4a55db"><span style="color:#4A55DB">]</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">y0</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <font color="#4a55db"><span style="color:#4A55DB">[</span></font><font color="rosybrown"><span style="color:rosybrown">0</span></font> <font color="#5c5c5c"><span style="color:#5C5C5C">-</span></font><font color="rosybrown"><span style="color:rosybrown">1</span></font> <font color="rosybrown"><span style="color:rosybrown">0</span></font> <font color="rosybrown"><span style="color:rosybrown">2</span></font> <font color="rosybrown"><span style="color:rosybrown">1</span></font><font color="#4a55db"><span style="color:#4A55DB">]</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">x</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <font color="rosybrown"><span style="color:rosybrown">0</span></font><font color="#ffaa00"><span style="color:#FFAA00">:</span></font><font color="rosybrown"><span style="color:rosybrown">0.05</span></font><font color="#ffaa00"><span style="color:#FFAA00">:</span></font><font color="rosybrown"><span style="color:rosybrown">4</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">y1</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif;mso-ansi-language:ES" lang="ES"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <u><font color="#ae5cb0"><span style="color:#AE5CB0">interp1</span></font></u><font color="#4a55db"><span style="color:#4A55DB">(</span></font><font color="black"><span style="color:black">x0,y0,x,</span></font><font color="rosybrown"><span style="color:rosybrown">'spline'</span></font><font color="#4a55db"><span style="color:#4A55DB">)</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">d</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <span class="SpellE"><font color="#32b9b9"><span style="color:#32B9B9">splin</span></font></span><font color="#4a55db"><span style="color:#4A55DB">(</span></font><font color="black"><span style="color:black">x0,</span></font> <font color="black"><span style="color:black">y0,</span></font><font color="rosybrown"><span style="color:rosybrown">"not_a_knot"</span></font><font color="#4a55db"><span style="color:#4A55DB">)</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><span class="SpellE"><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">e</span></font></span><font face="Monospaced"><span style="font-family:"Monospaced",serif;mso-ansi-language:ES" lang="ES"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <span class="SpellE"><font color="#32b9b9"><span style="color:#32B9B9">splin</span></font></span><font color="#4a55db"><span style="color:#4A55DB">(</span></font><font color="black"><span style="color:black">x0,</span></font> <font color="black"><span style="color:black">y0,</span></font><font color="rosybrown"><span style="color:rosybrown">"natural"</span></font><font color="#4a55db"><span style="color:#4A55DB">)</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">y2</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif;mso-ansi-language:ES" lang="ES"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <span class="SpellE"><font color="#32b9b9"><span style="color:#32B9B9">interp</span></font></span><font color="#4a55db"><span style="color:#4A55DB">(</span></font><font color="black"><span style="color:black">x,</span></font> <font color="black"><span style="color:black">x0,</span></font> <font color="black"><span style="color:black">y0,</span></font> <font color="black"><span style="color:black">d</span></font><font color="#4a55db"><span style="color:#4A55DB">)</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><font face="Monospaced" color="black" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">y3</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif;mso-ansi-language:ES" lang="ES"> <font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font> <span class="SpellE"><font color="#32b9b9"><span style="color:#32B9B9">interp</span></font></span><font color="#4a55db"><span style="color:#4A55DB">(</span></font><font color="black"><span style="color:black">x,</span></font> <font color="black"><span style="color:black">x0,</span></font> <font color="black"><span style="color:black">y0,</span></font> <font color="black"><span style="color:black">e</span></font><font color="#4a55db"><span style="color:#4A55DB">)</span></font><font color="black"><span style="color:black">;</span></font><o:p></o:p></span></font></pre>
<pre><span class="SpellE"><u><font face="Monospaced" color="#ae5cb0" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:#AE5CB0;mso-ansi-language:ES" lang="ES">clf</span></font></u></span><font face="Monospaced" color="#4a55db"><span style="font-family:"Monospaced",serif;color:#4A55DB;mso-ansi-language:ES" lang="ES">()</span></font><font face="Monospaced" color="black"><span style="font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif;mso-ansi-language:ES" lang="ES"><o:p></o:p></span></font></pre>
<pre><span class="SpellE"><u><font face="Monospaced" color="#ae5cb0" size="2"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:#AE5CB0;mso-ansi-language:ES" lang="ES">plot</span></font></u></span><font face="Monospaced" color="#4a55db"><span style="font-family:"Monospaced",serif;color:#4A55DB;mso-ansi-language:ES" lang="ES">(</span></font><font face="Monospaced" color="black"><span style="font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">x,y1,</span></font><font face="Monospaced" color="rosybrown"><span style="font-family:"Monospaced",serif;color:rosybrown;mso-ansi-language:ES" lang="ES">'blue'</span></font><font face="Monospaced" color="black"><span style="font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">,x,y2,</span></font><font face="Monospaced" color="rosybrown"><span style="font-family:"Monospaced",serif;color:rosybrown;mso-ansi-language:ES" lang="ES">'--green'</span></font><font face="Monospaced" color="black"><span style="font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">,x,y3,</span></font><font face="Monospaced" color="rosybrown"><span style="font-family:"Monospaced",serif;color:rosybrown;mso-ansi-language:ES" lang="ES">'red'</span></font><font face="Monospaced" color="black"><span style="font-family:"Monospaced",serif;color:black;mso-ansi-language:ES" lang="ES">,x0,y0,</span></font><font face="Monospaced" color="rosybrown"><span style="font-family:"Monospaced",serif;color:rosybrown;mso-ansi-language:ES" lang="ES">'Xblack'</span></font><font face="Monospaced" color="#4a55db"><span style="font-family:"Monospaced",serif;color:#4A55DB;mso-ansi-language:ES" lang="ES">)</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif;mso-ansi-language:ES" lang="ES"><o:p></o:p></span></font></pre>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-ansi-language:ES"
lang="ES"><o:p> </o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-ansi-language:ES"
lang="ES"><o:p> </o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif">What
do you mean by your result is “OK, but it's not perfect,
compared to a Fortran script</span></font><font
face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">”?<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">Could
you provide a snapshot?<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">Regards,<o:p></o:p></span></font></p>
<p class="MsoNormal"><span class="SpellE"><font face="Calibri"
color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">Rafae</span></font></span><font
face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">;<o:p></o:p></span></font></p>
<p class="MsoNormal"><font face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></font></p>
<div>
<div style="border:none;border-top:solid #E1E1E1
1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><font face="Calibri" color="black"
size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-font-family:"Times
New Roman";color:windowtext;font-weight:bold">From:</span></font></b><font
face="Calibri" color="black" size="2"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-font-family:"Times
New Roman";color:windowtext"> Claus <span
class="SpellE">Futtrup</span><br>
<b><span style="font-weight:bold">Sent:</span></b>
Friday, September 09, 2016 8:34 PM<br>
<b><span style="font-weight:bold">To:</span></b>
International users mailing list for Scilab.
<a class="moz-txt-link-rfc2396E" href="mailto:users@lists.scilab.org"><users@lists.scilab.org></a><br>
<b><span style="font-weight:bold">Subject:</span></b>
[Scilab-users] Cubic spline<o:p></o:p></span></font></p>
</div>
</div>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">Hi there<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">In Scilab I've used interp1 to
calculate a cubic spline interpolation, like this:<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">a</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font><font
color="#ffaa00"><span style="color:#FFAA00">:</span></font>,<font
color="rosybrown"><span style="color:rosybrown">1</span></font><font
color="#4a55db"><span style="color:#4A55DB">)</span></font>
<font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font>
<u><font color="#ae5cb0"><span style="color:#AE5CB0">interp1</span></font></u><font
color="#4a55db"><span style="color:#4A55DB">(</span></font><font
color="#32b9b9"><span style="color:#32B9B9">log</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font>f3<font
color="#4a55db"><span style="color:#4A55DB">)</span></font>,a1,<font
color="#32b9b9"><span style="color:#32B9B9">log</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font>f<font
color="#4a55db"><span style="color:#4A55DB">)</span></font>,<font
color="rosybrown"><span style="color:rosybrown">'spline'</span></font><font
color="#4a55db"><span style="color:#4A55DB">)</span></font>;<o:p></o:p></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">a = amplitude (magnitude). f3 is
a frequency (27000 linear spaced data), a1 is the original
data, f is the resampled 1200 frequencies (log-spaced),
where I need the spline to interpolate some data for me.<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">Above should work OK, but it's
not perfect, compared to a Fortran script. The
<span class="SpellE">fortran</span> script calculates with
its own cubic spline routine, utilizing
<span class="SpellE">LAPACK</span> (<span class="SpellE">DGTTRF</span>
and <span class="SpellE">
DGTTRS</span>) to solve for polynomial coefficients.</span></font><font
color="black"><span style="color:black"><o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">The question is - above code line
with the interp1 spline, which kind of spline is it?<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">Digging into
<span class="SpellE">interp</span>, I see multiple
options. Digging into <span class="SpellE">
splin</span>, I also see multiple options.<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">I looks like the interp1 is using
"natural" spline - is this correct?<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">It's strange because the
<span class="SpellE">splin</span> help documentation
doesn't recommend this. It says: Don't use the natural
type unless the underlying function have zero second end
points derivatives.<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">This might be my problem.<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">The Scilab help for interp1
doesn't give any examples, but does mention I can add an "<span
class="SpellE">extrap</span>" method. Could this be any
of the suggestions in the
<span class="SpellE">splin</span> documentation. For
example, could I write:<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">a</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font><font
color="#ffaa00"><span style="color:#FFAA00">:</span></font>,<font
color="rosybrown"><span style="color:rosybrown">1</span></font><font
color="#4a55db"><span style="color:#4A55DB">)</span></font>
<font color="#5c5c5c"><span style="color:#5C5C5C">=</span></font>
<u><font color="#ae5cb0"><span style="color:#AE5CB0">interp1</span></font></u><font
color="#4a55db"><span style="color:#4A55DB">(</span></font><font
color="#32b9b9"><span style="color:#32B9B9">log</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font>f3<font
color="#4a55db"><span style="color:#4A55DB">)</span></font>,a1,<font
color="#32b9b9"><span style="color:#32B9B9">log</span></font><font
color="#4a55db"><span style="color:#4A55DB">(</span></font>f<font
color="#4a55db"><span style="color:#4A55DB">)</span></font>,<font
color="rosybrown"><span style="color:rosybrown">'<span
class="SpellE">spline'<font color="black"><span
style="color:black">,</span></font>'not</span>-a-knot'</span></font><font
color="#4a55db"><span style="color:#4A55DB">)</span></font>;<o:p></o:p></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">?<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">P.S. Since not-a-knot is
mentioned as the default for the
<span class="SpellE">splin</span> function, I think it
should also be made the default for interp1 ... just my
two cents.<o:p></o:p></span></font></p>
<p><font face="Times New Roman" color="black" size="3"><span
style="font-size:12.0pt">/Claus<o:p></o:p></span></font></p>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users@lists.scilab.org">users@lists.scilab.org</a>
<a class="moz-txt-link-freetext" href="http://lists.scilab.org/mailman/listinfo/users">http://lists.scilab.org/mailman/listinfo/users</a>
</pre>
</blockquote>
<p><br>
</p>
</body>
</html>