Rechenprogramm

[Michael II]

Hallo Profi

weisst du wie man deine hlp in Win10 öffnet? Wenn ich das tun will, dann wird eine m$ Seite angezeigt mit Downloads für Systeme bis und mit Win 8.1. Ich habe dann mal 8.1 geladen - aber das führt einfach nur zu einer weiteren Fehlermeldung.

Wichtigere Frage... Wieso kann ich mit deinem Rechner "nur" bis 1754! rechnen?

Eine "tief" Funktion ( n tief k = n!/(k!*(n-k)!) ) wäre auch gut. Evt. hast du's ja bereits eingebaut und ich seh's einfach nicht.

Gruss
M

[Amateurprofi]

Zitat von Michael II:

Hallo Profi

weisst du wie man deine hlp in Win10 öffnet? Wenn ich das tun will, dann wird eine m$ Seite angezeigt mit Downloads für Systeme bis und mit Win 8.1. Ich habe dann mal 8.1 geladen - aber das führt einfach nur zu einer weiteren Fehlermeldung.

Wichtigere Frage... Wieso kann ich mit deinem Rechner "nur" bis 1754! rechnen?

Eine "tief" Funktion ( n tief k = n!/(k!*(n-k)!) ) wäre auch gut. Evt. hast du's ja bereits eingebaut und ich seh's einfach nicht.

Gruss
M

Zitat:

Eine "tief" Funktion ( n tief k = n!/(k!*(n-k)!) ) wäre auch gut.

Gib ein Varis(50,3) oder Varis(50,3,0)
Variationen ohne Wiederholung
Varis(50,3,1)
Variationen mit Wiederholung
Combis = Kombinationen mit oder ohne Wiederholung. Gleiche Syntax wie Varis.

Zitat:

weisst du wie man deine hlp in Win10 öffnet?

Nein, weiss ich leider nicht, ich arbeite mit Windows 7.

[Amateurprofi]

Zitat von Michael II:

Hallo Profi
Wichtigere Frage... Wieso kann ich mit deinem Rechner "nur" bis 1754! rechnen?
Gruss
M

Hab ich mit PN beantwortet, hier noch mal für andere.

Das Maximum 1754 für N! (im Normal-Modus) resultiert aus dem größten in einem Extended darstellbaren Wert (1.1E+4932).
1755! ergäbe 3.473E+4933 und würde den Maximalwert überschreiten.

Übrigens:
Im Normal-Modus werden in meinem Programm alle Werte in Extended-Variablen gespeichert, was auch der Grund ist, warum es keine 64Bit Version gibt (Da ist Extended = Double).

Tatsächlich rechne ich N! mit der FPU. Hier der Code

Delphi-Quellcode:

PROCEDURE MonOpFacul(var V:TValue);
asm // EAX=@TValue
// CheckMinMax(v.v,0,1754,pact,116);
// n:=Trunc(v.v);
// v.v:=1;
// while n>1 do begin
// v.v:=v.v*n;
// dec(n);
// end;
          // ASM Version braucht 1/6 der Zeit und ist kürzer, auch wenns
          // nicht so aussieht.
                   mov byte [eax].TValue.Vt,vtNone
                   fld tbyte[eax].TValue.V
                   fldz
                   fcomip st,st(1)
                   ja @Err
                   fild dword [@1754]
                   fcomip st,st(1)
                   jae @Ok
 @Err: mov eax,117
                   mov edx,pact
                   jmp ASMError1
 @1754: dd 1754
 @Ok: sub esp,8
                   fnstcw word [esp+4]
                   fnstcw word [esp+6]
                   fwait
                   or word [ESP+4], $0F00 // trunc toward zero, full precision
                   fldcw word [ESP+4]
                   fistp dword [esp]
                   fwait
                   fldcw word [esp+6]
                   fld1
                   mov edx,[esp]
                   sub edx,1
                   jbe @Fin
 @Loop: fimul dword [esp]
                   mov [esp],edx
                   sub edx,1
                   ja @Loop
 @Fin: fstp tbyte[eax].TValue.V
                   add esp,8
end;

Delphi-Quellcode:

   // Vom Parser benutzter Werte-Record
    // CanDeformat : Info, ob der von NumStr erzeugte Ausgabe-String
    // wieder deformatiert werden kann
    // Vt : Wofür der Wert steht
    // V : Der numerische Wert
    // Numerator : Zähler, wenn Vt=vtFraction
    // Denominator : Nenner, wenn Vt=vtFraction
    // Latitude : Breitengrad, wenn Vt=vtGeoPos
    // Longitude : Längengrad, wenn Vt=vtGeoPos
    // ClrVal : Fabbdaten, wenn Vt=vtColor
    TValue=packed Record
       CanDeformat:Boolean;
       Vt:TValueType;
       V:Extended;
       case Integer of
          1 : (Numerator:Extended;
                 Denominator:Extended);
          2 : (Latitude:Extended;
                 Longitude:Extended);
          3 : (LatLon:Array[Boolean] of Extended);
          4 : (ClrVal:TColorValue);
    end; 
  // Wofür ein Wert steht (vtNone=einfache Zahl)
    TValueType=(vtNone,vtBool,vtColor,vtFraction,
                vtSeconds,vtMinutes,vtHours,vtDays,vtWeeks,vtMonths,vtYears,
                vtTime,vtDate,vtDateTime,
                vtGeoPos);

[Michael II]

Besten Dank!

Du könntest
x = e^ln(x) bzw. hier besser
x = 10^lg(x) und
log(ab)=log(a)+log(b)
ausnutzen.

Dann kannst du x! schreiben als x! = 10^lg(1*2*3...*x) = 10^(lg(1)+...lg(x)).

Fürs Coden: Du berechnest zuerst
s = lg(1)+...lg(x)

Sei s = e + m, wobei e = trunc(s) und m = frac(s)

x! = 10^s = 10^(e+m) = 10^m * 10^e.

Output:
10^m rechnest du aus, e lässt du als Exponent stehen.

Ich fühl mich gerade zurückversetzt in die Zeiten des hp41CV... .

[Amateurprofi]

Zitat von Michael II:

Besten Dank!

Du könntest
x = e^ln(x) bzw. hier besser
x = 10^lg(x) und
log(ab)=log(a)+log(b)
ausnutzen.

Dann kannst du x! schreiben als x! = 10^lg(1*2*3...*x) = 10^(lg(1)+...lg(x)).

Fürs Coden: Du berechnest zuerst
s = lg(1)+...lg(x)

Sei s = e + m, wobei e = trunc(s) und m = frac(s)

x! = 10^s = 10^(e+m) = 10^m * 10^e.

Output:
10^m rechnest du aus, e lässt du als Exponent stehen.

Ich fühl mich gerade zurückversetzt in die Zeiten des hp41CV... .

Das wäre aber eine schöne Verschlimmbesserung.
Statt (für N=1000) hätte ich statt 1000 Multiplikationen 1000 mal Log10(), und Log10() braucht seine Zeit.

[Michael II]

Du hast natürlich Recht. Das von mir vorgeschlagene "Verfahren" ist mathematischer Bocks-Mist. Es führt zwar aufs richtige Resultat. Aber der Weg.... [Ich war soeben beim Ponys Misten und obschon es auf 800müM bereits ziemlich kalt und windig ist, war die Jacke plötzlich völlig überflüssig als mir einfiel, dass das was ich geschrieben.... Du kannst natürlich direkt Mantisse und Exponent separieren und über die 1754! Grenze hinaus multiplizieren... ich lasse es zur Belustigung stehen... ]


Ich schlage dieses Verfahren ja nicht für "kleine" Werte von x (also nicht für x <= 1754) vor. Mir ist vollkommen klar, dass der Aufwand via log10 weit grösser ist. Wenn ich dich aber richtig verstehe, kann dein Programm Näherungswerte "nur" bis 1754! berechnen, dann ist aus.

Wenn du aber über den Exponenten zur Basis 10 rechnest, kommst du an Orte, wo nie jemand zuvor gewesen ist.

Du könntest deine Superschnellmethode bis 1754 nutzen und ab dann via log die Näherungswerte berechnen.

Für grosse Werte - kommen in der Statistik durchaus vor - könnte dein Programm zum Beispiel jeden 100sten ! Wert tabellieren; dann müsstest du nicht allzu viele log10 addieren.

War nur ein Vorschlag - mehr nicht.

Besten Dank für deine Tipps zu "tief" und Co (Kombinatorik Teil in deinem Calculator).

[Amateurprofi]

Zitat von Michael II:

Ich schlage dieses Verfahren ja nicht für "kleine" Werte von x (also nicht für x <= 1754) vor. Mir ist vollkommen klar, dass der Aufwand via log10 weit grösser ist. Wenn ich dich aber richtig verstehe, kann dein Programm Näherungswerte "nur" bis 1754! berechnen, dann ist aus.

In dem Programm werden im "Normal"-Modus zur Speicherung der Daten Extended-Variablen verwendet, und da gibt es nun einmal die Obergrenze 1.1E+4932

Anders ist es im "BigDecimal"-Modus.
Da kannst du die Fakultät bis 10000 berechnen. Exact übrigens und nicht nur näherungsweise.
Die 10000 als Obergrenze für Fakultät ist willkürlich festgelegt.

Allerdings können diese im Ausgabefeld nicht in voller Länge angezeigt werden.
Das hatte ich hier https://www.delphipraxis.net/195060-...extlaenge.html einmal angesprochen.

Aber mache einmal folgendes:
Stelle auf "BigDecimal"-Modus (2tes Panel in der Statusbar)
Stelle "Auch Integerzahlen in technischer Darstellung" auf "Nein" (9tes Panel in der Statusbar)
Gib 10000! ein.
Als Rechenergebnis erhältst du "Zahl zu lang" in grüner Schrift.
Rechtsklicke den Text.
In einem PopUp-Menu wird angeboten, die Zahl ins Clipboard zu stellen oder in einer Datei zu speichern.
Wenn du in Datei speichern wählst, wird in der Datei "Log.txt" im Programmverzeichnis gespeichert.
Oder rechtsklicke den Text mit gedrückter Ctrl-Taste.
Dann wird der "Dialog für lange Zahlen" geöffnet, in dem du die Zahl (in Worten) ansehen kannst.
Im Anhang findest du 2 Dateien, die Zahl in Ziffern, die Zahl in Worten.

[jaenicke]

Zitat von Michael II:

weisst du wie man deine hlp in Win10 öffnet?

Die .hlp Hilfedateien wurden vor über 20 Jahren mit Windows 98 abgekündigt und durch .chm ersetzt. Mit Vista wurde die Unterstützung dann ganz entfernt und konnte bis 8.1 noch per Download ergänzt werden.
Unter Windows 10 geht es schließlich nur noch umständlich:
https://entwickler-ecke.de/viewtopic...=697956#697956

[KodeZwerg]

Zitat von jaenicke:

Zitat von Michael II:

weisst du wie man deine hlp in Win10 öffnet?

Die .hlp Hilfedateien wurden vor über 20 Jahren mit Windows 98 abgekündigt und durch .chm ersetzt. Mit Vista wurde die Unterstützung dann ganz entfernt und konnte bis 8.1 noch per Download ergänzt werden.
Unter Windows 10 geht es schließlich nur noch umständlich:
https://entwickler-ecke.de/viewtopic...=697956#697956

Da ich keine Windows XP/Vista/7 Dateien besitze...
habe ich gerade eben das hier entdeckt: winhlp32_for_windows_10
Das ist ein Installer der WinHlp32 nachrüstet, funktioniert auch wenn man einen Doppelklick auf eine .hlp Datei macht.

Blöd ist aber es funktioniert nicht innerhalb der Applikation, da wird immer noch diese Microsoft Seite geladen.
Hat dazu jemand einen Tipp?

[jaenicke]

Zitat von KodeZwerg:

Da ich keine Windows XP/Vista/7 Dateien besitze...
habe ich gerade eben das hier entdeckt: winhlp32_for_windows_10
Das ist ein Installer der WinHlp32 nachrüstet, funktioniert auch wenn man einen Doppelklick auf eine .hlp Datei macht.

Blöd ist aber es funktioniert nicht innerhalb der Applikation, da wird immer noch diese Microsoft Seite geladen.
Hat dazu jemand einen Tipp?

Na, jetzt hast du die Datei aus dem Zusatzinstaller von Windows 8.1 Update 1 ja, denn die steckt ja in dem Installer drin. (Das siehst du, wenn du den auspackst und dann das darin enthaltene .cab Archiv auspackst.)

Ich vermute einmal der Installer legt die Datei wegen des Rechteproblems irgendwo anders hin und ersetzt die Originaldatei nicht. Du kannst also diese winhlp32.exe nehmen und damit meiner Anleitung folgen. Das sollte reichen.