static char rcsid[] = "$Header: mkfs.c,v 10.0 86/07/21 14:52:54 bog Exp $";
static char sccsid[]="%W% %Y% %Q% %G%";

/*
 * make file system for cylinder-group style file systems
 *
#ifndef STANDALONE
 * usage: mkfs special size [ nsect ntrak bsize fsize cpg minfree rps nbpi ]
#else   STANDALONE
 * Boot this from CRM and it will recreate file systems on your disk..how nice!
#endif  STANDALONE
 */

/*
 * The following constants set the defaults used for the number
 * of sectors (fs_nsect), and number of tracks (fs_ntrak).
 */
#define DFLNSECT	32
#define DFLNTRAK	16

/*
 * The following two constants set the default block and fragment sizes.
 * Both constants must be a power of 2 and meet the following constraints:
 *	MINBSIZE <= DESBLKSIZE <= MAXBSIZE
 *	DEV_BSIZE <= DESFRAGSIZE <= DESBLKSIZE
 *	DESBLKSIZE / DESFRAGSIZE <= 8
 */
#define DESBLKSIZE	4096
#define DESFRAGSIZE	512

/*
 * Cylinder groups may have up to MAXCPG cylinders. The actual
 * number used depends upon how much information can be stored
 * on a single cylinder. The default is to used 16 cylinders
 * per group.
 */
#define	DESCPG		16	/* desired fs_cpg */

/*
 * MINFREE gives the minimum acceptable percentage of file system
 * blocks which may be free. If the freelist drops below this level
 * only the superuser may continue to allocate blocks. This may
 * be set to 0 if no reserve of free blocks is deemed necessary,
 * however throughput drops by fifty percent if the file system
 * is run at between 90% and 100% full; thus the default value of
 * fs_minfree is 10%.
 */
#ifdef s32
#define MINFREE		4
#else  s32
#define MINFREE		10
#endif s32

/*
 * ROTDELAY gives the minimum number of milliseconds to initiate
 * another disk transfer on the same cylinder. It is used in
 * determining the rotationally optimal layout for disk blocks
 * within a file; the default of fs_rotdelay is 4ms.
 */
#define ROTDELAY	2

/*
 * MAXCONTIG sets the default for the maximum number of blocks
 * that may be allocated sequentially. Since UNIX drivers are
 * not capable of scheduling multi-block transfers, this defaults
 * to 1 (ie no contiguous blocks are allocated).
 */
#define MAXCONTIG	1

/*
 * MAXBLKPG determines the maximum number of data blocks which are
 * placed in a single cylinder group. This is currently a function
 * of the block and fragment size of the file system.
 */
#define MAXBLKPG(fs)	((fs)->fs_fsize / sizeof(daddr_t))

/*
 * Each file system has a number of inodes statically allocated.
 * We allocate one inode slot per NBPI bytes, expecting this
 * to be far more than we will ever need.
 */
#define	NBPI		1536

/*
 * Disks are assumed to rotate at 60HZ, unless otherwise specified.
 */
#define	DEFHZ		60

/*
 * STANDALONE seems to cause everything to be done by default numbers -- 
 * something we definitely do not want. The postulated Valid standalone
 * mkfs/format/install/update package definitely will take parametrized
 * arguments and bandy them about for various disks. I think.
 *
 * Historical note for those who follow: at this point, we are obviously hazy
 * on the functionality and implementation of said standalone whizzy thing.
 *
 * (bog 0586)
 *
 * New understanding:  The ``meta disktab'', a file describing disks used in
 * Valid systems in a high-level fashion, is now used to generate:
 *	1. /etc/disktab
 *	2. /etc/dkconfig
 *	3. /usr/include/sa/dk_types.h  ( Maybe a different path name )
 * These three files are detailed descriptions of disk parameters.
 * Dk_types.h contains data structures, dk_names and dk_models, defining all
 * disks.  mkfs uses this data structure.
 */
#ifndef STANDALONE
#include <stdio.h>
#include <a.out.h>
#else   not STANDALONE

#define	STRUCT_INIT		/* Initialize dk_types structures. */
#include "dk_types.h"		/* For disk types structures. */

#ifndef FORMAT_FROM_DKMODELS
#include "../h/param.h"
#include "../h/fs.h"
#include "../h/inode.h"
#include "../salib/saio.h"	/* For ioctls */
#include "../s32/dklabel.h"
#endif FORMAT_FROM_DKMODELS

#include "../h/time.h"	/* For getting the time.      */
#endif  not STANDALONE

#include <sys/param.h>
#include <sys/inode.h>
#include <sys/fs.h>
#include <sys/dir.h>

#define UMASK		0755
#define MAXINOPB	(MAXBSIZE / sizeof(struct dinode))
#define POWEROF2(num)	(((num) & ((num) - 1)) == 0)

union {
	struct fs fs;
	char pad[MAXBSIZE];
} fsun;
#define	sblock	fsun.fs

#ifndef STANDALONE
struct	csum *fscs;
#else   STANDALONE
/*
 * ncg = ncyl / cpg.  cpg <= MAXCPG.
 * Say the max ncyl we're likely to get is 5000 (we're UNlikely to get that,
 * but what the hell). To be paranoid, let's say we might get 4000 struct
 * csums in a MONSTER disk.  Or one with a heinous set of frag, block sizes...
 *
 * This leaves us with 16000 bytes for the fscs[] structure.  What the hell.
 */
#define MAXCSUMS	4000
struct csum fscs[MAXCSUMS];
/*
 * Version string comes from vers.c -- created by vers.sh on the fly.
 */
extern char version[];
#endif  STANDALONE

union {
	struct cg cg;
	char pad[MAXBSIZE];
} cgun;
#define	acg	cgun.cg

struct	dinode zino[MAXIPG];

#ifdef	STANDALONE
char	fsys[LITTLEBSIZE];
#else	STANDALONE
char	*fsys;
#endif	STANDALONE
time_t	uutime;	/* For timestamping the file system. */
int	fsi;
int	fso;
daddr_t	alloc();

#ifdef STANDALONE
main()
#else  STANDALONE
main(argc, argv)
	int argc;
	char *argv[];
#endif STANDALONE
{
	long cylno, rpos, blk, i, j, inos, fssize, warn = 0;
#ifdef STANDALONE
	/*
	 * Extra variables for standalone use.
	 */
	struct timeval tval;	/* To get time of day. */
	register struct model *disk, *label;
	int part;
#endif STANDALONE

#ifndef STANDALONE
	argc--, argv++;
	time(&uutime);
	if (argc < 2) {
		printf("usage: mkfs special size [ nsect ntrak bsize fsize cpg minfree rps nbpi ]\n");
		exit(1);
	}
	fsys = argv[0];
	fssize = atoi(argv[1]);
	fso = creat(fsys, 0666);
	if(fso < 0) {
		printf("%s: cannot create\n", fsys);
		exit(1);
	}
	fsi = open(fsys, 0);
	if(fsi < 0) {
		printf("%s: cannot open\n", fsys);
		exit(1);
	}
#else   not STANDALONE
#ifndef notdef
	printf(version);
#else  notdef
	printf("\nMkfs\n");
#endif notdef
	gettimeofday(&tval);
	uutime = tval.tv_sec;
	do
	{
		printf("- Possible disks -\n");
		dump_dkmodels();
		printf("\n\nWhat sort of disk do you wish to format? ");
		gets(fsys);
		for (disk = dk_models; (disk) ; disk = disk->m_next) {
			register name_list_t *nl;
			
			for (nl = disk->m_name; nl && strcmp(nl->nl_name,fsys);
				nl = nl->nl_next) ;
			if (nl) break;
		}
		if (disk == 0)
		{
			printf("Invalid disk type %s\n",fsys);
			fso = fsi = -1;
		}
		else
		{
			/* 
			 * Using the Valid standalone library, you would open
			 * the partition you wish to do a mkfs upon, E.G.
			 * ``id(1,6)'' would be opened if you wished to create
			 * a file system on partition 6 of disk 1.
			 */
			printf("Unit to create file system on (ex. id(0,6)): ");
			gets(fsys);
			fso = open(fsys, 1);	/* write */
			fsi = open(fsys, 0);	/* read  */
		}
	} while (fso < 0 || fsi < 0);
	part = atoi(index(fsys,',')+1); /* The number after the comma. */
	/*
	 * Slap label parameters on top of retrieved disk model,
	 * replacing possibly wrong things like partition sizes.
	 */
	if (ioctl(fsi, SAIOGPARM, label))
		exit(perror(fsys));
	if (label->m_nCyl != disk->m_nCyl || label->m_aCyl != disk->m_aCyl ||
	    label->m_nHead != disk->m_nHead || 
	    label->m_nSectors != disk->m_nSectors || 
	    label->m_Gap1 != disk->m_Gap1 || label->m_Gap2 != disk->m_Gap2) {
		printf ("label->m_nCyl(%d) != disk->m_nCyl(%d)\n",
			label->m_nCyl, disk->m_nCyl);
		printf ("label->m_aCyl(%d) != disk->m_aCyl(%d)\n",
			label->m_aCyl, disk->m_aCyl);
    		printf ("label->m_nHead(%d) != disk->m_nHead(%d)n",
			label->m_nHead, disk->m_nHead);
    		printf ("label->m_nSectors(%d) != disk->m_nSectors(%d)\n",
			label->m_nSectors, disk->m_nSectors);
    		printf ("label->m_Gap1(%d) != disk->m_Gap1(%d)\n",
			label->m_Gap1, disk->m_Gap1);
		printf ("label->m_Gap2(%d) != disk->m_Gap2(%d)\n",
			label->m_Gap2, disk->m_Gap2);

		exit (printf ("Label mismatch error, reformat\n"));
	}
	for (i = 0; i < DK_NSIZES; i++)
	{
		/* If the partition size is 0 then we don't care
		 * about the start partition 
		 *
		 * KLUDGE for XCRM 			12/4/86	  BW
		 */
		if( label->m_PartSizes[i].sz_nblk == 0 )
			continue;


		if (label->m_PartSizes[i].sz_cyloff != 
		    disk->m_PartSizes[i].sz_cyloff ||
		    label->m_PartSizes[i].sz_nblk !=
		    disk->m_PartSizes[i].sz_nblk)
			break;
	}
	if (i != DK_NSIZES)
		if(strcmp ("eagle", disk->m_name->nl_name))
			exit (printf ("Label partition size error\n"));
		else {
			struct dk_label lab;

			if (ioctl (fsi, SAIORLAB, &lab))
				exit (perror (fsys));
			for (i = 0; i < DK_NSIZES; i++) {
				lab.dkl_sizes[i].sz_cyloff = 
				    disk->m_PartSizes[i].sz_cyloff;
				lab.dkl_sizes[i].sz_nblk =
				    disk->m_PartSizes[i].sz_nblk;
			}
			if (ioctl (fso, SAIOSPARM, disk))
				exit (perror (fsys));
			if (ioctl (fso, SAIOWLAB, &lab))
				exit (printf("Can't update eagle label\n"));
		}
			
#ifdef	DEBUG
	printf("Partition 0x%x:\n",part);
	dump_dkname(disk);
	dump_disk(disk);
#endif	DEBUG
	fssize = disk->m_PartSizes[part].sz_nblk;
#endif  not STANDALONE
	/*
	 * Validate the given file system size.
	 * Verify that its last block can actually be accessed.
	 */
	if (fssize <= 0)
	{
		printf("Negative or Zero size %d\n", fssize);
		exit(1);
	}
	/* 
	 * Write the last block.
	 * Sblock here is only a handy source of bytes.
	 * wtfs() calls exit() directly if it cannot write.
	 */
	wtfs(fssize - 1, DEV_BSIZE, (char *)&sblock);
	/*
	 * collect and verify the sector and track info
	 */
#ifndef STANDALONE
	if (argc > 2)
		sblock.fs_nsect = atoi(argv[2]);
	else
		sblock.fs_nsect = DFLNSECT;
	if (argc > 3)
		sblock.fs_ntrak = atoi(argv[3]);
	else
		sblock.fs_ntrak = DFLNTRAK;
#else  STANDALONE
	sblock.fs_nsect = disk->m_nSectors;
	sblock.fs_ntrak = disk->m_nHead;
#endif STANDALONE
	if (sblock.fs_ntrak <= 0)
		printf("negative/zero ntrak %d.\n", sblock.fs_ntrak), exit(1);
	if (sblock.fs_nsect <= 0)
		printf("negative/zero nsect %d.\n", sblock.fs_nsect), exit(1);
	sblock.fs_spc = sblock.fs_ntrak * sblock.fs_nsect;
	/*
	 * collect and verify the block and fragment sizes
	 */
#ifndef STANDALONE
	if (argc > 4)
		sblock.fs_bsize = atoi(argv[4]);
	else
		sblock.fs_bsize = DESBLKSIZE;
	if (argc > 5)
		sblock.fs_fsize = atoi(argv[5]);
	else
		sblock.fs_fsize = DESFRAGSIZE;
#else   STANDALONE
	sblock.fs_bsize = disk->m_FSBlockSize;
	sblock.fs_fsize = disk->m_FragSize;
#endif  STANDALONE
	if (!POWEROF2(sblock.fs_bsize)) {
		printf("block size must be a power of 2, and %d isn't.\n",
		    sblock.fs_bsize);
		exit(1);
	}
	if (!POWEROF2(sblock.fs_fsize)) {
		printf("fragment size must be a power of 2, and %d isn't.\n",
		    sblock.fs_fsize);
		exit(1);
	}
	if (sblock.fs_fsize < DEV_BSIZE) {
		printf("fragment size %d is too small, minimum is %d\n",
		    sblock.fs_fsize, DEV_BSIZE);
		exit(1);
	}
	if (sblock.fs_bsize < MINBSIZE) {
		printf("block size %d is too small, minimum is %d\n",
		    sblock.fs_bsize, MINBSIZE);
		exit(1);
	}
	if (sblock.fs_bsize < sblock.fs_fsize) {
		printf("block size (%d) cannot be smaller than fragment size (%d)\n",
		    sblock.fs_bsize, sblock.fs_fsize);
		exit(1);
	}
	sblock.fs_bmask = ~(sblock.fs_bsize - 1);
	sblock.fs_fmask = ~(sblock.fs_fsize - 1);
	for (sblock.fs_bshift = 0, i = sblock.fs_bsize; i > 1; i >>= 1)
		sblock.fs_bshift++;
	for (sblock.fs_fshift = 0, i = sblock.fs_fsize; i > 1; i >>= 1)
		sblock.fs_fshift++;
	sblock.fs_frag = numfrags(&sblock, sblock.fs_bsize);
	for (sblock.fs_fragshift = 0, i = sblock.fs_frag; i > 1; i >>= 1)
		sblock.fs_fragshift++;
	if (sblock.fs_frag > MAXFRAG) {
		printf("fragment size %d is too small, minimum with block size %d is %d\n",
		    sblock.fs_fsize, sblock.fs_bsize,
		    sblock.fs_bsize / MAXFRAG);
		exit(1);
	}
	sblock.fs_nindir = sblock.fs_bsize / sizeof(daddr_t);
	sblock.fs_inopb = sblock.fs_bsize / sizeof(struct dinode);
	sblock.fs_nspf = sblock.fs_fsize / DEV_BSIZE;
	for (sblock.fs_fsbtodb = 0, i = sblock.fs_nspf; i > 1; i >>= 1)
		sblock.fs_fsbtodb++;
	sblock.fs_sblkno =
	    roundup(howmany(BBSIZE + SBSIZE, sblock.fs_fsize), sblock.fs_frag);
	sblock.fs_cblkno = (daddr_t)(sblock.fs_sblkno +
	    roundup(howmany(SBSIZE, sblock.fs_fsize), sblock.fs_frag));
	sblock.fs_iblkno = sblock.fs_cblkno + sblock.fs_frag;
	sblock.fs_cgoffset = roundup(
	    howmany(sblock.fs_nsect, sblock.fs_fsize / DEV_BSIZE),
	    sblock.fs_frag);
#ifndef notdef
	/* ``word full of ones'' = ffffffff is word length dependent. */
	for (sblock.fs_cgmask = 0xffffffff, i = sblock.fs_ntrak; i > 1; i >>= 1)
		sblock.fs_cgmask <<= 1;
#else   notdef
	/* ~0 is not word length dependent. */
	for (sblock.fs_cgmask = ~(0), i = sblock.fs_ntrak; i > 1; i >>= 1)
		sblock.fs_cgmask <<= 1;
#endif  notdef
	if (!POWEROF2(sblock.fs_ntrak))
		sblock.fs_cgmask <<= 1;
	for (sblock.fs_cpc = NSPB(&sblock), i = sblock.fs_spc;
	     sblock.fs_cpc > 1 && (i & 1) == 0;
	     sblock.fs_cpc >>= 1, i >>= 1)
		/* void */;
	if (sblock.fs_cpc > MAXCPG) {
		printf("maximum block size with nsect %d and ntrak %d is %d\n",
		    sblock.fs_nsect, sblock.fs_ntrak,
		    sblock.fs_bsize / (sblock.fs_cpc / MAXCPG));
		exit(1);
	}
	/* 
	 * collect and verify the number of cylinders per group
	 */
#ifndef STANDALONE
	if (argc > 6) {
		sblock.fs_cpg = atoi(argv[6]);
		sblock.fs_fpg = (sblock.fs_cpg * sblock.fs_spc) / NSPF(&sblock);
	} else {
		sblock.fs_cpg = MAX(sblock.fs_cpc, DESCPG);
		sblock.fs_fpg = (sblock.fs_cpg * sblock.fs_spc) / NSPF(&sblock);
		while (sblock.fs_fpg / sblock.fs_frag > MAXBPG(&sblock) &&
		    sblock.fs_cpg > sblock.fs_cpc) {
			sblock.fs_cpg -= sblock.fs_cpc;
			sblock.fs_fpg =
			    (sblock.fs_cpg * sblock.fs_spc) / NSPF(&sblock);
		}
	}
#else  STANDALONE
	sblock.fs_cpg = MAX(sblock.fs_cpc, DESCPG);
	sblock.fs_fpg = (sblock.fs_cpg * sblock.fs_spc) / NSPF(&sblock);
	while (sblock.fs_fpg / sblock.fs_frag > MAXBPG(&sblock) &&
		sblock.fs_cpg > sblock.fs_cpc) {
		sblock.fs_cpg -= sblock.fs_cpc;
		sblock.fs_fpg = (sblock.fs_cpg * sblock.fs_spc) / NSPF(&sblock);
	}
#endif  STANDALONE
	if (sblock.fs_cpg < 1) {
		printf("cylinder groups must have at least 1 cylinder\n");
		exit(1);
	}
	if (sblock.fs_cpg > MAXCPG) {
		printf("cylinder groups are limited to %d cylinders\n", MAXCPG);
		exit(1);
	}
	if (sblock.fs_cpg % sblock.fs_cpc != 0) {
		printf("cylinder groups must have a multiple of %d cylinders\n",
		    sblock.fs_cpc);
		exit(1);
	}
	/*
	 * Now have size for file system and nsect and ntrak.
	 * Determine number of cylinders and blocks in the file system.
	 */
	sblock.fs_size = fssize = dbtofsb(&sblock, fssize);
	sblock.fs_ncyl = fssize * NSPF(&sblock) / sblock.fs_spc;
	if (fssize * NSPF(&sblock) > sblock.fs_ncyl * sblock.fs_spc) {
		sblock.fs_ncyl++;
		warn = 1;
	}
	if (sblock.fs_ncyl < 1) {
		printf("file systems must have at least one cylinder\n");
		exit(1);
	}
	/*
	 * determine feasability/values of rotational layout tables
	 */
	if (sblock.fs_ntrak == 1) {
		sblock.fs_cpc = 0;
		goto next;
	}
	if (sblock.fs_spc * sblock.fs_cpc > MAXBPC * NSPB(&sblock) ||
	    sblock.fs_nsect > (1 << NBBY) * NSPB(&sblock)) {
		printf("%s %s %d %s %d.%s",
		    "Warning: insufficient space in super block for\n",
		    "rotational layout tables with nsect", sblock.fs_nsect,
		    "and ntrak", sblock.fs_ntrak,
		    "\nFile system performance may be impared.\n");
		sblock.fs_cpc = 0;
		goto next;
	}
	/*
	 * calculate the available blocks for each rotational position
	 */
	for (cylno = 0; cylno < MAXCPG; cylno++)
		for (rpos = 0; rpos < NRPOS; rpos++)
			sblock.fs_postbl[cylno][rpos] = -1;
	blk = sblock.fs_spc * sblock.fs_cpc / NSPF(&sblock);
	for (i = 0; i < blk; i += sblock.fs_frag)
		/* void */;
	for (i -= sblock.fs_frag; i >= 0; i -= sblock.fs_frag) {
		cylno = cbtocylno(&sblock, i);
		rpos = cbtorpos(&sblock, i);
		blk = i / sblock.fs_frag;
		if (sblock.fs_postbl[cylno][rpos] == -1)
			sblock.fs_rotbl[blk] = 0;
		else
			sblock.fs_rotbl[blk] =
			    sblock.fs_postbl[cylno][rpos] - blk;
		sblock.fs_postbl[cylno][rpos] = blk;
	}
next:
	/*
	 * Validate specified/determined cpg.
	 */
	if (sblock.fs_spc > MAXBPG(&sblock) * NSPB(&sblock)) {
		printf("too many sectors per cylinder (%d sectors)\n",
		    sblock.fs_spc);
		while(sblock.fs_spc > MAXBPG(&sblock) * NSPB(&sblock)) {
			sblock.fs_bsize <<= 1;
			if (sblock.fs_frag < MAXFRAG)
				sblock.fs_frag <<= 1;
			else
				sblock.fs_fsize <<= 1;
		}
		printf("nsect %d, and ntrak %d, requires block size of %d,\n",
		    sblock.fs_nsect, sblock.fs_ntrak, sblock.fs_bsize);
		printf("\tand fragment size of %d\n", sblock.fs_fsize);
		exit(1);
	}
	if (sblock.fs_fpg > MAXBPG(&sblock) * sblock.fs_frag) {
		printf("cylinder group too large (%d cylinders); ",
		    sblock.fs_cpg);
		printf("max: %d cylinders per group\n",
		    MAXBPG(&sblock) * sblock.fs_frag /
		    (sblock.fs_fpg / sblock.fs_cpg));
		exit(1);
	}
	sblock.fs_cgsize = fragroundup(&sblock,
	    sizeof(struct cg) + howmany(sblock.fs_fpg, NBBY));
	/*
	 * Compute/validate number of cylinder groups.
	 */
	sblock.fs_ncg = sblock.fs_ncyl / sblock.fs_cpg;
	if (sblock.fs_ncyl % sblock.fs_cpg)
		sblock.fs_ncg++;
	if ((sblock.fs_spc * sblock.fs_cpg) % NSPF(&sblock)) {
		printf("mkfs: nsect %d, ntrak %d, cpg %d is not tolerable\n",
		    sblock.fs_nsect, sblock.fs_ntrak, sblock.fs_cpg);
		printf("as this would would have cyl groups whose size\n");
		printf("is not a multiple of %d; choke!\n", sblock.fs_fsize);
		exit(1);
	}
	/*
	 * Compute number of inode blocks per cylinder group.
	 * Start with one inode per NBPI bytes; adjust as necessary.
	 */
	inos = MAX(NBPI, sblock.fs_fsize);
#ifndef STANDALONE
	if (argc > 9) {
		i = atoi(argv[9]);
		if (i <= 0)
			printf("%s: bogus nbpi reset to %d\n", argv[9], inos);
		else
			inos = i;
	}
#else   STANDALONE
	/*
	 * We don't support changing # inodes per cylinder group
	 * from Standalone mode - yet.
	 */
#endif  STANDALONE
	i = sblock.fs_iblkno + MAXIPG / INOPF(&sblock);
	inos = (fssize - sblock.fs_ncg * i) * sblock.fs_fsize / inos /
	    INOPB(&sblock);
	if (inos <= 0)
		inos = 1;
	sblock.fs_ipg = ((inos / sblock.fs_ncg) + 1) * INOPB(&sblock);
	if (sblock.fs_ipg > MAXIPG)
		sblock.fs_ipg = MAXIPG;
	sblock.fs_dblkno = sblock.fs_iblkno + sblock.fs_ipg / INOPF(&sblock);
	i = MIN(~sblock.fs_cgmask, sblock.fs_ncg - 1);
	if (cgdmin(&sblock, i) - cgbase(&sblock, i) >= sblock.fs_fpg) {
		printf("inode blocks/cyl group (%d) >= data blocks (%d)\n",
		    cgdmin(&sblock, i) - cgbase(&sblock, i) / sblock.fs_frag,
		    sblock.fs_fpg / sblock.fs_frag);
		printf("number of cylinders per cylinder group must be increased\n");
		exit(1);
	}
	j = sblock.fs_ncg - 1;
	if ((i = fssize - j * sblock.fs_fpg) < sblock.fs_fpg &&
	    cgdmin(&sblock, j) - cgbase(&sblock, j) > i) {
		printf("Warning: inode blocks/cyl group (%d) >= data blocks (%d) in last\n",
		    (cgdmin(&sblock, j) - cgbase(&sblock, j)) / sblock.fs_frag,
		    i / sblock.fs_frag);
		printf("    cylinder group. This implies %d sector(s) cannot be allocated.\n",
		    i * NSPF(&sblock));
		sblock.fs_ncg--;
		sblock.fs_ncyl -= sblock.fs_ncyl % sblock.fs_cpg;
		sblock.fs_size = fssize = sblock.fs_ncyl * sblock.fs_spc /
		    NSPF(&sblock);
		warn = 0;
	}
	if (warn) {
		printf("Warning: %d sector(s) in last cylinder unallocated\n",
		    sblock.fs_spc -
		    (fssize * NSPF(&sblock) - (sblock.fs_ncyl - 1)
		    * sblock.fs_spc));
	}
	/*
	 * fill in remaining fields of the super block
	 */
	sblock.fs_csaddr = cgdmin(&sblock, 0);
	/*
	 * fragroundup(fs, size) calculates roundup(size, fs->fs_fsize).
	 * fs_fsize >= DEV_BSIZE, and must be a power of 2.
	 */
	sblock.fs_cssize =
	    fragroundup(&sblock, sblock.fs_ncg * sizeof(struct csum));
	i = sblock.fs_bsize / sizeof(struct csum);
	sblock.fs_csmask = ~(i - 1);
	for (sblock.fs_csshift = 0; i > 1; i >>= 1)
		sblock.fs_csshift++;
	i = sizeof(struct fs) +
		howmany(sblock.fs_spc * sblock.fs_cpc, NSPB(&sblock));
	sblock.fs_sbsize = fragroundup(&sblock, i);
#ifndef STANDALONE
	fscs = (struct csum *)calloc(1, sblock.fs_cssize);
#else   STANDALONE
	if (sblock.fs_cssize > MAXCSUMS * sizeof(struct csum))
	{
		printf("fs_cssize 0x%x > sizeof(fscs) 0x%x\n",
			sblock.fs_cssize, sizeof(fscs));
		exit(1);
	}
#endif  STANDALONE
	sblock.fs_magic = FS_MAGIC;
	sblock.fs_rotdelay = ROTDELAY;
#ifndef STANDALONE
	if (argc > 7) {
		sblock.fs_minfree = atoi(argv[7]);
		if (sblock.fs_minfree < 0 || sblock.fs_minfree > 99) {
			printf("%s: bogus minfree reset to %d%%\n", argv[7],
				MINFREE);
			sblock.fs_minfree = MINFREE;
		}
	} else
		sblock.fs_minfree = MINFREE;
#else   STANDALONE
	sblock.fs_minfree = MINFREE;
#endif  STANDALONE
	sblock.fs_maxcontig = MAXCONTIG;
	sblock.fs_maxbpg = MAXBLKPG(&sblock);
#ifndef STANDALONE
	if (argc > 8)
		sblock.fs_rps = atoi(argv[8]);
	else
		sblock.fs_rps = DEFHZ;
#else   STANDALONE
	sblock.fs_rps = 60;
#endif  STANDALONE
	sblock.fs_cgrotor = 0;
	sblock.fs_cstotal.cs_ndir = 0;
	sblock.fs_cstotal.cs_nbfree = 0;
	sblock.fs_cstotal.cs_nifree = 0;
	sblock.fs_cstotal.cs_nffree = 0;
	sblock.fs_fmod = 0;
	sblock.fs_ronly = 0;
	/*
	 * Dump out summary information about file system.
	 */
	printf("%s:\t%d sectors in %d cylinders of %d tracks, %d sectors\n",
	    fsys, sblock.fs_size * NSPF(&sblock), sblock.fs_ncyl,
	    sblock.fs_ntrak, sblock.fs_nsect);
#ifdef	STANDALONE
	printf("\t%d.%dMb in %d cyl groups (%d c/g, %d.%dMb/g, %d i/g)\n",
	    (sblock.fs_size * sblock.fs_fsize) / 1000000,
	    ((sblock.fs_size * sblock.fs_fsize) % 1000000) / 100000,
	    sblock.fs_ncg, sblock.fs_cpg,
	    (sblock.fs_fpg * sblock.fs_fsize) / 1000000,
	    ((sblock.fs_fpg * sblock.fs_fsize) / 1000000) / 100000,
	    sblock.fs_ipg);
#else	STANDALONE
	printf("\t%.1fMb in %d cyl groups (%d c/g, %.2fMb/g, %d i/g)\n",
	    (float)sblock.fs_size * sblock.fs_fsize * 1e-6, sblock.fs_ncg,
	    sblock.fs_cpg, (float)sblock.fs_fpg * sblock.fs_fsize * 1e-6,
	    sblock.fs_ipg);
#endif	STANDALONE
	/*
	 * Now build the cylinders group blocks and
	 * then print out indices of cylinder groups.
	 */
	printf("super-block backups (for fsck -b#) at:");
	for (cylno = 0; cylno < sblock.fs_ncg; cylno++) {
		initcg(cylno);
		if (cylno % 10 == 0)
			printf("\n");
		printf(" %d,", fsbtodb(&sblock, cgsblock(&sblock, cylno)));
	}
	printf("\n");
	/*
	 * Now construct the initial file system,
	 * then write out the super-block.
	 */
	fsinit();
	sblock.fs_time = uutime;
	wtfs(SBLOCK, SBSIZE, (char *)&sblock);
	for (i = 0; i < sblock.fs_cssize; i += sblock.fs_bsize)
		wtfs(fsbtodb(&sblock, sblock.fs_csaddr + numfrags(&sblock, i)),
			sblock.fs_cssize - i < sblock.fs_bsize ?
			    sblock.fs_cssize - i : sblock.fs_bsize,
			((char *)fscs) + i);
	/* 
	 * Write out the duplicate super blocks
	 */
	for (cylno = 0; cylno < sblock.fs_ncg; cylno++)
		wtfs(fsbtodb(&sblock, cgsblock(&sblock, cylno)),
		    SBSIZE, (char *)&sblock);
#ifndef STANDALONE
	exit(0);
#endif  not STANDALONE
}

/*
 * Initialize a cylinder group.
 */
initcg(cylno)
	int cylno;
{
	daddr_t cbase, d, dlower, dupper, dmax;
	long i, j, s;
	register struct csum *cs;

	/*
	 * Determine block bounds for cylinder group.
	 * Allow space for super block summary information in first
	 * cylinder group.
	 */
	cbase = cgbase(&sblock, cylno);
	dmax = cbase + sblock.fs_fpg;
	if (dmax > sblock.fs_size)
		dmax = sblock.fs_size;
	dlower = cgsblock(&sblock, cylno) - cbase;
	dupper = cgdmin(&sblock, cylno) - cbase;
	cs = fscs + cylno;
	acg.cg_time = uutime;
	acg.cg_magic = CG_MAGIC;
	acg.cg_cgx = cylno;
	if (cylno == sblock.fs_ncg - 1)
		acg.cg_ncyl = sblock.fs_ncyl % sblock.fs_cpg;
	else
		acg.cg_ncyl = sblock.fs_cpg;
	acg.cg_niblk = sblock.fs_ipg;
	acg.cg_ndblk = dmax - cbase;
	acg.cg_cs.cs_ndir = 0;
	acg.cg_cs.cs_nffree = 0;
	acg.cg_cs.cs_nbfree = 0;
	acg.cg_cs.cs_nifree = 0;
	acg.cg_rotor = 0;
	acg.cg_frotor = 0;
	acg.cg_irotor = 0;
	for (i = 0; i < sblock.fs_frag; i++) {
		acg.cg_frsum[i] = 0;
	}
	for (i = 0; i < sblock.fs_ipg; ) {
		for (j = INOPB(&sblock); j > 0; j--) {
			clrbit(acg.cg_iused, i);
			i++;
		}
		acg.cg_cs.cs_nifree += INOPB(&sblock);
	}
	if (cylno == 0)
		for (i = 0; i < ROOTINO; i++) {
			setbit(acg.cg_iused, i);
			acg.cg_cs.cs_nifree--;
		}
	while (i < MAXIPG) {
		clrbit(acg.cg_iused, i);
		i++;
	}
	lseek(fso, fsbtodb(&sblock, cgimin(&sblock, cylno)) * DEV_BSIZE, 0);
	if (write(fso, (char *)zino, sblock.fs_ipg * sizeof (struct dinode)) !=
	    sblock.fs_ipg * sizeof (struct dinode))
		printf("write error 0x%x\n", numfrags(&sblock, tell(fso)));
	for (i = 0; i < MAXCPG; i++) {
		acg.cg_btot[i] = 0;
		for (j = 0; j < NRPOS; j++)
			acg.cg_b[i][j] = 0;
	}
	if (cylno == 0) {
		/*
		 * reserve space for summary info and Boot block
		 */
		dupper += howmany(sblock.fs_cssize, sblock.fs_fsize);
		for (d = 0; d < dlower; d += sblock.fs_frag)
			clrblock(&sblock, acg.cg_free, d/sblock.fs_frag);
	} else {
		for (d = 0; d < dlower; d += sblock.fs_frag) {
			setblock(&sblock, acg.cg_free, d/sblock.fs_frag);
			acg.cg_cs.cs_nbfree++;
			acg.cg_btot[cbtocylno(&sblock, d)]++;
			acg.cg_b[cbtocylno(&sblock, d)][cbtorpos(&sblock, d)]++;
		}
		sblock.fs_dsize += dlower;
	}
	sblock.fs_dsize += acg.cg_ndblk - dupper;
	for (; d < dupper; d += sblock.fs_frag)
		clrblock(&sblock, acg.cg_free, d/sblock.fs_frag);
	if (d > dupper) {
		acg.cg_frsum[d - dupper]++;
		for (i = d - 1; i >= dupper; i--) {
			setbit(acg.cg_free, i);
			acg.cg_cs.cs_nffree++;
		}
	}
	while ((d + sblock.fs_frag) <= dmax - cbase) {
		setblock(&sblock, acg.cg_free, d/sblock.fs_frag);
		acg.cg_cs.cs_nbfree++;
		acg.cg_btot[cbtocylno(&sblock, d)]++;
		acg.cg_b[cbtocylno(&sblock, d)][cbtorpos(&sblock, d)]++;
		d += sblock.fs_frag;
	}
	if (d < dmax - cbase) {
		acg.cg_frsum[dmax - cbase - d]++;
		for (; d < dmax - cbase; d++) {
			setbit(acg.cg_free, d);
			acg.cg_cs.cs_nffree++;
		}
		for (; d % sblock.fs_frag != 0; d++)
			clrbit(acg.cg_free, d);
	}
	for (d /= sblock.fs_frag; d < MAXBPG(&sblock); d ++)
		clrblock(&sblock, acg.cg_free, d);
	sblock.fs_cstotal.cs_ndir += acg.cg_cs.cs_ndir;
	sblock.fs_cstotal.cs_nffree += acg.cg_cs.cs_nffree;
	sblock.fs_cstotal.cs_nbfree += acg.cg_cs.cs_nbfree;
	sblock.fs_cstotal.cs_nifree += acg.cg_cs.cs_nifree;
	*cs = acg.cg_cs;
	wtfs(fsbtodb(&sblock, cgtod(&sblock, cylno)),
		sblock.fs_bsize, (char *)&acg);
}

/*
 * initialize the file system
 */
struct inode node;
#define PREDEFDIR 3
struct direct root_dir[] = {
	{ ROOTINO, sizeof(struct direct), 1, "." },
	{ ROOTINO, sizeof(struct direct), 2, ".." },
	{ LOSTFOUNDINO, sizeof(struct direct), 10, "lost+found" },
};
struct direct lost_found_dir[] = {
	{ LOSTFOUNDINO, sizeof(struct direct), 1, "." },
	{ ROOTINO, sizeof(struct direct), 2, ".." },
	{ 0, DIRBLKSIZ, 0, 0 },
};
char buf[MAXBSIZE];

fsinit()
{
	int i;

	/*
	 * initialize the node
	 */
	node.i_atime = uutime;
	node.i_mtime = uutime;
	node.i_ctime = uutime;
	/*
	 * create the lost+found directory
	 */
	(void)makedir(lost_found_dir, 2);
	for (i = DIRBLKSIZ; i < sblock.fs_bsize; i += DIRBLKSIZ)
		bcopy(&lost_found_dir[2], &buf[i], DIRSIZ(&lost_found_dir[2]));
	node.i_number = LOSTFOUNDINO;
	node.i_mode = IFDIR | UMASK;
	node.i_nlink = 2;
	node.i_size = sblock.fs_bsize;
	node.i_db[0] = alloc(node.i_size, node.i_mode);
	node.i_blocks = btodb(fragroundup(&sblock, node.i_size));
	wtfs(fsbtodb(&sblock, node.i_db[0]), node.i_size, buf);
	iput(&node);
	/*
	 * create the root directory
	 */
	node.i_number = ROOTINO;
	node.i_mode = IFDIR | UMASK;
	node.i_nlink = PREDEFDIR;
	node.i_size = makedir(root_dir, PREDEFDIR);
	node.i_db[0] = alloc(sblock.fs_fsize, node.i_mode);
	node.i_blocks = btodb(fragroundup(&sblock, node.i_size));
	wtfs(fsbtodb(&sblock, node.i_db[0]), sblock.fs_fsize, buf);
	iput(&node);
}

/*
 * construct a set of directory entries in "buf".
 * return size of directory.
 */
makedir(protodir, entries)
	register struct direct *protodir;
	int entries;
{
	char *cp;
	int i, spcleft;

	spcleft = DIRBLKSIZ;
	for (cp = buf, i = 0; i < entries - 1; i++) {
		protodir[i].d_reclen = DIRSIZ(&protodir[i]);
		bcopy(&protodir[i], cp, protodir[i].d_reclen);
		cp += protodir[i].d_reclen;
		spcleft -= protodir[i].d_reclen;
	}
	protodir[i].d_reclen = spcleft;
	bcopy(&protodir[i], cp, DIRSIZ(&protodir[i]));
	cp += DIRSIZ(&protodir[i]);
	return (cp - buf);
}

/*
 * allocate a block or frag
 */
daddr_t
alloc(size, mode)
	int size;
	int mode;
{
	int i, frag;
	daddr_t d;

	rdfs(fsbtodb(&sblock, cgtod(&sblock, 0)), sblock.fs_cgsize,
	    (char *)&acg);
	if (acg.cg_magic != CG_MAGIC) {
		printf("cg 0: bad magic number\n");
		return (0);
	}
	if (acg.cg_cs.cs_nbfree == 0) {
		printf("first cylinder group ran out of space\n");
		return (0);
	}
	for (d = 0; d < acg.cg_ndblk; d += sblock.fs_frag)
		if (isblock(&sblock, acg.cg_free, d / sblock.fs_frag))
			goto goth;
	printf("internal error: can't find block in cyl 0\n");
	return (0);
goth:
	clrblock(&sblock, acg.cg_free, d / sblock.fs_frag);
	acg.cg_cs.cs_nbfree--;
	sblock.fs_cstotal.cs_nbfree--;
	fscs[0].cs_nbfree--;
	if (mode & IFDIR) {
		acg.cg_cs.cs_ndir++;
		sblock.fs_cstotal.cs_ndir++;
		fscs[0].cs_ndir++;
	}
	acg.cg_btot[cbtocylno(&sblock, d)]--;
	acg.cg_b[cbtocylno(&sblock, d)][cbtorpos(&sblock, d)]--;
	if (size != sblock.fs_bsize) {
		frag = howmany(size, sblock.fs_fsize);
		fscs[0].cs_nffree += sblock.fs_frag - frag;
		sblock.fs_cstotal.cs_nffree += sblock.fs_frag - frag;
		acg.cg_cs.cs_nffree += sblock.fs_frag - frag;
		acg.cg_frsum[sblock.fs_frag - frag]++;
		for (i = frag; i < sblock.fs_frag; i++)
			setbit(acg.cg_free, d + i);
	}
	wtfs(fsbtodb(&sblock, cgtod(&sblock, 0)), sblock.fs_cgsize,
	    (char *)&acg);
	return (d);
}

/*
 * Allocate an inode on the disk
 */
iput(ip)
	register struct inode *ip;
{
	struct dinode buf[MAXINOPB];
	daddr_t d;
	int c;

	c = itog(&sblock, ip->i_number);
	rdfs(fsbtodb(&sblock, cgtod(&sblock, 0)), sblock.fs_cgsize,
	    (char *)&acg);
	if (acg.cg_magic != CG_MAGIC) {
		printf("cg 0: bad magic number\n");
		exit(1);
	}
	acg.cg_cs.cs_nifree--;
	setbit(acg.cg_iused, ip->i_number);
	wtfs(fsbtodb(&sblock, cgtod(&sblock, 0)), sblock.fs_cgsize,
	    (char *)&acg);
	sblock.fs_cstotal.cs_nifree--;
	fscs[0].cs_nifree--;
	if(ip->i_number >= sblock.fs_ipg * sblock.fs_ncg) {
		printf("fsinit: inode value out of range (%d).\n",
		    ip->i_number);
		exit(1);
	}
	d = fsbtodb(&sblock, itod(&sblock, ip->i_number));
	rdfs(d, sblock.fs_bsize, buf);
	buf[itoo(&sblock, ip->i_number)].di_ic = ip->i_ic;
	wtfs(d, sblock.fs_bsize, buf);
}

/*
 * read a block from the file system
 */
rdfs(bno, size, bf)
	daddr_t bno;
	int size;
	char *bf;
{
	int n;

	if (lseek(fsi, bno * DEV_BSIZE, 0) < 0) {
		printf("seek error: %ld\n", bno);
		perror("rdfs");
		exit(1);
	}
	n = read(fsi, bf, size);
	if(n != size) {
		printf("read error: %ld\n", bno);
		perror("rdfs");
		exit(1);
	}
}

/*
 * write a block to the file system
 */
wtfs(bno, size, bf)
	daddr_t bno;
	int size;
	char *bf;
{
	int n;

	if (lseek(fso, bno * DEV_BSIZE, 0) < 0) {
		printf("seek error: %ld\n", bno);
		perror("wtfs");
		exit(1);
	}
	n = write(fso, bf, size);
	if(n != size) {
		printf("write error: 0x%x\n", bno);
		perror("wtfs");
		exit(1);
	}
}

/*
 * check if a block is available
 */
isblock(fs, cp, h)
	struct fs *fs;
	unsigned char *cp;
	int h;
{
	unsigned char mask;

	switch (fs->fs_frag) {
	case 8:
		return (cp[h] == 0xff);
	case 4:
		mask = 0x0f << ((h & 0x1) << 2);
		return ((cp[h >> 1] & mask) == mask);
	case 2:
		mask = 0x03 << ((h & 0x3) << 1);
		return ((cp[h >> 2] & mask) == mask);
	case 1:
		mask = 0x01 << (h & 0x7);
		return ((cp[h >> 3] & mask) == mask);
	default:
#ifdef	STANDALONE
		printf("isblock bad fs_frag %d\n", fs->fs_frag);
#else	STANDALONE
		fprintf(stderr, "isblock bad fs_frag %d\n", fs->fs_frag);
#endif	STANDALONE
		return;
	}
}

/*
 * take a block out of the map
 */
clrblock(fs, cp, h)
	struct fs *fs;
	unsigned char *cp;
	int h;
{
	switch ((fs)->fs_frag) {
	case 8:
		cp[h] = 0;
		return;
	case 4:
		cp[h >> 1] &= ~(0x0f << ((h & 0x1) << 2));
		return;
	case 2:
		cp[h >> 2] &= ~(0x03 << ((h & 0x3) << 1));
		return;
	case 1:
		cp[h >> 3] &= ~(0x01 << (h & 0x7));
		return;
	default:
#ifdef	STANDALONE
		printf("clrblock bad fs_frag %d\n", fs->fs_frag);
#else	STANDALONE
		fprintf(stderr, "clrblock bad fs_frag %d\n", fs->fs_frag);
#endif	STANDALONE
		return;
	}
}

/*
 * put a block into the map
 */
setblock(fs, cp, h)
	struct fs *fs;
	unsigned char *cp;
	int h;
{
	switch (fs->fs_frag) {
	case 8:
		cp[h] = 0xff;
		return;
	case 4:
		cp[h >> 1] |= (0x0f << ((h & 0x1) << 2));
		return;
	case 2:
		cp[h >> 2] |= (0x03 << ((h & 0x3) << 1));
		return;
	case 1:
		cp[h >> 3] |= (0x01 << (h & 0x7));
		return;
	default:
#ifdef	STANDALONE
		printf("setblock bad fs_frag %d\n", fs->fs_frag);
#else	STANDALONE
		fprintf(stderr, "setblock bad fs_frag %d\n", fs->fs_frag);
#endif	STANDALONE
		return;
	}
}

#ifdef STANDALONE

dump_disk(diskmodel)
	register model_t *diskmodel;
{
	register int i;

	printf ("\n%d Cylinders %d Heads %d Sectors Per Track\n",
		diskmodel->m_nCyl, diskmodel->m_nHead, diskmodel->m_nSectors);
#ifdef notdef
	printf ("Gaps - 1 (%d), 2 (%d), 3 (%d) {Not for SMD}\n",
		diskmodel->m_Gap1, diskmodel->m_Gap2, diskmodel->m_Gap3);
	printf ("%s Interface, Disk Type %d, Skew %d ",
		(diskmodel->m_Interface == smd) ? "SMD" :
		(diskmodel->m_Interface == esdi) ? "ESDI" :
		"ST506", diskmodel->m_DiskType, diskmodel->m_Skew);
	printf ("Grouping %d, Head Load Time %d, Seek Time %d\n",
		diskmodel->m_Grouping, diskmodel->m_HeadLoad,
		diskmodel->m_SeekTime);
	printf ("Cache %s, Zero Latency %s Buffered Step",
		diskmodel->m_Cache?"ON":"OFF",
		diskmodel->m_ZeroLat?"ON":"OFF");
	if (diskmodel->m_BufStep)
		printf (" ON, Width %d, Interval %d\n",
			diskmodel->m_StepWidth, diskmodel->m_StepInt);
	else
		printf (" OFF\n");
	printf ("Write Precomp Cyl %d, Reduced Write Cyl %d\n",
		diskmodel->m_Precomp, diskmodel->m_RedWrite);
#endif notdef
	printf ("Root is %d Cyl, Small Swap is %d Cyl, Large Swap is %d Cyl\n",
		diskmodel->m_Root, diskmodel->m_Swap, diskmodel->m_LgSwap);
	printf ("Small User is %d Cyl, Large User is %d Cyl\n",
		diskmodel->m_Root, diskmodel->m_User, diskmodel->m_LgUser);
#ifndef FORMAT_FROM_DKMODELS
	printf("Part.\tCyloff\tNblk\n");
	for (i=0; i<DK_NSIZES; i++)
		printf("%x\t%x\t%x\n",i,diskmodel->m_PartSizes[i].sz_cyloff,
			diskmodel->m_PartSizes[i].sz_nblk);
#endif  FORMAT_FROM_DKMODELS
	printf("Frag Size is 0x%x, F.S. Block Size is 0x%x\n",
		diskmodel->m_FragSize, diskmodel->m_FSBlockSize);

}

dump_dkname(diskmodel)
	model_t *diskmodel;
{
	register struct name_list *np;

	printf ("Model ");
	np = diskmodel->m_name;
	if (!np) printf ("UNKNOWN");
	else while (np) {
		printf ("%s", np->nl_name);
		np = np->nl_next;
		if (np) printf (", ");
	}
	printf("\n");
}

dump_dkmodels ()
{
	register model_t *mop;

	mop = dk_models;
	while (mop)
	{
		dump_dkname(mop);
		/* dump_disk(mop); /* Blab, blab, blab */
		mop = mop->m_next;
	}
	
}
#endif STANDALONE
