PDA

View Full Version : کمک برای بازکردن کدی با فایل sif.gz.*



vahid javani
شنبه 04 آذر 1391, 18:25 عصر
سلام
چندین فایل به این شکل دارم که باید به کد c تبدیلش کنم ولی نمی دونم چطوری باز میشه و می تونم به کدهاش دسترسی داشته باشم.
تا حالا به اینجا رسیدم!
http://filext.com/file-extension/SIF

***************************
* SET UP THE INITIAL DATA *
***************************

NAME TORSIONF

* Problem :
* *********

* The quadratic elastic torsion problem

* The problem comes from the obstacle problem on a square.

* The square is discretized into (px-1)(py-1) little squares. The
* heights of the considered surface above the corners of these little
* squares are the problem variables, There are px**2 of them.

* The dimension of the problem is specified by Q, which is half the
* number discretization points along one of the coordinate
* direction.
* Since the number of variables is P**2, it is given by 4Q**2

* This is a variant of the problem stated in the report quoted below.
* It corresponds to the problem as distributed in MINPACK-2.

* Source: problem (c=20, starting point Z = origin) in
* J. More' and G. Toraldo,
* "On the Solution of Large Quadratic-Programming Problems with Bound
* Constraints",
* SIAM J. on Optimization, vol 1(1), pp. 93-113, 1991.

* SIF input: Ph. Toint, Dec 1989.
* modified by Peihuang Chen, according to MINPACK-2, Apr 1992.

* classification QBR2-MY-V-0

* Q is half the number of discretized points along the X axis

*IE Q 2 $-PARAMETER n= 16 original value
*IE Q 5 $-PARAMETER n= 100
*IE Q 11 $-PARAMETER n= 484
*IE Q 16 $-PARAMETER n= 1024
IE Q 37 $-PARAMETER n= 5476
*IE Q 50 $-PARAMETER n= 10000
*IE Q 61 $-PARAMETER n= 14884

* The force constant

RE C 20.0 $-PARAMETER force connstant

* Define a few helpful parameters

IA Q+1 Q 1
I+ P Q Q
IA P-1 P -1
RI 1/H P-1
RD H 1/H 1.0
R* H2 H H
R* C0 H2 C
RM LC C0 -1.0

IE 1 1
IE 2 2

VARIABLES

* Define one variable per discretized point in the unit square

DO J 1 P
DO I 1 P
X X(I,J)
ND

GROUPS


DO I 2 P
DO J 2 P
ZN GL(I,J)
ND

DO I 1 P-1
DO J 1 P-1
ZN GR(I,J)
ND

DO I 2 P-1
DO J 2 P-1
ZN G X(I,J) LC
ND


BOUNDS

* Fix the variables on the lower and upper edges of the unit square

DO J 1 P
XX TORSIONF X(1,J) 0.0
XX TORSIONF X(P,J) 0.0
ND

* Fix the variables on the left and right edges of the unit square

DO I 2 P-1
XX TORSIONF X(I,P) 0.0
XX TORSIONF X(I,1) 0.0
ND

* Define the upper and lower bounds from the distance to the
* boundary

* Lower half of the square

DO I 2 Q

DO J 2 I
IA J-1 J -1
RI RJ-1 J-1
R* UPPL RJ-1 H
RM LOWL UPPL -1.0
ZL TORSIONF X(I,J) LOWL
ZU TORSIONF X(I,J) UPPL
OD J

IM MI I -1
I+ P-I P MI
IA I-1 I -1
RI RI-1 I-1
R* UPPM RI-1 H
RM LOWM UPPM -1.0
IA P-I+1 P-I 1

DO J I P-I+1
ZL TORSIONF X(I,J) LOWM
ZU TORSIONF X(I,J) UPPM
OD J

DO J P-I+1 P-1
IM MJ J -1
I+ P-J P MJ
RI RP-J P-J
R* UPPR RP-J H
RM LOWR UPPR -1.0
ZL TORSIONF X(I,J) LOWR
ZU TORSIONF X(I,J) UPPR

ND

* Upper half of the square

DO I Q+1 P-1

IM MI I -1
I+ P-I P MI
IA P-I+1 P-I 1

DO J 2 P-I+1

IA J-1 J -1
RI RJ-1 J-1
R* UPPL RJ-1 H
RM LOWL UPPL -1.0

ZL TORSIONF X(I,J) LOWL
ZU TORSIONF X(I,J) UPPL

OD J

RI RP-I P-I
R* UPPM RP-I H
RM LOWM UPPM -1.0

DO J P-I+1 I

ZL TORSIONF X(I,J) LOWM
ZU TORSIONF X(I,J) UPPM

OD J

DO J I P-1

IM MJ J -1
I+ P-J P MJ
RI RP-J P-J
R* UPPR RP-J H
RM LOWR UPPR -1.0

ZL TORSIONF X(I,J) LOWR
ZU TORSIONF X(I,J) UPPR

ND

ELEMENT TYPE

EV ISQ V1 V2
IV ISQ U

ELEMENT USES



DO I 1 P-1

IA I+1 I 1

DO J 1 P-1

IA J+1 J 1

XT A(I,J) ISQ
ZV A(I,J) V1 X(I+1,J)
ZV A(I,J) V2 X(I,J)

XT B(I,J) ISQ
ZV B(I,J) V1 X(I,J+1)
ZV B(I,J) V2 X(I,J)

ND

DO I 2 P

IA I-1 I -1

DO J 2 P

IA J-1 J -1

XT C(I,J) ISQ
ZV C(I,J) V1 X(I-1,J)
ZV C(I,J) V2 X(I,J)

XT D(I,J) ISQ
ZV D(I,J) V1 X(I,J-1)
ZV D(I,J) V2 X(I,J)

ND

GROUP USES

DO I 1 P-1
DO J 1 P-1
XE GR(I,J) A(I,J) 0.25 B(I,J) 0.25
ND
*
DO I 2 P
DO J 2 P
XE GL(I,J) C(I,J) 0.25 D(I,J) 0.25
ND


OBJECT BOUND

* Solution

*LO SOLTN(2) -2.518500
*LO SOLTN(5) -2.798400
*LO SOLTN(11) -2.840600
*LO SOLTN(16) -2.846100
*LO SOLTN(37) -2.850200
*LO SOLTN(50) -2.850700
*LO SOLTN(61) -2.850800

ENDATA

***********************
* SET UP THE FUNCTION *
* AND RANGE ROUTINES *
***********************

ELEMENTS TORSIONF

INDIVIDUALS

T ISQ
R U V1 1.0 V2 -1.0
F U * U
G U U + U
H U U 2.0

ENDATA


این می تونه یه کد باشه؟ میشه راهنماییم کنید چطوری باید به c تبدیلش کنم؟؟
پیشاپیش ممنونم از پاسختون...!

shahmohammadi
شنبه 11 آذر 1391, 17:44 عصر
سلام.
از توضيحاتتون چيزي متوجه نشدم. ولي شما يكي از فايل ها رو اينجا بگذاريد تا كمكتون بكينم.

vahid javani
یک شنبه 12 آذر 1391, 12:57 عصر
ممنون از پاسختون دوست عزیز
http://www.rodfile.com/ic87dxisk44u/Javani.rar.html