Re: dependency diagram for druckprag

[email protected] Sun, 26 Feb 2012 14:11:23 +0800
Newsgroups gmane.comp.mathematics.tochnog.user
Organization HKUST
Message-ID <[email protected]>
This is a multi-part message in MIME format.
--===============7027815528544014108==
Content-Type: multipart/alternative;
	boundary="------------070804000506060400090208"

This is a multi-part message in MIME format.
--------------070804000506060400090208
Content-Type: text/plain; charset=windows-1252; format=flowed
Content-Transfer-Encoding: quoted-printable
X-MIME-Autoconverted: from 8bit to quoted-printable by mx3.ust.hk id
	q1Q6ADcw036335

On Sunday, February 26, 2012 10:35 AM,=20
[email protected] wrote:

Hi Fernando,

Thank you very much for the detailed explanations and kind suggestions.=20
As you pointed out, I wrongly used the dependency_diagram for multiple=20
parameters. With your corrected script, the program works smoothly and fa=
st.

For options_solver, I cannot test the Lapack solver in my Ubuntu desktop=20
(`I do not know : -matrix_lapack`). How to install extra solvers and=20
make them accessible by tochnog?

Best regards,
Ning

> Ning,
>
> I read your comment and still believe that you are not entering the dep=
endency diagram items the proper way.
> I checked the code and the open source version allows to enter more tha=
n one property. See the attached note
> from the manual on dependency_item:
>
> "=85.
> Temperature dependent Young's moduli in two maxwell chains of element_g=
roup 1 ( for the first chain the
> Moduli 1.e10, 1.e9,...  for the second chain the moduli 1.e12, 1.e11,..=
.  all relaxation times are 1.e-2 . ):
>
> dependency_item 1 -group_materi_maxwell_chain 1 -temp 4
> dependency_diagram 1
> 1. 2. 3. 4.
> 1.e10 1.e9 1.e8 3.e5
> 1.e-2 1.e-2 1.e-2 1.e-2
> 1.e12 1.e11 1.e10 3.e7
> 1.e-2 1.e-2 1.e-2 1.e-2
> =85"
>
> As you can see, you enter the list of unknowns, in the example above, t=
he values of -temp, then you enter the list
> of values of the
> first chain moduli, and then the list of the  relaxation time values, t=
hen the second list of the second chain moduli,
> and then the list of the second relaxation times for the secon chain.
>
> for your case, upon reviewing your file, you have
>
> "=85
> group_materi_plasti_druckprag 1 0.611 1.5e6 0.175
> group_integration_points 1 -maximal
>
> dependency_item 1 -group_materi_plasti_druckprag 1 -kap 2
> dependency_diagram 1 0. 0.1 0.611 1.5e6 0.175 0.611 5246.89 0.175
> =85"
>
> Therefore, I understand that  you are using:
> group_materi_plasti_druckprag index  phi c phiflow
>
> and that your inital your initial values are phi=3D 0.611, c =3D 1.5e6 =
and phiflow=3D 0.175
>
> You are keeping phi and phiflow constant and are modifying c. Your depe=
ndency diagram should therefore be as follows:
>
> ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
>
> dependency_item 1 -group_materi_plasti_druckprag 1 -kap 2
> dependency_diagram 1 0. 0.1  0.611 0.611 1.5e6 5246.89 0.175 0.175
>
> ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
>
> So I am giving two values of -kap, two values of phi, two values of c, =
and two values of phiflow.
>
> I changed the values and it runs well and very fast. I am attaching a c=
opy of the modified file for your review.
>
> Just a couple of reminders:
>
> 1.	The updated, that is current version of tochnog allows the use of th=
e Lapack solver which is faster than the
> 	default biconjugate solver. I built the version for Mac OSX but curren=
tly do not have access to Linux or Windows compilers.
>
> 2.	Rather than using "group_integration_points 1 -maximal" you can use =
the modern version of tochnog: "number_of_integration_points N" ,
> 	where N is any number of integration points that you want to use. If y=
ou read the manual you will learn that  -maximal has a limit,
> 	 but "number_of_integration_points" does not. Just be reasonable and u=
nderstand how many points you will need.
> 3.	If you use "control_print_gid "number"  -yes"  and use the same numb=
er or index for n as for control_timestep, then you can print
> 	and plot (at least in GiD) the results for any time step and not only =
at the end of the calculations.
> 4.	Finally, I keep suggesting to use the automatic iterations option in=
 tochnog, as it allows you to control the reslution or error  of the
> 	problem. I included that in the modified file attached.
>
> Good luck and please let me know if this works.
>
> Best regards
>
>
>
>
> Fernando Lorenzo
> [email protected]
>
>
>
>
>
> -----------------------------------------------------------------------=
-------
> Virtualization&  Cloud Management Using Capacity Planning
> Cloud computing makes use of virtualization - but cloud computing
> also focuses on allowing computing to be delivered as a service.
> http://www.accelacomm.com/jaw/sfnl/114/51521223/
>
>
> _______________________________________________
> Tochnog-users mailing list
> [email protected]
> https://lists.sourceforge.net/lists/listinfo/tochnog-users


--=20
Geotechnical Group
Department of Civil and Environmental Engineering
Hong Kong University of Science and Technology
Clear Water Bay, Kowloon, Hong Kong


--------------070804000506060400090208
Content-Type: text/html; charset=windows-1252
Content-Transfer-Encoding: quoted-printable
X-MIME-Autoconverted: from 8bit to quoted-printable by mx3.ust.hk id
	q1Q6ADcw036335

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content=3D"text/html; charset=3Dwindows-1252"
      http-equiv=3D"Content-Type">
    <title></title>
  </head>
  <body text=3D"#000000" bgcolor=3D"#ffffff">
    On Sunday, February 26, 2012 10:35 AM,
    <a class=3D"moz-txt-link-abbreviated" href=3D"mailto:tochnog-users@li=
sts.sourceforge.net">[email protected]</a> wrote:<br>
    <br>
    Hi Fernando,<br>
    <br>
    Thank you very much for the detailed explanations and kind
    suggestions. As you pointed out, I wrongly used the
    dependency_diagram for multiple parameters. With your corrected
    script, the program works smoothly and fast.<br>
    <br>
    For options_solver, I cannot test the Lapack solver in my Ubuntu
    desktop (`I do not know : -matrix_lapack`). How to install extra
    solvers and make them accessible by tochnog?<br>
    <br>
    Best regards,<br>
    Ning<br>
    <br>
    <blockquote
      cite=3D"mid:[email protected]"
      type=3D"cite">
      <pre wrap=3D"">Ning,

I read your comment and still believe that you are not entering the depen=
dency diagram items the proper way.=20
I checked the code and the open source version allows to enter more than =
one property. See the attached note=20
from the manual on dependency_item:

"=85.
Temperature dependent Young's moduli in two maxwell chains of element_gro=
up 1 ( for the first chain the=20
Moduli 1.e10, 1.e9,...  for the second chain the moduli 1.e12, 1.e11,... =
 all relaxation times are 1.e-2 . ):

dependency_item 1 -group_materi_maxwell_chain 1 -temp 4=20
dependency_diagram 1=20
1. 2. 3. 4.=20
1.e10 1.e9 1.e8 3.e5=20
1.e-2 1.e-2 1.e-2 1.e-2=20
1.e12 1.e11 1.e10 3.e7=20
1.e-2 1.e-2 1.e-2 1.e-2=20
=85"

As you can see, you enter the list of unknowns, in the example above, the=
 values of -temp, then you enter the list=20
of values of the
first chain moduli, and then the list of the  relaxation time values, the=
n the second list of the second chain moduli,
and then the list of the second relaxation times for the secon chain.

for your case, upon reviewing your file, you have=20

"=85
group_materi_plasti_druckprag 1 0.611 1.5e6 0.175
group_integration_points 1 -maximal

dependency_item 1 -group_materi_plasti_druckprag 1 -kap 2
dependency_diagram 1 0. 0.1 0.611 1.5e6 0.175 0.611 5246.89 0.175
=85"

Therefore, I understand that  you are using:
group_materi_plasti_druckprag index  phi c phiflow

and that your inital your initial values are phi=3D 0.611, c =3D 1.5e6 an=
d phiflow=3D 0.175

You are keeping phi and phiflow constant and are modifying c. Your depend=
ency diagram should therefore be as follows:

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

dependency_item 1 -group_materi_plasti_druckprag 1 -kap 2
dependency_diagram 1 0. 0.1  0.611 0.611 1.5e6 5246.89 0.175 0.175

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

So I am giving two values of -kap, two values of phi, two values of c, an=
d two values of phiflow.

I changed the values and it runs well and very fast. I am attaching a cop=
y of the modified file for your review.

Just a couple of reminders:=20

1.	The updated, that is current version of tochnog allows the use of the =
Lapack solver which is faster than the=20
	default biconjugate solver. I built the version for Mac OSX but currentl=
y do not have access to Linux or Windows compilers.

2.	Rather than using "group_integration_points 1 -maximal" you can use th=
e modern version of tochnog: "number_of_integration_points N" ,
	where N is any number of integration points that you want to use. If you=
 read the manual you will learn that  -maximal has a limit,
	 but "number_of_integration_points" does not. Just be reasonable and und=
erstand how many points you will need.
3.	If you use "control_print_gid "number"  -yes"  and use the same number=
 or index for n as for control_timestep, then you can print=20
	and plot (at least in GiD) the results for any time step and not only at=
 the end of the calculations.
4.	Finally, I keep suggesting to use the automatic iterations option in t=
ochnog, as it allows you to control the reslution or error  of the=20
	problem. I included that in the modified file attached.

Good luck and please let me know if this works.

Best regards</pre>
      <pre wrap=3D"">
<fieldset class=3D"mimeAttachmentHeader"></fieldset>


Fernando Lorenzo
<a class=3D"moz-txt-link-abbreviated" href=3D"mailto:[email protected]=
">[email protected]</a>



</pre>
      <pre wrap=3D"">
<fieldset class=3D"mimeAttachmentHeader"></fieldset>
-------------------------------------------------------------------------=
-----
Virtualization &amp; Cloud Management Using Capacity Planning
Cloud computing makes use of virtualization - but cloud computing=20
also focuses on allowing computing to be delivered as a service.
<a class=3D"moz-txt-link-freetext" href=3D"http://www.accelacomm.com/jaw/=
sfnl/114/51521223/">http://www.accelacomm.com/jaw/sfnl/114/51521223/</a><=
/pre>
      <pre wrap=3D"">
<fieldset class=3D"mimeAttachmentHeader"></fieldset>
_______________________________________________
Tochnog-users mailing list
<a class=3D"moz-txt-link-abbreviated" href=3D"mailto:Tochnog-users@lists.=
sourceforge.net">[email protected]</a>
<a class=3D"moz-txt-link-freetext" href=3D"https://lists.sourceforge.net/=
lists/listinfo/tochnog-users">https://lists.sourceforge.net/lists/listinf=
o/tochnog-users</a>
</pre>
    </blockquote>
    <br>
    <br>
    <pre class=3D"moz-signature" cols=3D"72">--=20
Geotechnical Group
Department of Civil and Environmental Engineering
Hong Kong University of Science and Technology
Clear Water Bay, Kowloon, Hong Kong</pre>
  </body>
</html>

--------------070804000506060400090208--


--===============7027815528544014108==
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

------------------------------------------------------------------------------
Virtualization & Cloud Management Using Capacity Planning
Cloud computing makes use of virtualization - but cloud computing 
also focuses on allowing computing to be delivered as a service.
http://www.accelacomm.com/jaw/sfnl/114/51521223/
--===============7027815528544014108==
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

_______________________________________________
Tochnog-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/tochnog-users

--===============7027815528544014108==--