prettyscavenger.com

  

Beste Artikel:

  
Main / Verlet-Algorithmus in fortran 90 wobei

Verlet-Algorithmus in fortran 90 wo

In GitHub arbeiten über 36 Millionen Entwickler zusammen, um Code zu hosten und zu überprüfen, Projekte zu verwalten und gemeinsam Software zu erstellen. Zum Inhalt springen. Permalink Dismiss Werden Sie heute Mitglied bei GitHub In GitHub arbeiten über 36 Millionen Entwickler zusammen, um gemeinsam Code zu hosten und zu überprüfen, Projekte zu verwalten und Software zu erstellen.

Anmelden. Datei suchen Kopierpfad. Derzeit können keine Mitwirkenden abgerufen werden. Rohe Schuld Geschichte. Lennard-Jones-Potential, 'Geschwindigkeits'-Verlet-Zeitintegrationsalgorithmus.

Berechnet kinetische und potentielle Energie, Dichte, Druck. Von diesem Programm verwendete Dateien: Ausgabe, die geeignet ist, direkt an 'gnuplot' übergeben zu werden, um Plots als! Enthält alle Datenstrukturen mit Informationen!

Enthält die Parameter, die die Simulation steuern. Name der Datei, die die Ganzzahl des Ausgabebeispiels enthalten soll :: Anzahl der Zeitschritte, die mit doppelter Genauigkeit ausgeführt werden sollen :: Zeitschritte mit reduzierter Einheit doppelte Genauigkeit :: Gewünschte Dichte. Enthält Parameter des Lennard-Jones-Potentials! Enthält statistische Größen, die sich während des Laufs angesammelt haben. Alle Mengen werden auf Null zurückgesetzt. Dies ist der Haupttreiber für das Molekulardynamikprogramm.

Initialisierungsvorgang zu Beginn einmal aufgerufen, vor dem! Lesen Sie die Benutzeranweisungen zur Steuerung der Simulation! Lesen Sie die Eingabeprobe mit den anfänglichen Partikelkoordinaten! Drucken Sie Informationen zum Lauf in die Standardausgabedatei!

Liest das erste Beispiel aus Dateieinheit 1! N ist ', N halt endif! Die Energien werden anfänglich auf Null gesetzt. Für "neue" Samples, dh Samples, die nur durch Definieren des Atoms erstellt wurden! Für Muster, die durch frühere Simulationen hergestellt wurden, haben wir auch! Die logische Variable 'VelAcc' unterscheidet zwischen diesen beiden Fällen. Lesen Sie die Parameter, die die Simulation steuern, von der Standardeingabe. Lesen Sie die Eingabeparameter, die die Simulation steuern, aus dem Standard! Druckt Informationen zu den Laufparametern auf der Standardausgabe!

Führende sind, die Ausgabe direkt als Gnuplot-Eingabe zu verwenden. Jetzt Kopfzeilen drucken! Dies ist die Hauptunterroutine, die die zeitliche Entwicklung des Systems steuert. Wir müssen die Anfangstemperatur bereithalten, falls wir gehen! Allen und Tildesley Buch, p.

Einfache Geschwindigkeitsskalierung für Geschwindigkeiten im vorherigen Schritt! Informationen zur Druckberechnung finden Sie im Allen and Tildesley-Buch, Abschnitt 2. Partikel, die die Box verlassen haben, werden von! Berechnen Sie mit dem Lennard-Jones die Kräfte auf Atome aus den Positionen! Sij, Rij doppelte Präzision :: Rsqij, phi, dphi doppelte Präzision :: Auf null potentielle Energien, Kräfte, Virialterm zurücksetzen!

Schleife über alle Partikelpaare! Ausgehend von den aktuell in vel gespeicherten Geschwindigkeiten wird es aktualisiert! Kündigungsverfahren einmal am Ende aufgerufen, nach zeitlicher Entwicklung! Drucken Sie eine Zeile mit Durchschnittswerten usw.! Schreiben Sie die endgültigen Probenpositionen, Geschwindigkeiten und Beschleunigungen in die Datei! Geben Sie alle dynamischen Arrays frei, um die "ökologische" Speicherpraxis zu bereinigen!

Drucken Sie den während des Laufs gemittelten Mittelwert der Statistik! Schreibt das endgültige Sample in Dateieinheit 2! Multiplizieren Sie die Koordinaten, die mit der Boxgröße skaliert sind, mit BoxSize in!

Machen Sie dasselbe für Geschwindigkeiten und Beschleunigungen. Sie haben sich mit einem anderen Tab oder Fenster angemeldet. Neu laden, um Ihre Sitzung zu aktualisieren. Sie haben sich in einem anderen Tab oder Fenster abgemeldet. Ausgabe geeignet, um direkt an 'gnuplot' weitergeleitet zu werden, um Diagramme als zu erzeugen.

Enthält alle Datenstrukturen mit Informationen. Enthält die Parameter, die die Simulation steuern. Beliebige 'Titel'-Zeichenfolge. Name der Datei mit dem Eingabebeispiel. Name der Datei, die das Ausgabebeispiel enthalten soll. Anzahl der zu erledigenden Zeitschritte. Zeitschritte reduzierten Einheiten. Enthält Parameter des Lennard-Jones-Potentials. Die Initialisierungsprozedur wird zu Beginn einmal vor dem aufgerufen. Lesen Sie die Benutzeranweisungen, die die Simulation steuern.

Lesen Sie die Eingabeprobe ab, die die anfänglichen Partikelkoordinaten enthält. Drucken Sie Informationen zum Lauf in die Standardausgabedatei. Liest das erste Sample aus Dateieinheit 1. Für "neue" Samples, dh Samples, die gerade durch Definieren des Atoms erstellt wurden. Für Proben, die durch frühere Simulationen hergestellt wurden, haben wir auch.

Lesen Sie die Eingabeparameter, die die Simulation steuern, aus dem Standard. Druckt Informationen zu den Laufparametern in der Standardausgabe. Drucken Sie nun die Überschriften der Spalten. Wir müssen die Anfangstemperatur bereithalten, falls wir gehen. Einfache Geschwindigkeitsskalierung für Geschwindigkeiten im vorherigen Schritt. Partikel, die die Box verlassen haben, werden von in die Box zurückgefaltet. Berechnen Sie mit dem Lennard-Jones die Kräfte auf Atome aus den Positionen. Auf null potentielle Energien, Kräfte, Virialterm zurücksetzen.

Schleife über alle Partikelpaare. Ausgehend von den aktuell in vel gespeicherten Geschwindigkeiten wird es aktualisiert. Das Kündigungsverfahren wird am Ende nach der zeitlichen Entwicklung einmal aufgerufen. Drucken Sie eine Linie mit Durchschnittswerten usw. Schreiben Sie die endgültigen Probenpositionen, Geschwindigkeiten und Beschleunigungen in die Datei.

Geben Sie alle dynamischen Arrays frei, um die "ökologische" Speicherpraxis zu bereinigen. Drucken Sie den während des Laufs gemittelten Mittelwert der Statistik. Schreibt das endgültige Sample in Dateieinheit 2. Multiplizieren Sie die Koordinaten, die mit der Boxgröße skaliert sind, mit BoxSize in.

(с) 2019 prettyscavenger.com