Skip to content

Commit

Permalink
apply ttws
Browse files Browse the repository at this point in the history
  • Loading branch information
thorade committed Dec 9, 2014
1 parent 1afaa4f commit e9bfbc7
Show file tree
Hide file tree
Showing 3 changed files with 18 additions and 18 deletions.
20 changes: 10 additions & 10 deletions .gitattributes
Original file line number Diff line number Diff line change
@@ -1,15 +1,15 @@
# Explicitly declare all file extensions to be text or binary
*.csv text
*.m text
*.md text
*.csv text
*.m text
*.md text
*.mo text
*.mos text
*.order text
*.pdf binary
*.rfp text
*.svg text
*.tex text
*.txt text
*.mos text
*.order text
*.pdf binary
*.rfp text
*.svg text
*.tex text
*.txt text

# Explicitly declare languages, see https://github.com/github/linguist#overrides
*.m linguist-language=Matlab
Expand Down
6 changes: 3 additions & 3 deletions Extras/LaTeX/tables/RewriteDerivatives.m
Original file line number Diff line number Diff line change
Expand Up @@ -18,9 +18,9 @@
dP_D = T*R*(1 + 2*frd + frdd);
dP_T = D*R*(1 + frd - frtd);
dU_T = R*(-fitt - frtt);
%
%
d2P_D2 = T*R/D*(2*frd + 4*frdd + frddd);
d2P_T2 = D*R/T*(frttd);
d2P_T2 = D*R/T*(frttd);
d2P_TD = R *(1 + 2*frd + frdd - 2*frtd - frtdd);
d2U_T2 = R/T*(fittt + frttt + 2*fitt + 2*frtt);

Expand Down Expand Up @@ -89,7 +89,7 @@

%% PVT second order derivatives
d2T_D2_P = -(d2P_D2*dP_T - dP_D*d2P_TD)/dP_T^2 ...
+(d2P_TD*dP_T - dP_D*d2P_T2)*dP_D/dP_T^3
+(d2P_TD*dP_T - dP_D*d2P_T2)*dP_D/dP_T^3
d2T_D2_P = simplify(d2T_D2_P)
d2T_D2_P = expand(d2T_D2_P)

Expand Down
10 changes: 5 additions & 5 deletions Extras/ReferenceState/ReferenceState.m
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
% manually delete ReferenceState.csv
% uncomment print header in ReferenceState.mo
% set idealPower[1,1] and idealPower[2,1] to 0
% set idealPower[1,1] and idealPower[2,1] to 0
% run ReferenceState.mo
% comment print header
% change idealPower[1,1] and [2,1] to some value different from 0
Expand All @@ -15,20 +15,20 @@
% ------------------------------------------------------------------------
% y=m*x+b
% for y=0 x=-b/m
% for y=200
% for y=200
% ------------------------------------------------------------------------

% Entropy
% find linear fit for reverse function
% find linear fit for reverse function
[polyS,StructureS,muS] = polyfit(sref,idealPower1,1)
% evaluate reverse function at desired value
sref0=1e3;
[idealPower1_ref,deltaS] = polyval(polyS,sref0,StructureS,muS)


% Enthalpy
% find linear fit for reverse function
% find linear fit for reverse function
[polyH,StructureH,muH] = polyfit(href,idealPower2,1)
% evaluate reverse function at desired href value
href0=2e5;
[idealPower2_ref,deltaH] = polyval(polyH,href0,StructureH,muH)
[idealPower2_ref,deltaH] = polyval(polyH,href0,StructureH,muH)

0 comments on commit e9bfbc7

Please sign in to comment.