[room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Aleksander N. Gorohovski]

On Wed, 02 Nov 2005 13:51:19 -0000, Michael Shigorin wrote:

> On Wed, Nov 02, 2005 at 03:33:11PM -0000, Aleksander N. Gorohovski wrote:
>> >>У него свой язык, типа программирования, но, по моему, лучше
>> >>уж писать на С.
>> >Математику?  Ну-ну.
>> Сурьезную математику, а не беллетристику.
>
> Сурьёзную математику на Це писать сильно дорого выходит.
>
> Разберёте такую мелочь, как аттачик (в смысле "что оно делает")
Миша,
я до конца не понял ваш язык изопа?
What is "аттачик"?
:-)

Re: [room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Mike Lykov]

В сообщении от Четверг 03 Ноябрь 2005 13:43 Aleksander N. Gorohovski написал:

> я до конца не понял ваш язык изопа?

язык из опа (ирц оператора)?
язык из зопа (..) ? ;)

-- 
Mike Lykov
Samara, "Vesna" parfum company, System administrator 

Re: [room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Genix]

Mike Lykov wrote:

>>я до конца не понял ваш язык изопа?
> 
> 
> язык из опа (ирц оператора)?
> язык из зопа (..) ? ;)

язык и зопа!



-- 
У каждого в башке свои тараканы...

Re: [room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Michael Shigorin]

On Thu, Nov 03, 2005 at 03:46:00PM +0400, Mike Lykov wrote:
> > я до конца не понял ваш язык изопа?
> язык из опа (ирц оператора)?
> язык из зопа (..) ? ;)

Питон это. :]

-- 
 ---- WBR, Michael Shigorin <mike@altlinux.ru>
  ------ Linux.Kiev http://www.linux.kiev.ua/

[room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Michael Shigorin]

[-- Attachment #1: Type: text/plain, Size: 622 bytes --]

On Thu, Nov 03, 2005 at 09:43:30AM -0000, Aleksander N. Gorohovski wrote:
> >>>>У него свой язык, типа программирования, но, по моему, лучше
> >>>>уж писать на С.
> >>>Математику?  Ну-ну.
> >>Сурьезную математику, а не беллетристику.
> >Сурьёзную математику на Це писать сильно дорого выходит.
> >Разберёте такую мелочь, как аттачик (в смысле "что оно делает")
> Миша, я до конца не понял ваш язык изопа?

Вообще-то это был язык C. :]

> What is "аттачик"?

an attach.  Цепляю ещё раз (хотя и тогда вроде не забыл?).

> :-)

-- 
 ---- WBR, Michael Shigorin <mike@altlinux.ru>
  ------ Linux.Kiev http://www.linux.kiev.ua/

[-- Attachment #2: alkanol.c --]
[-- Type: text/plain, Size: 4243 bytes --]

/*
 * ALCANOL -- the alcanol isomeres number calculations
 * Version 0.20 of May 10 2000 by Michael Shigorin <michael.shigorin@usa.net>
 *
 * THIS IS FREE SOFTWARE - USE FREELY BUT ON YOUR OWN RISK (CISK;)
 * MAY BE REDISTRIBUTED/MODIFIED UNDER THE TERMS OF GNU GPL ][
 *
 */
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <stdio.h>
#include <string.h>
#include "gmp.h"

#define counter_t unsigned int		/* used for counters -- 32bit for me */

//#define TIMER				/* timing support */
//#define DEBUG				/* more verbosity */
#define MAX_DEGREE	4096		/* CPU/RAM/time's the limit :) */
#define MAX_LEN		255		/* max length of iofile name */
#define FORMAT		"%d:\t"		/* ^data file format excl. HUGE numbers$ */
#ifdef TIMER
  #define LOGFORMAT	"%d\t%f"	/* number <tab> seconds */
  #define TIMEXT	".log"		/* extension to add for log */
  #include <sys/time.h>
  #include <unistd.h>
#endif

mpz_t source [MAX_DEGREE], factor, buffer;
mpz_t temp;

FILE *iofile;
char filename [MAX_LEN];
counter_t max;

#ifdef TIMER
  char logname [MAX_LEN];
  FILE *logfile;
  float timer;
  struct timeval raw_timer;
  struct timezone tz;		/* stub */

  float alk_time (void);
#endif

void init (void);
void parse (void);
void calculate (void);
void print (void);
void usage (void);
void shutdown (void);

int main (int argc, char *argv[]) {
	if( argc != 2 ) usage ();
	strncpy (filename, argv [1], strlen (argv [1])+1);

#ifdef TIMER
	snprintf (logname, strlen (filename)+strlen (TIMEXT)+1, "%s%s", filename, TIMEXT);
#endif
	init ();
	parse ();
	calculate ();
	print ();
	shutdown ();
	return (0);
}

#ifdef TIMER
float alk_time (void) {
	gettimeofday (&raw_timer, &tz);
	return (float)raw_timer.tv_sec + (float)raw_timer.tv_usec/1000000;
}
#endif

void usage (void) {
	fprintf (stderr, "I need one argument (a file to proceed)!\n");
	exit (1);
}

void init (void) {
#ifdef DEBUG
	fprintf (stderr, "\nReading data from <%s>...\n", filename);
#endif

	if( !(iofile = fopen (filename, "r")) ) {
		fprintf (stderr, "Error: cannot Open Source file <%s>\n", filename);
		exit (1);
	}
}

void parse (void) {
	mpz_init (temp);

#ifdef DEBUG
	fprintf (stderr, "Parsing...");
#endif
	while( !feof(iofile) ) {
		if( fscanf(iofile, FORMAT, &max) ) {
			if( mpz_inp_str (temp, iofile, 10) ) mpz_set (source [max], temp);
		} else {
			fprintf (stderr, "Error: malformed line #%d", max);
			exit (2);
		}
	}
#ifdef DEBUG
	fprintf (stderr, " ok, maximal degree is %d.\n", max);
#endif
}

void calculate (void) {
	counter_t i, j;

#ifdef TIMER
	timer = alk_time ();
#endif
#ifdef DEBUG
	fprintf (stderr, "Calculating factor by degree %d.\n", max+1);
#endif

	mpz_init_set_ui (factor, 0);
	/* first part -- C^3(x). Take only those parts of degree == max */
	for( i=0; i<=max; i++ )
		for( j=0; j<=max-i; j++ ) {
			mpz_mul (temp, source [i], source [j]);
			mpz_mul (temp, source [max-i-j], temp);
			mpz_add (factor, temp, factor);
		}
	/* next part -- 3*C(x)*C(x^2). A bit tricky... */
	for( i=max%2; i<=max; i+=2 ) {			/* evens and odds */
		mpz_mul_ui (temp, source [i], 3);
		mpz_mul (temp, source [(max-i)>>1], temp);
		mpz_add (factor, temp, factor);
	}
	/* last part -- 2*C(x^3). */
	if( !(max % 3) ) {
		mpz_mul_ui (temp, source [max/3], 2);
		mpz_add (factor, temp, factor);
	}
	/* and finally, divide factor by 6. If it doesn't........ */
	mpz_tdiv_qr_ui (factor, temp, factor, 6);
	if( mpz_cmp_ui (temp, 0) ) {
		fprintf (stderr, "INTERNAL ERROR OR BAD DATA!\n");
		fprintf (stderr, " is not divisible by 6 (#%d)\n", max+1);
		exit (3);				/* ...... */
	}
#ifdef TIMER
	timer = alk_time () - timer;
#endif
}

void print (void) {
	if( !(iofile = fopen (filename, "a")) ) {
		fprintf (stderr, "Error: cannot open output file <%s>\n", filename);
		fclose (iofile);
		exit (1);
	}

#ifdef TIMER
	if( !(logfile = fopen (logname, "a")) ) {
		fprintf (stderr, "Warning: cannot open log file <%s>\n", logname);
	} else {
		fprintf (logfile, LOGFORMAT, max+1, timer);
		putc ('\n', logfile);
	}
#endif
#ifdef DEBUG
	fprintf (stderr, "Writing to <%s>.\n", filename);
#endif

	fprintf (iofile, FORMAT, max+1);
	mpz_out_str (iofile, 10, factor);
	putc ('\n', iofile);
}

void shutdown (void) {
	fclose (iofile);
#ifdef TIMER
	fclose (logfile);
#endif
}

Re: [room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Andrey Rahmatullin]

[-- Attachment #1: Type: text/plain, Size: 318 bytes --]

On Thu, Nov 03, 2005 at 11:45:43PM +0200, Michael Shigorin wrote:
> an attach.  Цепляю ещё раз (хотя и тогда вроде не забыл?).
Забыл ;)

-- 
WBR, wRAR (ALT Linux Team)
Powered by the ALT Linux fortune(8):

<raorn> наш jpegtran не корёжит exif, родной (или таки бсдшный?) - корёжит
<gvy> raorn, вешай багу :)

[-- Attachment #2: Digital signature --]
[-- Type: application/pgp-signature, Size: 189 bytes --]

Re: [room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Aleksey Korotkov]

Вот вам MATLAB. И даже open-source :)

http://www.nixp.ru/news/6769

-- 
С уважением,
Алексей Коротков

mailto:ziga@uni.udm.ru

Re: [room] Re: [JT] Re: [Comm] Re: MathCAD, MathLab

[Andrey Rahmatullin]

[-- Attachment #1: Type: text/plain, Size: 436 bytes --]

On Sat, Nov 05, 2005 at 03:23:02AM +0400, Aleksey Korotkov wrote:
> Вот вам MATLAB. И даже open-source :)
Где видно, что опенсорс?
http://www.nixp.ru/cdrom/req/306

-- 
WBR, wRAR (ALT Linux Team)
Powered by the ALT Linux fortune(8):

Моим преподавателем астрономии был проф. М.М.Дагаев, в конце 70-х ему
было уже, наверное, за 70, он курил "Казбек" и ненавидел постоянную
Хаббла за ее непостоянство.
		-- aen in talk-room@

[-- Attachment #2: Digital signature --]
[-- Type: application/pgp-signature, Size: 189 bytes --]

[room] Re: MathLab

[Aleksey Korotkov]

On Sat, 5 Nov 2005 12:43:59 +0500
Andrey Rahmatullin wrote:

AR> Где видно, что опенсорс?

Так по ссылке, которую я приводил, заглядывали?
http://www.nixp.ru/news/6769

Цитирую:

	В разделе распространения дисков с open-source ПО на
	nixp.ru пополнение: Scientific Linux 4.1 i386 (4 CD, 200
	рублей); Matlab R14 for Linux/Unix/MacOSX (3 CD, 300
	рублей); OpenOffice 2.0 RUS (1 CD, 50 рублей).

Читать внимательно первую строчку.

-- 
С уважением,
Алексей Коротков

mailto:ziga@uni.udm.ru

Re: [room] Re: MathLab

[Michael Shigorin]

On Sun, Nov 06, 2005 at 01:54:27AM +0400, Aleksey Korotkov wrote:
> AR> Где видно, что опенсорс?
> Так по ссылке, которую я приводил, заглядывали?
> http://www.nixp.ru/news/6769
> Цитирую:
> 	В разделе распространения дисков с open-source ПО на
> 
> Читать внимательно первую строчку.

Угу, там ещё Oracle, America's Army, HoMM3 неподалёку.
Боюсь, для выводопостроения сильно по касательной.

-- 
 ---- WBR, Michael Shigorin <mike@altlinux.ru>
  ------ Linux.Kiev http://www.linux.kiev.ua/